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Periodic  Precipitation  of  Liesegang  System 
under  Microwave  Radiation 
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Abstract:  We  studied  precipitation  patterns  in  a  Liesegang  system  under  MW  (microwave)  irradiation  in  order  to  investigate  metal 
salt  diffusion  in  an  electrolyte  gel.  The  gel  and  salt  concentrations  were  varied.  MW  irradiation  induced  periodic  patterns  of 
precipitation  because  polar  molecules  vibrate  and  rotate  in  an  electromagnetic  field.  For  example,  the  number  of  patterns  increased 
by  the  irradiation.  Accordingly,  microwave  irradiation  nonlinearly  accelerated  the  diffusion  of  ionic  molecules. 

Key  words:  MW,  Liesegang,  diffusion,  precipitation. 


1.  Introduction 

Liesegang  phenomena,  a  function  of  diffusion  and 
reaction  conditions,  form  attractive  and  unique 
periodic  stripe  patterns.  Such  pattern  formation 
frequently  inspires  microstructural  design  in  materials 
science.  Moreover,  morphological  pattern  transitions 
between  periodic  precipitation  and  tree-like  crystal 
aggregates  are  scientifically  interesting  [1,  2]. 
Liesegang  phenomena  usually  manifest  themselves 
when  a  concentrated  salt  solution  diffuses  through  a 
gel  that  contains  another  electrolyte  at  a  lower 
concentration.  Diffusion  and  chemical  reactions  cause 
little  precipitation  of  otherwise  soluble  salts  in  a 
periodic  manner.  The  precipitation  and  redissolution 
reactions  for  a  generic  metal  ion  M2+  are  as  follows: 
M2+  (aq)  +  2NH4OH  (aq)  ^  2NH4+  (aq)  +  M(OH)2(s)       ( 1 ) 

M(OH)2  (s)  +  6NH4+  (aq)  -> 
M(NH3)62+  (aq)  +  2H20  (1)  +  4H+  (aq) 
When  NH4OH  diffusion  causes  metal  ions  to 
precipitate  in  the  form  of  M(OH)2,  redissolution  at  the 
top  of  the  pattern  induces  a  Liesegang  band. 

Recent    research    indicates     that  spontaneous 


(2) 
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diffusion-related  self-organization  dynamics  and 
irreversible  reactions  or  precipitation  systems  based 
on  Ostwald  supersaturation  theory  explain  Liesegang 
pattern  formation  [3].  However,  nonlinearity  hinders 
elaboration  of  the  exact  mechanism. 

MW  (microwave)  radiation  may  be  helpful  in  this 
context.  MW  radiation  has  recently  been  a  very  active 
area  of  research  due  to  the  many  advantages  of 
nonthermal  effects  such  as  reaction  promotion, 
nucleation  induction,  and  diffusion  facilitated  by  the 
molecular  vibration  of  polar  molecules  [4,  5].  We 
hypothesized  that  the  reaction  and  diffusion 
characteristics  of  a  Liesegang  system  could  be 
modulated  or  explained  via  various  MW  properties.  In 
addition,  Liesegang  phenomena  may  suggest  a  new 
application  of  MW  radiation  and  suitable  operating 
conditions  for  highly  functional  inorganic 
nanomaterials.  Salt  and  gel  concentration  are  essential 
for  solving  the  diffusion  challenge  of  Liesegang 
phenomena.  In  this  study,  we  varied  gel  and  salt 
concentrations  to  understand  the  MW  effect.  We  used 
a  MW  reactor  equipped  with  an  in  situ  observation 
system  to  analyze  the  dynamic  growth  behavior  of  the 
periodic  pattern. 
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2.  Experiments 

To  form  a  Liesegang  pattern,  we  used  a  gel  because 
it  can  prevent  precipitate  sedimentation  and 
hydrodynamic  turbulence.  A  salt  was  added  to  a 
prepared  gelatin-based  gel  [6,  7].  In  this  study,  gelatin 
and  salt  concentrations,  Cg,  Cc  were  varied,  and  we 
used  Co(NC>3)2  (Table  1).  We  first  heated  a  solution 
(i.e.,  gel  and  salt)  to  90  °C.  We  then  poured  the 
solution  into  a  glass  test  tube  (<?:  13  mm;  height:  100  mm; 
volume:  5  mL)  and  gradually  cooled  it  to  room 
temperature  to  impart  a  homogeneous  gelatin  network 
structure.  Temperature  is  the  most  important  factor  for 
pattern  formation;  we  used  two  different  methods, 
MW  irradiation  and  conventional  heating  (OB,  oil 
bath)  for  maintaining  solution  temperature.  After  we 
poured  2  mL  of  a  28%  ammonium  solution  at  the  gel 
surface,  we  covered  it  with  a  silicon  rubber  lid. 
Ammonium  salt  precipitates  via  ammonium  ion 
diffusion  into  the  gel,  and  Liesegang  stripe-pattern 
formation  starts  from  the  top  surface  to  the  bottom  of 
the  tube.  We  measured  the  stripe  number  and  length 
for  the  two  heating  methods. 

2.1  MW  Heating 

Fig.  1  shows  a  tube-guide  MW  reactor  (Shikoku 
Instrumentation),  which  is  specially  designed  for 
preventing  MW  leakage  and  equipped  with  a 
microscope  camera  (Sigma  Koki,  model 
SK-TC202US B -  AT)  for  in  situ  observation.  We  hung 
a  test  tube  containing  the  solution  by  a  thread  at  the 


Camera 


■a 


Sample 


center  of  the  reactor.  Immediately  after  pouring  the 
ammonium  solution  onto  the  gel  surface,  we  initiated 
MW  irradiation  for  30  min.  Because  the  stripe  pattern 
grows  from  the  surface  to  the  bottom,  we  gradually 
lifted  up  the  test  tube  at  a  controlled  speed  to  acquire 
photographs  around  the  vertical  stripe  center,  through 
the  side  of  the  reactor,  with  the  fixed  microscope.  We 
measured  the  Liesegang  pattern  length  via  the  pixel 
number  between  the  top  and  bottom  of  the 
precipitation  region,  and  we  counted  the  stripes 
through  the  photographs  every  5  min.  This  approach 
enables  dynamic  evaluation  of  Liesegang  pattern 
propagation  with  respect  to  MW  radiation.  After  30 
min,  we  measured  the  final  temperature  with  an 
optical  fiber  (Anritsu  Meter)  because  the  solution 
temperature  increases  under  irradiation. 

2.2  Conventional  Heating 

We  used  an  OB  to  maintain  the  solution 
temperature.  The  bath  temperature  was  increased  over 
the  course  of  30  min  to  the  final  temperature  observed 
in  the  MW  experiment.  We  poured  ammonium 
solution  into  the  test  tube  immersed  in  the  OB,  and 
photographed  the  stripe  pattern  after  30  min  with  the 
microscope  camera. 

Table  1    Experimental  conditions. 


Gel  concentrate 

Cg  (wt%) 

Salt  concentrate 
Cc  (mol/L) 

MW  output  (W) 
or  temperature 

(°C) 

MW    0.5,  1,  2 

0.3,  0.6,  0.9 

25  (W) 

OB      0.5,  1,  2 

0.3,  0.6,  0.9 

40  (°C) 

(a)  Side  view 

Fig.  1    MW  reactor. 


(b)  3D  view 
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3.  Results  and  Discussion 

Figs.  2  and  3  show  photographs  of  the  stripe  pattern 
near  the  gel  surface  for  each  gel  concentration  under 
MW  or  OB  conditions  after  30  min.  The  experiment 
involves  precipitation  after  the  salt  concentration  in 
the  gel  exceeds  the  critical  values  (i.e.,  the  solubility 
due  to  ammonium  ion  diffusion).  Although  a 
stripe  pattern  was  observed,  an  expected  pattern 
transition  (e.g.,  a  tree-like  pattern  [1])  was  not 
observed  in  this  experiment.  On  the  other  hand,  the 
parallel  stripe  pattern  remained  even  under  MW 
irradiation. 

To  aid  in  understanding  the  mechanism  of 
stripe-pattern  formation,  Figs.  4  and  5  plot  the  stripe 
number  and  length,  respectively,  in  relation  to 
reaction  time.  Filled  and  open  symbols  indicate  MW 
and  OB  conditions,  respectively.  Although 
precipitation  starts  just  after  contact  with  the 
ammonium  solution,  either  the  stripes  do  not 
immediately  form  or  the  camera  resolution  limit 
precludes  sufficiently  fast  and  precise  image 
acquisition;  accordingly,  we  assigned  the  initial 
number  of  stripes  as  0.  We  counted  the  number  of 
stripes  after  approximately  15-25  min.  The  stripe 
length  linearly  increased  with  reaction  time  (Fig.  5). 
After  30  min  of  reaction  time,  nonthermal  effects 
induced  by  MW  irradiation  caused  the  stripe  number 
and  length  to  increase  more  rapidly  for  MW  relative 
to  OB.  This  is  attributable  to  enhanced  molecular 
diffusion  and  facilitated  supersaturation.  Molecular 
vibration  induced  by  MW  irradiation  accelerates  ionic 
diffusion.  Moreover,  such  vibration  readily  induces 
supersaturation  of  the  precipitate.  Consequently, 
cyclic  periods  of  supersaturation,  nucleation,  and 
depletion  become  faster.  However,  temperature  must 
be  carefully  controlled  to  modulate  such  cycling. 

Figs.  6  and  7  show  an  almost  linear  relation  of 
stripe  number  and  length,  after  30  min  of  irradiation, 
with  respect  to  Co(N03)2  concentration.  MW 
imparted  a  greater  effect  than  OB  in  case  of  the 


lower  salt  concentration.  Consequently,  MW 
irradiation  more  effectively  induced  precipitation  via 
supersaturation. 


(a)  Cg  =  0.5  wt%  (b)  Cg  =  1.0  wt%  (c)  Cg  =  2.0  wt% 
Fig.  2    Stripe  pattern  under  MW  at  Cc  =  0.3  mol/L. 


(a)  Cg  =  0.5  wt%     (b)  Cg  =  1 .0  wt%    (c)  Cg  =  2.0  wt% 


Fig.  3    Stripe  pattern  in  OB  at  Cc  =  0.3  mol/L. 
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Fig.  4    Stripe  number  for  the  time  at  Cc  =  0.9  mol/L. 
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5    Pattern  length  for  the  time  at  Cc  =  0.9  mol/L. 
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Fig.  6  Stripe  number  for  the  cobalt  concentration  after  30  min. 
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Cobalt  nitrate  concentration  (mol/L) 

Fig.  7   Pattern  length  for  cobalt  concentration  after  30  min. 

4.  Conclusions 

To  understand  the  mechanisms  of  precipitation  and 
diffusion  in  Liesegang  patterns,  experiments  varying 
the  salt  and  gel  concentrations  were  conducted  in  an 
oil  bath  and  MW  reactor  equipped  with  in  situ 
observation  system.  Stripe  length  and  number 
increased  more  quickly  with  respect  to  time  under 


MW  irradiation  than  oil  bath  heating,  and  also 
decreased  with  respect  to  salt  concentration.  In 
particular,  stronger  MW  effect  was  confirmed  at  lower 
salt  condition.  We  explained  these  findings  in  the 
context  of  nonthermal  effects  (i.e.,  higher  molecular 
diffusion  and  facilitated  supersaturation)  that  are 
accelerated  by  radiation-induced  molecular  vibration. 
MW  irradiation  can  modulate  the  diffusion  and 
manifestation  of  Liesegang  patterns,  and  may  be  a 
useful  tool  in  other  fine-scale,  nonequilibrium 
self-organization  processes  that  involve  precipitation 
and  nonlinearity. 
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Abstract:  Sodium  bis[N-2-oxyphenyl-5-substituted-salicylideneiminato-ONO]ruthenate(III)  hemitriethylamine  compounds  of 
general  formula  Na[Ru(N-R-5-X-salim)2]-0.5Et3N,  where,  R  =  C6H40,  X  =  CI,  Br,  N02,  salim  =  salicylideneiminato,  Et3N 
triethylamine,  were  synthesized  and  characterized  on  the  basis  of  elemental  analysis,  MALDI-TOF  (matrix-assisted  laser 
desorption/ionization/time-of-flight)  mass  spectra,  'H  NMR  (nuclear  magnetic  resonance),  electronic  spectra  and  cyclic 
voltammograms.  Ru(III)  is  chelated  with  two  02N  anionic  tridentate  Schiff  bases  derived  from  5-X-salicyladehyde  and 
2-aminophenol.  Spectroscopic  titrations  of  the  Ru(III)  compounds  with  CT  DNA  (calf  thymus  DNA)  in  the  LMCT  (ligand  metal 
charge  transfer)  region  have  shown  moderate  intercalation  properties  of  compounds  with  binding  constants  Kb  =  2.06-3.85  x  104  M"  . 
Electrochemical  evidence  of  intercalative  mode  of  binding  is  based  on  the  decrease  of  current  and  shift  of  anodic  and  cathodic  peaks 
towards  higher  values  of  potential  with  uE- increment  of  CT  DNA  in  Ru(III)  compound  solution. 

Keywords:  Ruthenium,  Schiff  bases,  salicylideneimine,  CT  DNA,  intercalation. 


1.  Introduction 

Ruthenium  complexes  are  the  subject  of  increasing 
interest  in  chemistry  and  interdisciplinary  fields  of 
science  due  to  many  significant  properties,  especially 
anticancer  and  antimetastatic  properties  [1-8]  and  the 
ability  to  catalyze  a  number  of  organic  reactions  [9-12] 
and  the  aging. 

A  number  of  ruthenium  complexes  showed  electron 
transfer  capability  as  a  key  feature  in  the  development 
of  new  sensors  [13-14],  and  appropriate  thermal  and 
electrical  properties  essential  for  applications  in  solar 
cells  [15].  In  the  development  of  new  compounds,  as 
potential  antitumor  agents,  DNA  is  still  considered  as 
the  primary  target.  Metal  compounds  are  known  to 
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bind  DNA  by  covalent  or  non-covalent  mode.  For 
covalent,  either  inter-  or  intra-strand  mode  of  binding, 
easily  leaving  groups,  e.g.,  chlorides,  are  generally 
required.  Non-covalent  binding  mainly  assumes 
groove  binding  or  intercalation  which  means  insertion 
of  fused  and  planar  aromatic  moiety  between  base 
pairs  of  DNA.  Many  metal  complex  species,  in 
particular  neutral  or  positively  charged,  with  some 
typical  intercalating  organic  molecules  containing 
extended  7i-electronic  system,  e.g.,  2,2-bipyridine 
(bipy)  or  1 , 1 0-phenanthroline  (phen)  have  shown 
ability  to  intercalate  DNA  [16-21].  Some  Ru(III) 
complexes  with  such  ligands  have  also  shown  the 
property  to  bind  DNA  by  intercalation  [22-24]. 

Schiff  bases  derived  from  salicyladehyde,  either  as 
neutral  or  anionic  species,  are  considered  to  be  very 
functional  ligands,  therefore  many  metal  complexes 
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with  this  type  of  ligands  were  reasonably  synthesized 
for  different  purposes.  Some  of  ruthenium  complexes 
which  contain  ON  (oxygen,  nitrogen)  or  ONO 
(oxygen,  nitrogen,  oxygen)  anionic  salicylidenimine 
as  ligands  showed  antifungal  and  antibacterial  activity 
[25-27].  Recently,  we  reported  the  DNA  binding 
properties  of  two  anionic  Ru(III)  compounds 
containing  anionic  bidentate  salicylidenimine  ligands 
and  two  chlorides  [28]. 

In  continuation  of  our  study  on  synthesis  and 
characterization  of  Ru(III)  compounds  with  Schiff 
bases  and  interaction  with  DNA,  we  report  herein  on 
Ru(III)  compounds  containing  two  anionic  O2N  Schiff 
bases  derived  from  5-substituted-salicylaldehyde  and 
2-aminophenol  and  study  of  interaction  with  CT  DNA 
(calf  thymus  DNA). 

2.  Experiments 

All  chemicals  of  analytical  grade  were  purchased 
from  Sigma  Aldrich  and  Merck  and  were  used  without 
further  purification.  Calf  thymus  DNA  type  I-fibrous 
(^260/^280  =  1.3)  was  purified  with  phenol- 
chloroform-isoamyl  alcohol  procedure  until  ratio 
^26o/^280  =  1-8  was  obtained  and  then  DNA  was 
precipitated  as  sodium  salt  by  adding  0.1  M  acetate 
buffer  (pH  4.60)  and  ethanol.  The  precipitate  was 
dried  at  50  °C  and  kept  at  4  °C  for  further  use.  Solid 
nucleic  acid  of  sufficient  purity  was  suspended  in  0. 1  M 
Tris-HCl  buffer,  pH  7.40,  and  left  one  day  for 
hydration  at  4  °C.  The  stock  solution  of  CT  DNA  was 
prepared  in  Tris-HCl  buffer  at  pH  7.40  and  stored  at  4  °C 
for  a  maximum  of  1-4  days.  The  concentration  of  CT 
DNA  solution  was  calculated  on  the  basis  of  DNA 
extinction  coefficient  s  =  6,600  IVr'-cm"1. 

2. 1  Synthesis 

2.1.1  General  Procedure  for  Synthesis  of 
Na[Ru(N-R-5-X-salim)2]-0.5  Et3N 

The  solution  of  0.76  mmol  of  appropriate 
N-2-hydroxyphenyl-5-X-salicylideneimine  (X  =  CI, 
Br,  NO2)  in  50  mL  absolute  ethanol  was  mixed  with 


0.76  mmol  triethylamine,  afterwards  0.38  mmol 
RuCl3-3H20  in  5  mL  absolute  ethanol  was  added.  The 
solution  was  heated  in  a  rotary  evaporator  at  60-70  °C 
during  4  h  and  then  was  cooled  to  room  temperature. 
Small  volume  of  water  solution  (0.5  mL)  containig 
0.38  mmol  NaCl  was  added  and  resulting  solution  was 
left  over  night  at  ambient  temperature.  The  precipitate 
was  filtered  off,  washed  with  cold  water,  ethanol  and 
petrol  ether.  Recrystallization  was  performed  from 
absolute  ethanol/dichlormetan  (1/1,  v/v).  The 
compounds  were  dried  at  80  °C.  Yield:  60%-71%. 

2.1.2  Na[Ru(N-R-5-Cl-salim)2]0.5Et3N 

Anal.  Calc.  for  NaRuC26Hi6Cl2N2O40.5C6H15N 
(%):  C  52.30,  H  3.56,  N  5.26.  Found  (%):  C  52.32,  H 
4.15,  N  5.55.  MALDI-TOF/TOF  MS  mlz:  591.9534 
([C26H16N204Cl2Ru]-,  100%).  IR  (KBr)  vmJcmA: 
1,598  vs  (C=N),  1,287  m  (C-Ophe„),  1,033  w  (C2N).  IR 
(Csl)  vjcm"1:  460  w  (Ru-N),  387  w  (Ru-O).  UV/Vis 
[CH2Cl2,/Ux/nm]:  344,  395,  580  (e/M^-cm1):  12,394, 
14,457,  5,320.  'H  NMR  (nuclear  magnetic  resonance) 
Na  [Ru(/V-R-5-Cl-salim)2]-0.5  C6H15N  (600  MHz, 
acetone-d6)  5  10.23  s  (4  Hd),  7.94  dd  (4  Hc,  J=  8.0  Hz, 
0.5  Hz),  7.82  td  (4  Hg,  J=  6.9  Hz,  1.3  Hz),  7.54  dd  (4 
He,  J=  8.8  Hz,  2.8  Hz),  7.34  dd  (4  Ha,  J=  7.9  Hz,  0.8 
Hz),  7.12  dd  (4  Hh,  J=  9.0  Hz,  2.7  Hz),  6.88  td  (4  Hf, 
J  =  7.5  Hz,  4.1  Hz),  3.10q(6Hj,/=4.7Hz),  1.18  t  (9 

H,  /=7.3  Hz). 

2.1.3  Na[Ru(N-R-5-Br-salim)2]-0.5Et3N 

Anal.  Calc.  for  NaRuC26Hi6Br2N2O4-0.5  C6Hi5N 
(%):  C  46.14,  H  3.14,  N  4.64.  Found  (%):  C  47.46,  H 
3.84,  N  5.04.  MALDI-TOF/TOF  MS  mlz:  681.8532 
([C26H16N204  Br2Ru]\  (100%).  IR  (KBr)  v^/cm"1: 

I,  598  vs.  (C=N),  1,286  m  (C-Ophen),  1,032  w  (C2N). 
IR  (Csl)  vmi/cm"':  469  w  (Ru-N),  393  w  (Ru-O). 
UV/Vis  [CH2CI2,  Unm]:  345,  391,  588  (e/M1  cm"1): 
13,060,  14,461,  5,150.  'h  NMR 
Na[Ru(/V-R-5-Br-salim)2]-0.5  C6H15N  (600  MHz, 
acetone-d6)  5  10.21  s  (4  Hd),  7.93  dd  (4  Hc,  J=  8.0  Hz, 
0.6  Hz),  7.82  td  (4  Hg,  J=  8.2  Hz,  1.1  Hz),  7.77  dd  (4 
Ha,  J=  2.5  Hz,  0.8  Hz),  7.65  dd  (4  He,  J=  8.8  Hz,  2.6 
Hz),  7.21  dd  (4  Hh,  J=  9.0  Hz,  2.5  Hz),  7.03  td  (4  Hf, 


Calf  Thymus  DNA  Intercalation  by  Anionic  Ru(lll)  Complexes  Containing  Tridentate  Schiff  Bases  337 
Derived  from  5-X-Substituted  Salicyladehyde  and  2-Aminophenol 


J  =  7.9  Hz,  0.9  Hz),  3.09  q  (6Hj,  J=  A  A  Hz),  1.18  t  (9 
H,  7=7.3  Hz). 

2. 1 .4  Na[Ru(N-R-5-NO2-salim)2]-0.5Et3N 
Anal.  Calc.  for  NaRuC26H16N4O8-0.5  C6Hi5N  (%): 
C  50.69,  H  3.45,  N  9.17.  Found  (%):  C  49.58,  H  4.05, 
N  8.63.  MALDI-TOF/TOF  MS  mlz:  614.0018 
([C26H16N408Ru]",  100%).  IR  (KBr)  vmJcml:  1,595 
vs.  (C=N),  1,300  m  (C-0  phen),  1,033  w  (C2N).  IR  (Csl) 
vmi/cm"':  471  w  (Ru-N),  384  w  (Ru-O).  UV/Vis 
[CH2C12,  Umn]:  335,  371,  589  (e/M"1  cm"1):  17,273, 
18,941,  5,137.  'H  NMR  Na[Ru(/V-R-5-N02- 
salim)2]-0.5  C6H15N  (600  MHz,  acetone-d6)  5  10.28  s 
(4  Hd),  9.68  dd  (4  Hc,  J=  1.1  Hz,  0.24  Hz),  8.16  dd  (4 
Ha,  J=  2.3  Hz,  0.5  Hz),  8.33  dd  (4  Hb,  J=  9.1  Hz,  2.8 
Hz),  8.03  dd  (4  He,  J=  8.3  Hz,  0.84  Hz),  7.92  td  (4  Hg, 
J  =  9.3  Hz,  2.7  Hz),7.38  dd  (4  Hh,  J=  8.52  Hz,  0.82 
Hz),  7.06  td  (4  Hf,  J  =  8.4  Hz,  1.3  Hz),  3.10  q  (6Hj,  J 
=  1.2  Hz),  1.18  t  (9  H;, /=  7.3  Hz). 

2.2  Physical  Measurements 

Mass  spectra  were  recorded  on  a  complex 
matrix-assisted  laser  desorption/iomzation-time-of-flight 
MALDI-TOF/TOF  mass  spectrometer  (4800  Plus 
MALDI  TOF/TOF  Analyzer,  Applied  Biosystems  Inc., 
Foster  City,  CA,  USA)  equipped  with  Nd:  YAG  laser, 
wavelength  of  355  nm.  The  sample  was  resuspended 
in  10  |iL  THAP/DAC  MALDI  matrix  (1.9  mg 
2,4,6-Trihydroxyacetophenone  and  2.3  mg  ammonium 
hydrogen  citrate  in  1  mL  mixed  solvent 
acetonitrile/H20,  1/1,  v/v)  and  1  |iL  was  set  on  the 
MALDI  plate.  Mass  spectra  were  recorded  at  200  Hz 
in  a  negative  mode,  in  the  range  mlz  values  10-1,000 
Da  with  focus  on  the  mass  500  Da.  As  an  internal 
standard  in  the  negative  mode,  riboflavin  and 
3 -aminosalicylic  acid  were  used. 

'h  NMR  spectra  were  recorded  on  600  MHz  Bruker 
BioSpin  GmbH  instrument  in  acetone -d6  using  TMS 
(tetramethylsilane)  as  internal  standard. 

Infrared  spectra  were  collected  on  Perkin  Elmer  BX 
FT-IR  (fourier  transform  infrared)  spectrophotometer 
as  KBr  pallets  in  the  region  4,000-400  cm"1  and  for 


detection  Ru-N  and  Ru-0  absorptions  CsJ  pallets  were 
used  in  the  region  600-380  cm"1. 

Electronic  spectra  were  recorded  on  Perkin  Elemer 
UV/visible  lambda  35  spectrophotometer  in  CH2C12 
solution  over  the  range  200-700  nm. 

Cyclic  voltammograms  of  Na[Ru(N-R-5-X- 
salim)2]-0.5Et3N  were  recorded  on  an  electrochemical 
workstation  Autolab  potentiostat/galvanostat  (PGSTAT 
12)  using  three  electrode  system:  glassy  carbon 
electrode  as  working,  Ag/AgCl  as  reference  electrode 
and  Pt  wire.  We  used  acetonitrile  as  solvent  and 
tetraethylammonium  perchlorate  as  supporting 
electrolyte.  Electrochemical  titrations  of 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N  with  CT  DNA  were 
performed  on  the  same  instrument  with  the  same 
three-electrode  system  in  Tris-HCl  buffered  water 
solutions  in  potential  range  -1.0  V  to  -0.1  V.  The 
compounds  were  initially  dissolved  in  DMSO 
(dimethyl  sulfoxide),  afterwards  diluted  with  Tris-HCl 
buffer,  pH  =  7.40,  in  the  presence  of  150  mM  NaCl, 
up  to  a  total  volume  of  2  mL.  DNA  is  added  in  25  uL 
portions  and  cyclic  voltammograms  were  recorded 
after  each  addition  of  DNA. 

Spectroscopic  titration  of  Ru(III)  compounds  with 
CT  DNA  was  performed  by  the  successive  addition  of 
DNA  (six  portions  of  1 0  (iL)  in  2  mL  of  the  compound 
solutions  at  pH  7.40  in  the  presence  of  150  mM  NaCl. 
Spectra  were  recorded  2  min  after  addition  of  the 
DNA  in  the  range  200-700  nm.  Each  addition  of  DNA 
was  compensated  in  the  blank. 

3.  Results  and  Discussion 

3.1  Synthesis 

N-2-hydroxyphenyl-5-substituted-salicylideneimines, 
N-RH-5-X-salimH,  were  prepared  from  absolute 
ethanol  solutions  containing  equimolar  quantities  of 
2-aminophenol  and  5-X-salicyladehyde  where  X  =  CI, 
Br,  N02.  Reactions  were  performed  at  50  °C  during 
20  min.  The  orange  solids  were  purified  by 
recrystallization  from  warm  absolute  ethanol  with 
good  yields  of  80%-90%. 


RuCI3+2  [| 


OH     4  h,  70°C,  NaCI 


Fig.  1 
N02). 


L 

Synthesis  of  Na[Ru(N-R-5-X-salim)2],  (X  =  CI,  Br, 
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Sodium  bis[N-2-oxyphenyl-5-substituted- 
salicylideneiminato-ONO]ruthenate(III), 
Na[Ru(N-R-5-X-  salim)2]-0.5Et3N,  where,  X  =  CI,  Br, 
have  been  synthesized  according  to  published 
procedure  for  complex  anions  [29],  but  for  the 
purpose  of  interaction  with  DNA,  the  complex  anions 
were  precipitated  as  sodium  salts.  The  compound 
containing  NO2  substituted  salicylidenimine  was 
prepared  as  an  analogue.  The  starting  compounds 
RuCl3-3H20  and  N-RH-5-X-salimH  were  used  in 
molar  ratio  1:2.  Deprotonation  of  phenolic  oxygen 
atoms  from  ligands  was  performed  by  adding 
triethylamine  in  molar  ratio  1 :2,  respectively.  Reactions 
were  performed  at  60-70  °C  for  4  h  in  rotary  evaporator. 
The  precipitation  was  performed  with  stoichiometric 
amounts  of  NaCI  in  water  solution.  The  scheme  of 
synthesis  is  shown  in  Fig.  1.  The  solids  are  stable  in 
air,  insoluble  in  water,  soluble  in  acetonitrile,  N,N-DMF 
(dimethylformamide),  acetone,  ethyl  alcohol,  DMSO, 
dichlormethane. 

On  the  basis  on  CHN  (carbon  nitrogen  hydrogen) 
elemental  analysis,  mass  spectra,  'h  NMR,  infrared 
and  UV/visible  spectroscopic  measurements  the 
compounds  were  formulated  as 

Na[Ru(N-R-5-X-salim)2]-0.5Et3N.  Mass  spectra 
showed  molecular  ions  (M")  at  mlz  (100%)  at 
591.9534,  681.8532,  614.0018  for 

[RuC26H16Cl2N204]",       [RuC26H16Br2N204]"  and 
[RuC26Hi6N408]",  respectively. 

Ru(III)  in  octahedral  coordination  is  chelated  with 
two  02N  tridentate  Schiff  bases  trough  azomethine 
nitrogen  and  deprotonated  phenolic  oxygen  atoms. 
Elemental  analysis,  IR  and  !H  NMR  spectra  are  in 
accordance  with  presence  of  one  solvate  molecule  of 
triethylamine  per  two  complex  anions. 

3.2  Spectroscopic  Measurements 
3.2.1  !H  NMR  Spectra 

!H  NMR  spectra  of  Na[Ru(N-R-5-X-salim)2]-0.5Et3N 
showed  absence  of  singlet  in  region  10  ppm  to  12  ppm 


corresponding  to  phenolic  hydrogen,  consequently 
confirming  that  N-2-oxyphenyl-5-X-salicylideneimines 
bind  Ru(III)  via  deprotonated  phenolic  oxygen  atoms. 
Singlets,  positioned  in  this  region  (10.21  ppm  to  10.28 
ppm),  with  small  variations  in  the  chemical  shift  (0.09 
ppm),  correspond  to  azomethine  hydrogen  which  is 
slightly  sensitive  to  5-X  substituent  on  salicylaldehyde. 
The  positions  of  azomethine  singlet  correspond  to 
hydrogen  bonding  of  azomethine  group  with 
triethylamine.  Four  doublets  of  doublets  and  two 
triplets  of  doublets  in  region  7  ppm  to  9  ppm 
correspond  to  six  aromatic  hydrogen  atoms  with 
ortho-coupling  constants  Jhh  between  7.9  Hz  and  9.3 
Hz  and  metha-coupling  constants  Jhh  0.8  Hz  to  2.6  Hz. 
Quartet  at  3.1  ppm  and  triplet  at  1.18  ppm  confirm 
methylene  and  methyl  hydrogen  respectively  in  ratio 
1:1.5.  'h  NMR  data  strongly  support  proposed 
structures  of  complexes.  Different  hydrogen  atoms  in 
the  structure  of  [Ru(N-R-5-X-salim)2]"-0.5Et3N  are 
shown  in  Fig.  2. 

3.2.2  Infrared  Spectra 

Infrared  spectra  of  free  Schiff  bases  and 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N  are  in  accordance 
with  coordination  trough  azomethine  nitrogen  and 

H  H 

H-C-C-N(C2H6)2 
Hi  Hj 


Fig.  2  The  structure  of  [Ru(N-R-5-X-salim)2]"-0.5Et3N 
with  different  hydrogen  atoms  which  correspond  to  'H 
NMR  spectra. 
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both  phenolic  oxygen  atoms  in  Schiff  base.  The 
characteristic  frequencies  of  compounds  and  free 
ligands  are  given  in  Table  1.  Azomethine  group 
appears  in  spectra  of  free  N-RH-5-X-salimH  at 
1,630-1,629  cm"1  while  after  coordination  the 
frequency  of  azomethine  is  significantly  shifted 
toward  lower  values  for  32-29  cm"1  and  appears  in 
corresponding  Na[Ru(N-R-5-X-salim)2]-0.5Et3N 
compounds  at  1,598-1,595  cm"1.  Coordination  through 
deprotoned  phenolic  oxygen  caused  the  shift  of  C-0 
absorption  towards  higher  values  of  wave  numbers  for 
15-19  cm"1  as  a  result  of  the  formation  of  C-O(Ru) 
instead  of  C-O(H)  group.  Phenolic  C-O(H) 
frequencies  appear  at  1,272-1,281  cm"1  in  free  ligands, 
while  in  Na[Ru(N-R-5-X-salim)2]-0.5Et3N  were  found 
at  1,287-1,300  cm"1.  The  presence  of  solvate 
triethylamine  molecule  is  in  accordance  with  weak 
frequencies  corresponding  to  C2N  asymmetric 
stretching  at  1,033-1,032  cm"1  and  also  with  broad 
absorptions  centred  on  3,430  cm"1  which  can  be 
attributed  to  Et3N...H  (azomethine)  secondary 
interaction.  Weak  Ru-N  and  Ru-0  frequencies  in 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N  were  found  at 
460-471  cm"1  and  387-384  cm"1,  respectively. 
3.2.3  Electronic  Spectra 

Electronic  spectra  of  free  N-RH-5-X-salimH  and 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N,  recorded  in 
dichloromethane,  are  in  good  agreement  with  the 
mode  of  coordination  of  the  ligands.  Electronic 
spectra  of  free  Schiff  bases  have  three  or  four  in  the 
case  of  X  =  N02  compound,  well-defined  absorption 
bands  (Table  2).  Coordination  of  Schiff  bases 
significantly  changes  the  position  of  the  n(O)— >;:*  and 
n(N)— >7r*  bands  by  moving  them  toward  higher  and 
lower  wavelength  respectively  due  to  bonding  trough 
deprotonated  phenolic  oxygen  and  lone  pair  on 


azomethine  nitrogen.  Two  new  bands  were  found  in 
the  spectra  of  complex  compounds,  LMCT  bands  in 
the  region  371-395  nm  and  very  weak,  broad  and 
poorly  defined  absorptions  in  the  region  580-589  nm 
corresponding  to  t2g5^t2g4eg1  spin-allowed  electron 
transition  of  low-spin  Ru(III). 

3.3  Cyclic  Voltammetry 

Cyclic  voltammograms  of  Na[Ru(N-R-5-X- 
salim)2]-0.5Et3N  were  recorded  in  acetonitrile  in  the 
presence  of  tetraethylammonium  perchlorate  as 
supporting  electrolyte  using  glassy  carbon  working 
electrode  (Fig.  3). 

The  ratio  i/Jia  and  Ek-Ea  separations  in  cyclic 
voltammograms  indicate  quasi-reversible  one-electron 
processes.  Cathodic  peaks  appear  in  the  potential 
range  -883  mV  to  -915  mV,  anodic  peaks  at  -712  mV 
to  -742  mV  with  peaks  separations  141-183  mV. 
Formal  potentials  E1/2,  assigned  to  Ru(III)/Ru(II) 
couple,  are  quite  negative  varying  from  -804  mV  to 
-824  mV.  The  characteristic  potentials  are  given  in 
Table  3.  The  values  are  consistent  with  increased 
stabilization  of  Ru(III)  through  [Ru04N2]  skeleton 
compared  to  [Ru02N2]  in  the  compounds  with  similar 
type  of  ligand  [28]. 

3.4  Calf  Thymus  DNA  Binding 

3.4.1  Spectroscopic  Evidence  for  Intercalation 
Study  of  the  interaction  of  metal  compounds  with 
DNA  is  an  important  step  in  the  evaluation  of  the 
biological  properties  of  the  compounds.  Many  metal 
complexes  with  fused  and  planar  aromatic  ligands,  or 
aromatic  heterocyclic  rings,  are  described  as  DNA 
intercalators.  In  addition,  external  minor  or  major 
groove  binding  or  external  binding  to  the  phosphate 
backbone  is  possible.  Spectrophotometric  titration  of 


Table  1    IR  frequencies  for  Na[Ru(N-R-5-X-salim)2]-0.5Et3N  and  N-RH-5-X-salimH. 


X 

v(CH=N)  (cm"1) 

v(C-O)  (cm"1) 

v(Ru-N)  (cm"1) 

v(Ru-O)  (cm"1) 

CI 

1,598  (1,630) 

1,287  (1,272) 

460  (-) 

387  (-) 

Br 

1,598  (1,629) 

1,286  (1,270) 

469  (-) 

387  (-) 

NQ2 

1,595  (1,629) 

1,300  (1,281) 

471  (-) 

384  (-) 

Values  in  parentheses  refer  to  the  free  ligand  N-RH-5-X-salimH. 
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Table  2    Electronic  spectra  of  Na[Ru(N-R-5-X-salim)2]-0.5Et3N  and  N-RH-5-X-salimH. 


X 

X[n— >n*]  (nm) 

l[n(0)^>n*]  (nm) 

A[n(N)->jt*]  (nm) 

A[LMCT]  (nm) 

CI 

230  (230) 

296  (269) 

344  (366) 

395  (-) 

Br 

230  (230) 

297  (269) 

345  (366) 

391  (-) 

NQ2 

233  (229) 

295  (266) 

335  (347) 

371  (-) 

Values  in  parentheses  refer  to  the  free  ligand  N-RH-5-X-salimH. 


-1100 


Fig.        3  Cyclic        voltammograms  of 

Na[Ru(N-R-5-X-salim)2]-0.5Et3N  (X  =  CI,  Br,  N02)  in 
acetonitrile/tetraethyamnioniuni  perchlorate;  working 
electrode  GCE  vs.  Ag/AgCl;  scan  rate  300  mV-s"1. 


Table  3  The 

Na  [Ru(N-R-5-X-salim)2]  -0.5Et3N. 


potentials 


of 


X 

Ek(mV) 

Ea  (mV) 

Em  (mV) 

AEp  (mV) 

CI 

-895 

-712 

-804 

183 

Br 

-883 

-742 

-812 

141 

NQ2 

-915 

-734 

-824 

181 

metal  complex  with  DNA  is  a  useful  method  to 
estimate  binding  mode  and  calculate  binding  constant 
Kb.  In  this  work,  spectroscopic  study  of  interaction  of 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N  with  CT  DNA  was 
carried  out  by  titrating  a  fixed  amount  of  the  complex 
with  increasing  amount  of  calf  thymus  CT  DNA.  The 
[DNA]/[complex]  ratio  was  between  0-2.28  for  chloro, 
0-1.54  for  bromo-  and  0-2.30  for  NO2  derivative.  The 
binding  constants  Kb  were  calculated  on  the  basis  of 
Eq.  (1)  [30]: 

[DNA]  =  [DNA]    |  1 

(Ea-F.f)          (Eb-£f)  Kb(Ea-Ef) 

where,  sa,  £/  and  sh  represent  extinction  coefficients  for 
particular  measurements  (.4obs/[DNA]),  free  complex 
and  completely  bound  form,  respectively.  The  binding 
constant  Kb  is  obtained  as  the  ratio  of  the  slope  and 
Spectroscopic  titration  and  graphical  calculation  of 


intercept  in  graph  [DNA]/(efl  -  sj)  vs.  [DNA].  binding 
constant  are  shown  in  Figs.  4  and  5.  Apparent 
hypochromism  in  the  region  of  LMCT  absorptions  and 
values  of  Kb  (2.06-3.85  x  104  M"1)  designate 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N  as  moderate 
DNA-intercalators.  The  binding  constants  are  given  in 
Table  4.  The  impact  of  5-X-substituent  on  intercalative 

0.201  \ 


500  550 
X  (nm) 

Fig.  4  The  spectroscopic  titration  of  31.5  uM 
Na[Ru(N-R-5-Br-salim)2]-0.5Et3N  with  increasing 
concentration  of  CT  DNA  (0-49  uM)  in  Tris-HCl  buffer,  pH 
7.40  and  150  mM  NaCl. 


E 

O 

2 

to 

I 


< 

Q 


y  =  8.52  x  10  x 
3.92xl0-8x 

r  =  0.986 


0 


10°  [DNA]  (M) 
Fig.  5    Plot  of  (£/-  £„)  vs.  [DNA]/[DNA]. 
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Table       4  The       binding       constants  of 

Na[Ru(N-R-5-X-salim)2]-0.5Et3N  to  CT  DNA. 


X 

10*%  (M"1) 

CI 

2.06 

Br 

2.18 

NQ2 

3.85 

^E  (mV) 

-1000 


-250 


i  (10  [iA) . 


Fig.  6  Electrochemical  evidence  for  intercalation  29  uM 
Na[Ru(N-R-5-Cl-salim)2]-0.5Et3N  in  the  presence  of 
increasing  concentration  of  CT  DNA  (0-50  uM)  under 
phisiological  condition:  Tris-HCl,  pH  7.40,  150  mM  NaCl; 
working  electrode  GCE  vs.  Ag/AgCl;  scan  rate  800  mV-s"1; 
step  potential  50  mV. 

activities  of  the  compounds  correlates  with  electron 
withdrawing  properties  and  ability  of  substituent  to 
decrease  electron  density  from  ^-system  making  it  less 
nucleophilic. 

3.4.2  Electrochemical  Evidence  for  Intercalation 
Cyclic  voltammetry  is  also  a  useful  technique 
available  for  studying  the  interaction  of  metal 
complexes  with  DNA.  Procedure  is  carried  out  by 
successive  adding  of  uL-volumes  of  DNA  to  the 
solution  containing  fixed  amount  of 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N.  Intercalation  is 
followed  by  the  current  decrease  and  shift  of  cathodic 
and  anodic  peaks  towards  more  positive  values  of 
potentials,  though  the  peaks  are  poorly  visible  in 
aqueous  solution  due  to  insolubility  of  compounds  in 
water  (Fig.  6).  The  molar  ratio  r  [DNA]/[complex] 
was  within  the  range  0-1.8  for  X  =  CI,  and  0-1.5  for 
X=  Br,  N02. 

4.  Conclusions 

We    reported    here    on    synthesis    of  Ru(III) 


compounds  with  5-X  substituted  salicylideneimine 
ligands  derived  from  salicylaldehyde  and 
2-aminophenol.  The  compounds  of  the  general 
formula  Na[Ru(N-R-5-X-salim)2]-0.5Et3N  were 
synthesized  and  characterized  by  elemental  analysis, 
different  spectroscopic  techniques  and  cyclic 
voltammetry.  Ruthenium  is  chelated  by  two  tridentate 
Schiff  base  in  Ru04N2  octahedral  skeleton.  The 
interaction  with  CT  DNA  was  studied  by 
spectrophotometric  titration  and  cyclic  voltammetry 
indicating  capability  of  compounds  to  intercalate  CT 
DNA.  The  quantification  of  affinity  of  each  compound 
toward  DNA  was  done  by  binding  constants,  Kh  = 
2.06-3.85  x  104  M"1  which  correspond  to  intercalative 
mode  of  binding.  The  absorption  spectral  bands  at 
LMCT  region  showed  hypochromism  in  the  presence 
of  increasing  amount  of  DNA.  The  lack  of  measurable 
evidence  for  redshift  is  in  accordance  with  values  of 
binding  constants  and  moderate  intercalating 
properties  of  the  compounds.  Substituent  X  in 
Na[Ru(N-R-5-X-salim)2]-0.5Et3N  affects  an  increase 
of  binding  constants  with  increasing  electron 
withdrawing  properties  of  substituent  in  order  CI  ~  Br 
<N02. 
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Abstract:  At  pH  =  8,  the  complex  formation  process  of  Ni(II)  ions  with  FA  was  studied  by  the  solubility  method.  FA  were  separated 
from  the  river  Mtkvari  by  the  adsorption-chromotographic  method.  The  charcoal  (BAU,  Russia)  was  used  as  a  sorbent.  The  old 
suspension  of  Ni(II)  hydroxide  was  used  as  the  solid  phase,  on  which  was  added  the  increasing  quantity  of  standard  solution  of  FA. 
In  diluted  solutions,  at  pH  =  8.0  the  dominant  form  of  Ni(II)  is  nickel  dihydroxocomplex  Ni(OH)2°.  It  was  established  that  in  the 
Ni(OH)2(solid)-Ni(OH)20(sofution)-FA2"-H2O  system  dominates  nickel  dyhydroxofulvate  complex  with  the  structure  1:1, 
[Ni(OH)2FA]2~.  The  average  stability  constant  of  nickel  dyhydroxofulvate  complex  was  calculated  based  on  experimental  data 
/^i(OH)2FA2"  =  5.3xl05. 

Key  words:  Nickel,  FA,  average  stability  constant,  nickel  fulvate,  water. 

1.  Introduction 

FA  (fulvic  acids)  have  the  functional  groups,  so 
they  take  an  active  part  in  complexing  and  sorption 
processes,  proceeding  in  natural  waters,  suspended 
substances,  bottom  sediment  and  soils.  They  form 
stable  complexes  with  heavy  metals  and  radionuclides 
and  stipulate  migration  forms  thereof  in  natural  waters 
and  soils  [1-12]. 

The  carboxylic  and  phenolic  groups  take  part  in 
fulvate  complex  formation  processes  [12-16]. 
P^h,cooh  changes  from  4.1  to  4.8,  and  pKRjh-0H 
changes  from  9.3  to  10.5  [12,  17,  18]. 

In  spite  of  researches,  experimental  data  on  stability 
constants  of  complex  compounds  of  FA  are 
heterogeneous.  One  of  the  reasons  to  this  is  application 
of  various  average  molecular  weights  of  FA  associates. 
In  its  turn,  this  puts  the  correctness  of  calculation  of 
stability  constants  of  fulvic  complexes  under  question. 
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Therefore,  it  is  difficult  to  investigate  complexing 
and  sorption  processes  taking  place  in  natural  waters 
and  soils,  identify  migration  fonns  of  heavy  metals 
and  radio  nuclides,  evaluate  and  assess 
chemical-ecological  condition  of  natural  waters  and 
soils. 

The  authors'  objective  was  to  obtain  pure  samples 
of  FA,  to  study  complex  formation  processes  between 
heavy  metals  and  FA,  to  detennine  the  structure  and 
average  stability  constants  of  fulvic  complexes. 

In  this  paper,  the  authors  discuss  the  complex 
formation  process  of  nickel  fulvate.  According  to  the 
data,  found  in  the  literature,  stability  constants  (J3)  of 
nickel  fulvic  complexes  differ  in  several  lines  from 
each  other.  Value  of  lg  /?  changes  from  3.3  to  7.1 
[19-22]. 

The  complex  formation  process  was  studied  at  pH  = 
8.0  by  the  solubility  method.  The  suspension  of  old 
nickel  hydroxide  was  used  as  the  solid  phase. 

2.  Experiment  Methodology 

To   obtain  pure   samples   of  FA,   the  authors 


The  Determination  of  Average  Stability  Constant  of  Nickel-Fulvic  Acids 
Complex  at  pH  =  8.0  by  the  Solubility  Method 


345 


concentrated  the  water  of  the  river  Mtkvari  by  the 
frozen  method.  Filtered  water  samples  were  acidified 
to  pH  2  and  were  put  for  2  h  on  water  bath  at  60  °C, 
for  coagulation  of  humin  acids.  Then,  the  solution  was 
centrifuged  for  10  min  at  8,000  rpm.  To  isolate  FA 
from  centrifugate,  the  adsorption-chromotographic 
method  was  used.  The  charcoal  (BAU,  Russia)  was 
used  as  a  sorbent.  Desorption  of  amino  acids  and 
carbohydrates  were  performed  by  means  of  0.1  N  HC1. 
For  the  desorption  of  polyphenols,  the  90%  acetone  water 
solution  was  used.  The  elution  of  FA  fraction  was 
performed  with  0.1  N  NaOH  solution  [23,  24].  The 
obtained  alkalic  solution  of  FA  for  the  purification 
was  passed  through  a  cation-exchanger  (KU-2,  Russia) 
and  dried  under  the  vacuum  until  the  constant  mass 
was  obtained.  Then,  model  solutions  of  FA  was 
prepared.  The  required  pH  was  adjusted  by  the 
addition  of  hydrochloric  acid  or  sodium  hydroxide. 

0.1  mL  old  suspension  of  nickel  hydroxide  and 
increasing  quantity  of  standard  solution  of  FA  were 
placed  in  1 5  mL  capacity  fluroplastic  cylinders.  pH  = 
8.0,  ju=  0.01  (KN03),  V=  10  mL.  Solution  containing 
the  sediment  was  placed  in  thermostat  at  25  °C.  Then, 
it  was  stirred  in  a  mechanical  mixer  for  60  h  (until  the 
balance  was  achieved)  and  then  suspension  was 
filtered  through  the  membrane  filters  (Sinpor  N6, 
diameter  of  phores  400  nm).  In  filtrates,  nickel  was 
measured  by  atomic  absorption  spectrophotometer 
(Perkin  aimer  200). 


3.  Results  and  Discussion 

The  data  (Table  1)  show  that  in  line  with  the 
increase  of  concentration  of  FA  in  the  solution,  the 
concentration  of  nickel  increases  for  several  times  as 
well  due  to  formation  of  fulvic  complex. 

During  the  calculation  of  molar  concentrations  of 
FA,  the  fact,  that  FA  form  associates  in  water 
solutions,  which  value  of  average  molecular  mass  (Mw) 
depends  on  pH:  Mw  =  l,350pH-4,540  (pH  >  4)  [24]  was 
taken  into  consideration. 

In  the  diluted  water  solutions  the  forms,  of  nickel 
in  %  were  calculated  through  Eqs.  (l)-(3). 

Ni2+%=  100/(1  +AiOH[OH>Al(OH)20[OFf]2)  (1) 

NiOH+%  =  /  WOH"]100/(l+/WOH>#,i 


(oh)2°[OH"]2) 


(2) 


Ni(OH)2  %  =  jB 

Ni(OH)2 

°[OlT]2 1 00/(  1 +A,oH[OFf ]+#,  i(OH)2  °[OFf]2)  (3) 
Al0H  =  5.5  x  104  ^i(oh)2°  =  5.5  x  1014  [25] 
As  the  calculations  show  ((Fig.  1),  in  balanced 
solutions  at  pH  =  8,  practically  the  whole  Ni(II)  is  in 
the  form  of  nickel  dihydroxocomplex  Ni(OH)2°. 

So  the  complex  formation  reaction  (pH  =  8.0)  may 
be  written  in  Eq.  (4). 

Ni(OH)2°  +  mFA2"  ±5  Ni(OH)2FA2m"  (4) 

In  Ref.  [26],  the  method  of  calculation  of  the  stability 
constant  of  fulvic  complexes  was  described  in  detail,  so 
in  this  paper  the  item  will  be  discussed  briefly.  The 
numeral  meaning  (m)  of  the  stechiometral  coefficient  or 
the  ligands    number  in  the  inner  coordination 


Table  1 

The  data  (mol/L)  required  for  calculating  stability  constants  of  nickel  fulvic  complex  (pH  =  8.0;  //  =  0.01  (KNO3),  t  = 

25  °C,  Mw 

(FA)  =  6,260,  C'N.  =  7.9  x  10"6). 

C  -C 

Si  ^  Ni 

[FA2"] 

/?=[CM-C'M]/{[C'M][FA2-]} 

6.2  x  1(V5 

5.0  x  10"5 

1.2  x  10~5 

5.3  x  10s 

7.2  x  1(V5 

5.9  x  10~5 

1.3  x  10'5 

5.8  x  105 

9.3  x  10"5 

7.4  x  10~5 

1.9  x  10-5 

4.9  x  10s 

1.1  x  10"4 

8.9  x  10'5 

2.1  x  10"5 

5.4  x  105 

1.2  x  10"4 

9.7  x  10"5 

2.3  x  10"5 

5.3  x  105 

1.7  x  10-4 

13.8  x  10-5 

3.2  x  10~5 

5.5  x  105 

2.0  x  1(V4 

16.0  x  10"5 

4.0  x  10'5 

5.1  x  10s 

/?Ni(OH)2FA2"  =  5.3x  105. 
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Fig.  1  The  distribution  diagram  of  the  forms  of  Ni(II)  in  diluted  water  solutions  (in  %). 
Table  2   The  data  to  calculate  the  meaning  of  m  by  the  least  square  method. 


Xi 

lg  C 


-4.2076 

-4.3010 

17.7038 

18.0968 

-4.1426 

-4.2291 

17.1611 

17.5194 

-4.0315 

-4.1307 

16.2529 

16.6529 

-3.  9586 

-4.0506 

15.6705 

16.0347 

-3.9208 

-4.0132 

15.3726 

15.7349 

-3.7695 

-3.8601 

14.2091 

14.5506 

-3.6989 

-3.7958 

13.6818 

14.0402 

ZXi  =  -27.7295,  Zyi  = 

-28.3805,  Ejc,2  =  1 10.05 18,  Ex,  >•,  = 

1 12.6295,  (Ik,-)2  =  768.9251,  Zx;  E>-,  = 

786.9770,  m  =  0.99. 

sphere  of  complex  equals  to  tangens  of  tilt  angle  of 
straight  line  built  in  coordinates: 

lg  [Cm-CNi]    mlg  CFA  (5) 

where,  C  m  is  the  initial  concentration  of  nickel,  CNi  is 
the  final  concentration  of  nickel,  received  after 
formation  of  complex  and  Cfa  is  concentration  of  FA 
in  solution. 

To  calculate  the  exact  value  of  tangens,  tilt  angle  of 
straight  line  was  used  the  least  square  method  (Table  2): 
m  =  tg a  =  (nl.xiyi  -  Xx.Zy,.)/  (raXx2  -  (1.x f)  (6) 
xi  =  lgCFA  and  y.  =  lg[CM  -  C  M] 
After  the  calculations  we  obtained  the  numeral 


meaning  of  stechiometral  coefficient  (m),  which 
equals  to  0.99,  which  demonstrates  that  complex  with 
the  structure  1:1  dominates  in  the  Ni(OH)2-H20-FA 
system.  Therefore,  the  formula  for  calculating  the 
stability  constant  can  be  expressed  in  the  following 
way: 


P=  [CM-CM]/{[CNi][FA2-]} 


where, 


[FA2"] 


(7) 


(8) 


Taking  into  account  the  average  molecular  mass  of 
associate  at  pH  =  8,  (Mw  =  1,350  pH-4,540  [24]),  the 
average  stability  constant  of  nickel  fulvic  complex 
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was  calculated  on  the  basis  of  experimental  data 
(Table  1). 

So,  in          Ni(OH)2  (solid)-Ni(OH)2° 

(solution)-FA  "-H20  system  at  pH  =  8  dominates  the 
nickel  dyhydroxofulvic  complex  with  the  structure  1:1, 
[Ni(OH)2FA]2",  which  average  stability  constant  is: 
M(OH)2FA2"  =  5.3  x  105. 

4.  Conclusions 

So,  in         Ni(OH)2  (solid)-Ni(OH)2° 

(solution)-FA2"-H20  system  at  pH  =  8  dominates  the 
nickel  dyhydroxofulvic  complex  with  the  structure  1:1, 
[Ni(OH)2FA]2",  which  average  stability  constant  is: 
/?Ni(OH)2FA2"  =  5.3  x  105. 
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Abstract:  This  study  reports  preparation  of  glass  composition  (54.50  wt.%)  Si02,  (10.80  wt.%)  B203,  (14.20  wt.%)  Na20,  (1.20 
wt.%)  K20,  (6.00  wt.%)  CaO,  (4.00  wt.%)  Fe2C>3  and  (9.30  wt.%)  Ti02  by  melt  quenching  method  using  direct  microwave  heating 
and  conventional  resistive  heating.  Study  of  dielectric  loss  factor  of  the  glass  as  function  of  temperature  illustrated  increasing  loss 
factor  above  370  °C,  550  °C,  650  °C  and  900  °C,  indicating  enhanced  microwave  absorption  by  the  glass  at  above  these  temperatures. 
Chemical  analysis  results  of  both  the  glasses  depicted  more  volatilization  loss  of  volatile  ingredients  in  conventional  heating.  The 
study  of  chemical  durability  was  performed  from  leachate  analysis  describing  less  leaching  of  Na20,  K20  and  other  constituents  from 
glass  melted  in  microwave  furnace.  Glass  transition  temperatures  (Tg)  were  found  to  be  576.3  °C  and  569.5  °C  for  glass  melted  in 
conventional  and  microwave  heating  route,  respectively.  Laboratory  experiment  of  glass  melting  utilizing  microwave  energy  as  an 
alternate  heating  source  demonstrated  70%-75%  electrical  power  saving. 

Key  words:  Microwave  heating,  energy  efficient  glass  melting,  dielectric  constant,  loss  factor,  chemical  durability. 

~  500  °C  before  direct  microwave  heating  for  melting  [4]. 
Feasibility  study  of  glass  melting  by  microwave 
heating  has  been  reported  in  different  glass  systems 
such  as,  niobium  phosphate  [5],  alumino-phosphate  [6], 
molybdo-phosphate  and  tungsto-phosphate  [7],  iron 
phosphate  [8]  and  inorganic  glasses  [9]  etc..  However, 
most  of  researches  were  carried  out  in  kitchen 
microwave  oven  in  which  proper  control  of  process 
parameter  could  not  be  done. 

Present  study  reports  preparation  of  glass 
composition  (54.47  wt.%)  Si02,  (10.81  wt.%)  B203, 
(14.16  wt.%)  Na20,  (1.19  wt.%)  K20,  (6.02  wt.%) 
CaO,  (4.05  wt.%)  Fe203  and  (9.30  wt.%)  Ti02  using 
microwave  heating.  Some  of  the  physical  and  optical 
properties  were  measured  and  compared  with  the  glass 
melted  from  identical  batch  in  resistive  heating.  The 
study  of  power  consumption  in  microwave  heating 
and  in  conventional  resistive  heating  were  presented 
and  discussed. 


1.  Introduction 

Microwave  heating  has  find  its  application  in 
material  processing  mainly  because  of  rapid 
volumetric  heating,  flexibility  in  furnace  operation 
(i.e.,  switching  off  at  maximum  temperature)  and 
comfortable  working  enviromnent  owing  to  only 
sample  heating  keeping  the  cavity  cold,  substantial 
reduction  in  energy  consumption,  improved 
mechanical  properties  and  many  others  [1,  2].  The 
difficulty  of  glass  melting  using  microwave  energy  is 
mainly  due  to  low  dielectric  loss  factor  of  major  raw 
material  at  room  temperature  making  the  material  as 
microwave  transparent  material.  However,  the  loss 
factor  of  glass  increases  with  temperature  enabling 
microwave  coupling  with  materials  above  "critical 
temperature"  and  thereby  generating  heat  within  it  [3]. 
Hence,  glass  batches  are  heated  conventionally  above 
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2.  Experiments 

2. 1  Conventional  Melting  of  Glass 

The  glass  was  prepared  using  high  purity  raw 
materials  like  quartz,  H3B03,  Na2C03,  K2C03,  CaC03, 
Fe203  and  TiC>2.  Properly  mixed  raw  materials  for  60  g 
glass  was  uniaxially  pressed  to  form  pellet  in  a 
hydraulically  operated  pellet  press  and  the  pellet  was 
placed  in  a  50  mL  alumina  crucible  which  was 
subjected  to  heating  in  electrical  resistive  furnace 
(raising  hearth  furnace  from  Bysakh  &  Co.,  India; 
heating  zone:  230  m  x  230  m  x  300  m).  B-type 
thermocouple  was  fixed  with  the  furnace  to  monitor 
the  temperature.  Power  consumption,  temperature  and 
time  were  recorded  to  generate  T-t-P 
(temperature-time-power)  profile  in  resistive  heating. 
Melting  was  carried  out  at  1,300  °C  for  2  h  with  a  heat 
up  rate  4  °C/min  and  manual  stirring  was  adopted 
twice  with  30  min  interval  using  silica  rod  to  yield 
homogeneous  glass  prior  to  pouring  into  a  preheated 
steel  mold.  The  hot  glass  block  was  annealed  in  a 
muffle  furnace  at  560  °C,  below  glass  transition 
temperatures  Tg,  for  2  h  followed  by  controlled 
cooling  to  room  temperature.  The  annealed  glass  was 
termed  as  SBNRH  (glass  prepared  in  resistive 
heating). 

2.2  Microwave  Melting  of  Glass 

Identical  pellet  of  glass  batch  for  60  g  glass  was 
placed  inside  an  alumina  crucible  insulated  by 
microwave  transparent  cage  insulating  box  and  placed 
inside  a  3  kW,  2.45  GHz,  multimode  microwave 
furnace  (from  M/S  Enerzi  Microwave  Systems  Pvt. 
Ltd.,  India,  microwave  cavity  size  367  mm  x  349  mm 
x  386  mm).  An  opening  of  20  mm  diameter  was  made 
at  top  cover  of  insulation  box  to  allow  temperature 
monitoring  through  a  non  contact  IR  (Infrared) 
pyrometer  (260-1,800  °C).  The  temperature  measurement 
accuracy  and  repeatability  was  ±  0.3%  of  the 
measured  value  +  1  °C  with  adjustable  emissivity  (e: 
0.1-1.0).  MW  power  (Microwave  output  power)  and 


temperature  were  recorded  by  DAQSOFT  software  in 
a  separate  computer.  A  separate  energy  meter  (kW-h 
meter)  was  also  connected  to  the  input  of  the  electrical 
circuit  to  monitor  mainly  consumption  of  instant 
electrical  power  and  total  power  consumption  for 
overall  system  including  utilities.  An  integrated  closed 
loop  water  cooling  system  was  used  to  cool  the 
microwave  cavity,  magnetron  housing  (model: 
2M265-M12WJ-Panasonic)  and  others.  The  melting 
of  the  studied  glass  was  carried  out  in  air  atmosphere. 
Fumes,  hot  gas  generated  during  melting  of  batch  was 
exhausted  out  intermittently  by  a  diaphragm  pump. 
Glass  melting  was  carried  out  at  1,350  °C  for  1  h. 
Molten  glass  was  cast  into  a  preheated  mold  and 
annealed  as  explained  earlier  section.  The  annealed 
glass  was  termed  as  SBNM  (glass  melted  in 
microwave  furnace). 

2.3  Characterization 

XRD  (X-ray  diffraction)  analysis  was  performed  on 
94  X'Pert,  PANalytical  using  Ni-filtered  CuKa 
radiation  with  wavelength  of  1.5406  A  to  confirm  the 
sample  amorphous  characteristic.  The  scanning  was 
carried  out  in  the  range  of  10°-85°  with  a  step  size  of 
0.05  7min. 

ICP-AES  (inductively  coupled  plasma  atomic 
emission  spectroscopy)  (model:  Spectro  Ciros  Vision, 
Germany)  was  used  for  chemical  analysis  of  Fe203, 
TiC>2,  Na20  and  K20;  volumetric  titrimetry  for  SiC>2 
and  B203  and  complexiometry  for  Al203  and  CaO. 
Chemical  durability  was  evaluated  from  leachate 
analysis  by  ICP-AES.  The  bulk  glasses  were 
pulverized  and  sheaved  (within  mesh  No.  40-60). 
Then,  the  glasses  were  washed  in  acetone  and  dried  at 
110  °C  for  2  h.  The  washed  and  dried  glasses  were 
cooled  down  inside  a  desiccator.  Each  sample  of 
weight  5  g  was  taken  into  a  Teflon  beaker  with  100  mL 
deionised  water  and  kept  at  90  °C  for  14  days.  During 
the  treatment  in  the  above  mentioned  condition,  some 
of  the  constituents  of  the  glasses  were  leached  out  into 
the  water  and  these  were  measured  by  ICP-AES 
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(model:  Spectro  Ciros  Vision,  Germany).  The  relative 
standard  deviation  of  three  successive  measurements 
for  all  the  constituents  was  in  the  range  from 
0.59%-0.74%. 

Dielectric  loss  factor  and  dielectric  constant  were 
measured  at  2.45  GHz  as  a  function  of  temperature 
(up  to  1,000  °C).  A  cylindrical  cavity,  a  vector 
network  analyzer,  a  microwave  amplifier,  an  infrared 
pyrometer  and  a  computer  were  employed  to  carry  out 
the  measurement.  The  cylindrical  rod  shape  sample 
was  placed  inside  a  quartz  tube  containers  and 
microwave  are  applied  for  heating  and  measuring 
simultaneously  inside  a  microwave  cavity  [10]. 
Typical  uncertainty  of  these  measurements  is  around  2% 
in  real  part  and  5%  in  imaginary  part.  For 
high-temperatures  (around  1,000  °C),  the  uncertainty 
can  increase  to  10%  because  the  dimensions  of  the 
sample  change  when  the  materials  are  melting.  The 
dimension  of  the  sample  is  a  variable  used  in  the 
dielectric  calculations. 

The  Tg  was  identified  by  DSC  (differential  scanning 
calorimeter)  (NETZSCH  STA  449F3)  at  a  heating  rate 
of  10  K/min  in  flowing  nitrogen  with  an  accuracy  ±  1  °C. 
The  refractive  indices  of  glass  samples  prepared  by 
both  methods  were  measured  at  five  different 
wavelengths  (473  nm,  532  nm,  633  nm,  1,064  nm  and 
1,552  nm)  on  Prism  Coupler  (Metricon  Model-20 1 0, 
NJ,  USA)  fitted  with  five  different  lasers  as 
illuminating  sources.  Three  measurements  were 
recorded  for  each  sample  and  the  standard  deviation 
was  reported  as  error  of  measurement  (±  0.0004). 

FTIR  (Fourier  transform  infrared  spectroscopy) 
measurements  of  both  the  glasses  were  performed  on  a 
FTIR  spectrophotometer  (model:  1615,  Perkin  Elmer 
Corporation,  USA)  with  a  spectral  resolution  of  4  cm"1 
over  a  spectral  range  of  400-2,000  cm"1  to  see  the 
structural  changes,  if  any. 

Longitudinal  velocity  (VL)  and  shear  velocity  (VT) 
were  measured  by  the  ultrasonic  pulse-echo  method, 
using  Ultrasonic  Flaw  Detector  (model  No.  EPOCH 
1000,  Blue  Star).  An  ultrasonic  pulse  of  20  MHz 


frequency  employed  through  longitudinal  and  shear 
transducers.  The  Poisson's  ratio  (v),  Young's  (£), 
shear  (G)  and  bulk  (K)  modulus  were  calculated 
applying  the  measured  bulk  density  (p),  VL,  and  VT  in 
Eqs.(l)-(4)[11,  12]. 


G=p-V*  (3) 
K  =  l  -OVl-  4Fr2)  (4) 

2.4  T-t-P  Profile 

Temperature  and  variation  of  MW  power  with  time 
for  melting  of  60  g  studied  glass  using  direct 
microwave  heating  was  presented  in  Fig.  1.  MW 
power  was  manually  controlled  to  see  the  microwave 
interaction  with  sample.  A  rapid  temperature  increase 
was  observed  after  30  min  with  slow  increase  of  MW 
power  at  0.6  kW.  A  sharp  increase  of  T-t  profile  at 
900  °C  was  seen  in  the  Fig.  1,  predicting  higher  dielectric 
loss  factor  of  the  glass  above  this  temperature.  Fig.  1 
also  depicted  maximum  MW  power  ~  0.60  kW  during 
the  melting  till  1,350  °C.  However,  efficiency  of 
microwave  generator  was  within  the  range  55%-65% 
of  input  electrical  power.  Thus  maximum  input 
electrical  power  would  be  ~  1  kW  during  the  melting 
of  studied  glass  in  microwave  furnace.  Fig.  1  further 
illustrated  that  total  time  was  within  3  h  during 
melting  including  30  min  initial  start  up  time  prior  to 
initiation  of  heating.  Further,  operation  flexibility  in 
microwave  furnace  enables  easy  switch  off  at  high 
temperature  reducing  further  time  and  energy. 

Temperature  and  variation  of  input  power  with  time 
for  the  studied  glass  melted  in  the  resistive  heating 
furnace  was  displayed  in  Fig.  2.  Heating  rate  in 
conventional  furnace  was  restricted  to  4  °C/min  to 
protect  the  life  of  heating  element  and  refractory 
lining  inside  the  furnace.  Soaking  time  at  maximum 
temperature  in  conventional  melting  was  almost 
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Fig.  1  Temperature  and  variation  of  microwave  output 
power  with  time  during  microwave  melting  of  glass. 
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Fig.  2  Temperature  and  electrical  power  profile  with  time 
for  conventional  heating  in  resistance  heating  furnace. 

double  the  time  needed  for  microwave  processing 
resulting  total  melting  time  was  around  ~  7  h.  Fig.  2 
also  illustrated  that  maximum  instant  electrical  power 
in  electrical  raising  hearth  furnace  was  ~  4.00  kW 
almost  four  times  than  that  of  microwave  heating. 

A  comparison  of  total  power  consumption  for  both 
the  heating  method  using  above  furnaces  for  the  equal 
amount  of  studied  glass  was  illustrated  in  Fig.  3.  It 
was  clearly  seen  from  the  Fig.  4  that  the  total  power 
consumption  was  ~  26  kWh  in  the  conventional 
laboratory  glass  melting  process  using  above  raising 
hearth  furnace.  Overall  melting  time  was  found  to  be  7  h1 
in  conventional  glass  melting  process.  In  microwave 
melting,  the  total  power  consumption  was  recorded  ~ 
6  kWh  incorporating  additional  power  load  for 


utilities,  i.e.,  chilled  water  system,  fume  exhaust 
system.  However,  total  power  requirement  could  be 
further  minimised  by  changing  the  design  of  the 
microwave  furnace  using  air-cooled  magnetron  and  by 
placing  a  small  amount  of  microwave  susceptor  like 
SiC  around  the  crucible,  which  preheats  the  raw 
material  before  directly  absorbing  microwave 
radiation.  At  higher  temperature,  microwaves  interact 
with  the  entire  volume  of  glass  at  the  same  rate 
causing  volumetric  heating,  which  further  reduces 
holding  time  at  maximum  temperature  in  glass 
melting  compared  to  conventional  glass  melting. 
Since  glass  melting  is  a  huge  energy  consuming 
process,  it  would  be  of  great  relief  if  alternate  method 
of  producing  glass  reduces  energy  requirement 
substantially.  In  this  context,  microwave  glass  melting 
could  be  a  potential  energy  efficient  technology 
saving  a  significant  amount  of  electrical  energy  and 
time,  needed  for  present  glass  making.  Further, 
microwave  melting  of  glass  would  eventually 
contribute  to  the  cleaner  environment.  However,  the 
difficulty  remains  in  microwave  glass  melting  with 
proper  melting  crucible  and  design  of  suitable 
microwave  equipment  for  larger  scale  glass  melting. 

3.  Results  and  Discussion 

3. 1  Microwave  Absorption  of  Material 

Microwave  heating  is  the  result  of  absorption  of 
electromagnetic  energy  by  the  material  exposed  to  the 
microwave  field  distributed  within  a  reflective  cavity. 
It  is  based  on  the  power  absorbed  per  unit  volume  (Eq. 
(5)  [13,  14]). 


P  =  a\E2\  =  Infs^s \  tan 5E2 


(5) 


where,  E  is  the  magnitude  of  the  internal  electric  field, 
e0  is  the  permittivity  of  free  space,  /  is  the  microwave 
frequency,  a  is  the  total  effective  conductivity,  s'r  is 
the  relative  dielectric  constant,  and  tand  is  the  energy 
loss  required  to  store  a  given  quantity  of  energy,  s"e/f 


'Furnace  cooling  down  to  room  temperature  was  not 
considered. 
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Fig.  3  Total  power  consumption  profile  with  temperature; 
for  resistance  heating  and  for  microwave  heating. 
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Fig.  4  Dielectric  constant  and  loss  factor  for  the  studied 
glass  was  measured  around  2.45  GHz  as  a  function  of 
temperature  up  to  1,000  °C. 

(=  e'/tarkS)  is  the  relative  effective  dielectric  loss 
factor.  It  is  clearly  seen  from  the  above  Eq.  (5)  that  the 
dielectric  properties  (e0,  e'r  and  tanr5)  assume  a 
significant  role  in  determining  the  extent  of  power 
absorbed  by  a  material. 

Fig.  4  portrayed  the  dielectric  constant  and 
dielectric  loss  factor  as  function  of  temperature  (up  to 
1 ,000  °C)  for  this  glass.  Dielectric  loss  factor,  one  of 
the  important  parameters,  indicates  the  degree  of 
microwave  absorption  by  the  material  to  generate 
heating  within  it  for  a  fixed  frequency  (2.45  GHz)  and 
a  given  electrical  field  within  the  material.  The  Fig.  4 
depicted  that  loss  factor  increased  at  ~  370  °C,  ~  560  °C, 
~  650  °C  and  ~  900  °C  suggesting  enhance  microwave 
absorption  and  heating  of  glass  at  these  temperatures. 
Sharp  increase  of  loss  factor  0  at  around  900  °C 
corroborated  the  rapid  heating  observed  from  the  T-t 


profile  in  microwave  heating  in  Fig.  1.  However, 
microwave  interaction  with  batch  material  at  lower 
temperature  was  found  different  due  to  the  presence  of 
different  raw  material  (i.e.,  carbonate,  H3BO3). 

3.2  Comparison  of  Properties 

Fig.  5  displayed  photograph  of  1mm  thick  polished 
glass  (15  mm  x  20  mm)  melted  in  a  conventional 
resistance-heating  (SBNRH)  and  in  microwave 
heating  (SBNM).  In  Fig.  5,  improved  homogeneity  of 
SBNM  was  seen,  apparently  compared  to  SBNRH. 
Though  the  SBNM  was  prepared  without  employing 
any  mechanical  stirring,  improved  homogeneity  was 
due  to  the  result  of  microwave  induced  "self-stirring" 
effect  [15]. 

The  XRD  patterns  for  the  SBNRH  and  SBNM 
samples,  shown  in  Fig.  6,  depicting  identical  broad 
humps  confirmed  the  amorphous  nature  of  the 
samples.  The  absence  of  sharp  Bragg  peak  in  the  XRD 
profile  of  SBNM,  similar  to  that  of  SBNRH, 
confirmed  no  unmelted  crystalline  phase  left  within 
the  glass  melted  using  microwave. 

Chemical  analysis  results  for  SBNM  and  SBNRH 
samples  were  listed  in  Table  1 .  B203  and  Na20  were 
found  less  in  SBNRH  glass  illustrating  higher 
evaporation  loss  of  these  ingredients  in  conventional 
melting.  A1203  was  found  6.47%  in  SBNRH  and  1 .56%  in 
SBNM  sample.  In  conventional  heating,  heat  travels 
from  the  heating  element  to  the  material  by  means  of 
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Fig.  5  Photographs  of  polished  SBNRH  and  SBNM  glass 
samples. 
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Fig.  6    XRD  curves  for  SBNM  and  SBNRH  samples. 

radiation,  conduction,  convection  mechanism  and  thus 
temperature  of  the  crucible  wall  remains  higher  than 
the  molten  glass.  This  higher  wall  temperature 
together  with  higher  melting  time  contributed  higher 
leaching  of  AI2O3  into  the  molten  glass  in 
conventional  melting.  In  microwave  heating,  heat 
generates  within  the  material  internally  resulting  in 
volumetric  heating  within  the  glass  without  heating 
the  alumina  crucible  owing  to  its  low  dielectric  loss 
factor.  Due  to  this  rapid  internal  heating,  temperature 
gradient  in  microwave  heating  exists  in  the  reverse 
direction  than  that  of  conventional  heating  resulting 
less  temperature  in  crucible  wall  and  in  the  top  surface 
of  molten  glass.  Thus,  less  crucible  wall  temperature 
along  with  less  melting  time  favour  low  contamination 
of  alumina  into  the  glass  melt  in  microwave 
heating.  This  reverse  temperature  gradient  also 
contributes  less  evaporation  of  volatile  ingredients  of 
glass  in  microwave  heating.  However,  slightly  higher 
melting  temperature  is  needed  in  case  of  microwave 


melting  due  to  less  surface  temperature  of  molten 
glass. 

The  chemical  durability  of  both  the  glasses  was 
evaluated  from  leachate  analysis  by  Inductively 
Coupled  Plasma-Atomic  Emission  Spectroscopy.  The 
results,  as  shown  in  the  Table  2,  illustrated  higher 
concentration  of  AI2O3,  CaO,  Fe203,  TiC>2,  Na20  and 
K20,  leached  out  into  water  from  the  SBNRH  sample 
signifying  higher  leaching  rate  of  above  constituents 
from  SBNRH  sample.  This  indicated  improved 
chemical  durability  of  glass  was  prepared  in 
microwave  heating.  This  may  be  explained  that  more 
evaporation  of  glass  former,  i.e.,  SiC>2  and  B2O3  the 
main  building  block  in  glass,  weaken  the  glass 
structure  in  conventional  heating  method. 

Fig.  7  exemplified  DSC  thermograph  for  both  the 
glasses  SBNRH  and  SBNM.  Tg  was  identified  569.5  °C 
and  576.3  °C  for  SBNM  and  SBNRH  samples, 
respectively.  Tg  value  for  SBNM  was  observed  less 
due  to  difference  in  the  composition  of  the  glasses 
resulting  from  low  loss  of  volatile  constituents  of 
glass.  Higher  presence  of  AI2O3  in  SBNRH  owing  to 
leaching  from  crucible  wall,  also  contributed  the 
higher  Tg  in  SBNRH. 

The  some  of  the  physical,  optical  and  elastic 
properties  of  SBNM  and  SBNRH  were  listed  in  Table  3. 
Optical  properties  for  SBNRH  and  SBNM  samples 
were  calculated  from  their  measured  refractive  indices 
using  relevant  expressions  and  were  presented  in 
Table  2. 

Fig.  8  portrayed  dispersion  curves,  obtained  by 
fitting  the  measured  refractive     indices  with  the 


Table  1    Chemical  analysis  results  for  SBNRH  and  SBNM  samples. 

Composition  (wt.%) 

SBNRH 

SBNM 

Si02 

56.36  (±  0.08) 

58.88  (±  0.08) 

B203 

7.92  (±  0.04) 

9.13  (±0.04) 

Na20 

10.92  (±  0.04) 

12.02  (±  0.04) 

K20 

0.68  (±  0.03) 

0.82  (±  0.03) 

CaO 

4.82  (±  0.04) 

5.03  (±  0.04) 

Fe203 

3.65  (±0.03) 

3.65  (±0.03) 

Ti02 

8.91  (±0.04) 

8.57  (±  0.04) 

A12Q3 

6.47  (±  0.03) 

1.56  (±0.03) 
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Table  2    Glass  transition  temperature,  optical  properties  and  elastic  properties  for  SBNRH  and  SBNM  are  presented. 

Property 

SBNRH 

SBNM 

Measured  Refractive  index  (n)  at 

473  nm  (±  0.0004) 

1.5969 

1.6006 

532  nm  (±  0.0004) 

1.5904 

1.594 

633  nm  (±  0.0004) 

1.5827 

1.586 

1,064  nm(±  0.0004) 

1.5691 

1.573 

1,552  nm  (±  0.0004) 

1.5624 

1.5656 

Estimated  Refractive  indices  at 

nd.587.6  (±  0.0005) 

1.5860 

1.5895 

nF.mA  (±  0.0005) 

1.5953 

1.5988 

nc-656.3  (±  0.0005) 

1.5814 

1.5848 

n<i-i 

0.586 

0.59 

Mean  dispersion  (iif-iic) 

0.014 

0.014 

Abbe  No. 

42.16 

42.10 

Reflection  loss  (%) 

5.17 

5.13 

Elastic  properties 

Longitudinal  velocity  (m/s) 

5,776 

5,813 

Shear  velocity  (m/s) 

3,448 

3,505 

Young's  modulus  (E,  GPa) 

78.90 

80.20 

The  Poisson's  ratio  (v) 

0.219 

0.235 

Shear  modulus  (G,  GPa) 

90.07 

93.76 

Bulk  modulus  (K,  GPa) 

46.90 

50.46 

Glass  transition  temperature,  Tg  (±  1  °C) 

576.3 

569.5 
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Fig.  7    DSC  thermographs  for  SBNRH  and  SBNM  sample. 

Cauchy  Eq.  (6).  The  refractive  indices  at  standard 
wavelengths  rid  (at  587.6  nm),  nF  (at  486.1  nm)  and  «c 
(at  656.3  nm)  were  estimated  from  the  curve. 


n(A)  =  A  +  ^r  +  ^r  +  - 

A2  A2 


(6) 


where,  n  is  refractive  index  at  wavelength  X.  The  Abbe 
number  V  of  a  material  is  defined  as  Eq.  (7): 


_  (nd-\) 
(nF-nc) 


(7) 


where,  rid,  nF  and  ric  are  the  refractive  indices  of  the 
samples  at  the  wavelengths  587.6  nm,  486.1  nm  and 
656.3  mn,  respectively.  Abbe  number  was  calculated 
and  found  very  similar  for  both  SBNRH  and  SBNM 
as  42.16  and  42.10,  respectively.  Mean  diffraction  and 
reflection  loss  for  SBNRH  and  SBNM  samples  were 
also  found  similar. 

Elastic  properties  like  v,  Young's  (E)-shear 
(G)-bulk  (K)  modulus  were  calculated  and  listed  in 
Table  1 .  The  identical  or  slightly  higher  value  of  these 
properties  in  the  microwave  prepared  glass  confirmed 
that  similar  glass  or  slightly  improved  property  of 
glass  could  be  prepared  in  microwave  heating. 

The  FTIR  reflectance  spectra  of  samples  SBNRH 
and  SBNM  were  presented  in  Fig.  9.  Both  the  spectra 
of  these  glasses  consisted  of  similar  broad  and  sharp 
bands  in  different  regions  (400-2,000  cm"1) 
confirming  identical  structure.  FTIR  spectra  for  both 
glasses  showed  the  presence  of  silicate  peak  centered 
at  1,022  cm"1  and  Si-O-Si  bend  at  around  443  cm"1  [16]. 
The  bands  at  around  1,020  cm"1  peak  were  attributed 
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Table  3  The  concentration  (in  ppm)  of  different  constituents  present  in  deionised  water  after  14  days  insertion  of  SBNRH 
and  SBNM  glasses  at  90  °C. 


Leachate  from  samples 

A1203 

CaO 

Fe203 

Ti02 

Na20 

K20 

SBNM 

31 

19 

10 

15 

131 

21 

SBNRH 

59 

68 

21 

35 

377 

74 

A -SBNRH 
B ■ SBNM 


400  600  800  1000  1200  1400  1600 

Wavelength  (nm) 

Fig.  8  RI  vs.  wavelength  fitted  by  Cauchy  curve  fitting 
equation. 


pi 


5= 

(2 


2000 

Wavenumber  (cm"1) 
Fig.  9    FTIR  reflectance  spectra  of  studied  samples  in  the 
wave  number  range  400-2,000  cm"1. 

to  a  stretching  vibration  of  B-O-Si  linkage  [17].  The 
peak  ~  700  cm"1  was  due  to  the  bending  of  B-O-B 
linkages  in  the  borate  glassy  network  [18].  The 
boroxol  ring  appeared  at  780  cm"1  [19].  The  band 
around  1,360  cm"  was  assigned  to  B-0  stretching 
vibration  of  trigonal  [BO3]  unites  [20,  21]. 

4.  Conclusions 

A   glass   containing  Si02-B203-Na20-K20-CaO- 


Fe203-Ti02  was  successfully  melted  using  direct 
microwave  heating.  Dielectric  constant  and  loss  factor 
as  function  temperature  depicted  increase  of  these 
properties  at  370  °C,  550  °C,  650  °C  and  900  °C, 
indicating  enhanced  microwave  absorption  with  this 
glass  composition  above  these  temperatures.  Similar 
glass  was  also  prepared  from  identical  batch  in 
conventional  resistive  heating.  Physical  properties  and 
structure  of  both  the  glasses  were  found  similar.  Glass 
melted  in  microwave  heating  was  found  closer 
chemical  composition  as  prescribed  in  batch, 
improved  chemical  durability  and  lower  glass 
transition  temperature.  Substantial  less  consumption 
of  energy  and  shorter  processing  time  was 
demonstrated  in  laboratory  scale  glass  melting  in 
microwave  furnace.  Thus,  microwave  processing  of 
glass  could  be  a  potential  method  for  reduction  of 
energy  consumption  and  time. 
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Electrochemical  Behavior  of  Electrodeposited  Sn  Films: 
Possible  Negative  Electrode  for  Na+  Rechargeable 
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Abstract:  Tin  films  on  copper  substrate,  obtained  by  electrodeposition  procedure,  were  structural  and  electrochemical  characterized. 
In  particular  to  investigate  the  possibility  to  use  such  metal  as  possible  negative  electrode  in  Na+  rechargeable  batteries,  EPS 
(electrochemical  potential  spectroscopy)  and  galvanostatic  charge/discharge  cycling  of  the  electrodes  were  investigated,  at  room 
temperature  in  organic  electrolyte.  Three  crystalline  and  one  amorphous  phases  were  identified  as  well  as  high  discharge  capacity 
(738  mAh/g)  was  obtained  after  4  cycles.  Unfortunately  material  fading,  due  to  the  internal  stress  during  sodiation/desodiation 
process,  causes  poor  cyclability. 

Key  words:  Sodium  ion  batteries,  anode  materials,  electrodeposed  tin,  potential  spectroscopy. 

to  only  active  materials  weigh,  for  such  a  Na  ion 
battery  would  be  around  275  Wh/kg,  which  is  about 
70%  of  commercially  available  Li-ion  battery,  so  far. 
This  limitation  is  mainly  due  to  two  factors:  The 
discharge  voltage  profiles  of  the  mixed  oxide  which 
would  reduce  the  operative  voltage  of  the  whole 
system  and  the  lower  specific  capacity  of  the  hard 
carbon/Na+  system  compared  to  the  lithiated  graphite. 
Moreover,  the  specific  power  of  the  battery  may  be 
limited  by  the  low  discharge  kinetic  of  hard  carbon. 
Similar  electrochemical  properties  of  Na  and  Li  show 
that  materials,  electrolytes,  and  cell  structures,  which 
have  been  used  in  Li-ion  batteries,  may  be  suitable  for 
Na  ion  rechargeable  batteries;  this  has  been  already 
proven.  Many  positive  electrodes,  constituted  of  some 
metal  oxides  (Na0.44MnO2  [7],  NaACo02  [8] 
NaJVlno.5Feo.5O2  [6]),  NASICON-type  materials 
(Na3V2(P04)3  [9]),  olivine  phases  (NaFeP04  [10]), 
and  others  (NaVPC^F  [11])  were  successfully  used. 
On  the  contrary,  even  if  there  are  some  studies  of 
possible  anode  materials,  a  good  negative  electrode  is 
still    one    of   the    fundamental    impediments  in 


1.  Introduction 

Since  first  commercial  production  of  Li-ion 
batteries  in  1991,  Li  row  material  has  been  highly 
demanded  to  power  not  only  portable  electronic 
devices  but  also  EVs  and  HEVs  (electronic  and  hybrid 
electric  vehicles).  This  issue  has  turned  Li  to  a 
strategic  and  expensive  material  [1]  which  is 
unfortunately  in  short  supply  and  it  may  not  meet  the 
industrial  needs  in  the  near  future.  This  concern  has 
motivated  scientists  to  develop  different  battery 
systems  and,  recently,  among  them  the  Na  ion 
rechargeable  batteries,  which  were  initially 
investigated  earlier  than  Li-ion  batteries  became 
commercially  available  [2-4].  Since  Na  is  one  of  the 
most  abundant  elements  in  the  Earth's  crust,  Na  -ion 
systems  can  really  become  alternative  for  Li-ion 
batteries.  Considering  carbon  [5]  and  Nax(x  = 
o.66)Mn0.5Feo.502  [6]  as  a  possible  anode  and  cathode, 
respectively,  the  energy  density,  calculated  according 
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progressing  Na-ion  batteries;  in  fact,  Na  does  not 
intercalate  between  typical  graphitic  carbon  sheets  [12]. 
Nevertheless,  other  different  materials  can  be  used  as 
negative  electrodes,  among  them  the  Na-metal  alloys 
showing,  in  phase  diagram,  wide  biphasic  regions  can 
play  an  important  role;  in  particular,  if  the 
composition  of  each  compound  has  high  amount  of 
sodium.  For  this  reason,  the  Na-Sn  system  was 
investigated;  in  fact,  as  reported  by  Sangster  and  Bale  [13], 
at  temperatures  lower  than  100  °C,  many  biphasic 
regions  with  compounds  having  large  amount  of 
sodium  (from  NaSn6  to  Nai5Sn4)  are  present.  In  such 
conditions,  invariant  and  high  chemical  potentials  of  Na 
have  to  be  present  and  electrodes  with  theoretical 
capacity  of  about  850  mAh/g  might  be  realized. 
However,  these  capacity  values  are  generally  observed 
for  a  very  few  cycles  in  conventional  composite 
electrodes.  Looking  at  the  wide  knowledge  coming 
from  the  investigation  of  such  class  of  materials  in 
lithium  ion  batteries,  the  first  source  of  electrode  failure 
is  the  decrepitation  due  to  the  large  volume  expansion 
and  contraction  during  the  alloying/dealloying  process, 
respectively.  In  fact,  the  commercial  available  binders 
are  not  able  to  accept  the  stress  induced  by  the  strain 
coming  from  the  reaction  with  the  consequent  electrode 
pulverization.  A  good  way  to  avoid  using  binder  is  to 
produce  Sn  layer  directly  onto  proper  current  collectors. 
This  can  be  made  by  several  methods,  among  them 
electro-deposition  is  one  of  the  most  convenient 
procedures  for  preparing  Sn  in  different  morphologies. 
In  addition,  this  technique  is  very  cheap  with  respect  to 
the  others,  and  it  does  not  need  any  advanced  facility. 
For  this  reason,  in  the  present  paper  the  charge-discharge 
characteristics  of  Sn  films  electrodeposited  onto  copper 
current-collectors  were  performed.  The  experiments 
were  carried  out,  at  room  temperatures,  using  three 
electrodes  T-shaped  Swagelok  cells  having  Na  as 
counter  and  reference  electrode  and  1  M  NaC104  in 
propylene  carbonate  as  electrolyte;  cyclability  and 
capacity  of  the  electrodes  were  investigated  between 
0.01  V  and  1  V  vs.  Na/Na+  electrode  potential. 


2.  Experiments 

2.1  Tin  Films  Preparation 

Tin  films  have  been  deposited  electrochemically  on 
copper  substrates  from  a  chloride  bath,  using 
SnCl2-2H20  and  tri-ammonium  citrate  [14]. 
Deposition  bath  was  made  of  50  g/L  tin  chloride(II) 
and  50  g/L  tri-ammonium  citrate,  dissolved  in  distilled 
water.  Platinum  net,  copper  foil  (1  cm  )  and  SCE 
(standard  calomelan  electrode)  were  used  as  counter, 
working  and  reference  electrodes,  respectively.  Before 
electro-deposition,  copper  substrates  were  degreased 
in  acetone  by  applying  ultrasonic  waves.  Depositions 
have  been  done  while  solution  was  stirring,  and 
current  density  was  24  mA/cm2  for  all  depositions. 

2.2  Characterization 

Microstructure  and  chemical  composition  of  Sn 
metal  films  have  been  characterized  by  SEM 
(scanning  electron  microscopy)  (Vega  TS5136  XM 
Tescan)  equipped  with  an  EDX  (energy  dispersive 
x-ray)  detector,  in  high  vacuum  configuration.  X-ray 
diffraction  measurement  was  carried  out  by  Bruker  D8 
advance  XRD  (X-ray  diffraction)  equipped  with 
secondary  beam  monochromator,  copper  source  (ka 
radiation),  and  operating  at  0.02  deg/step,  and  3 
sec/step. 

2.3  Electrochemical  Characterization 

Electrochemical  characterizations  were  carried  out 
using  three  electrodes  T-shaped  Swagelok  cells 
equipped  with  stainless  steel  current  collectors.  Sn 
films,  deposited  on  Cu  substrates,  were  sandwiched 
between  stainless  steel  and  Na  metal  disc,  the  latter 
pasted  on  the  top  of  current  collector.  Metallic  Na  was 
used  as  both  counter  and  reference  electrodes,  and 
trilayer  polypolprylene  membranes  were  used  as 
separators.  Separators  were  wetted  by  the  solution  of  1 
M  NaC104  in  propylene  carbonate,  which  was  used  as 
electrolyte.  In  order  to  remove  all  the  aqueous  remains 
present  in  the  Sn  films,  all  the  electro-deposited 
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samples  were  dried  in  vacuum  chamber  over-night, 
after  drying  at  130  °C  for  2  h.  All  potentials  are 
reported  vs.  Na/Na+  electrode  value,  and  the 
electrochemical  characterizations  were  performed 
using  Biologic  VMP3  multi-channel  battery  tester,  at 
room  temperature.  Typical  Sn  load  was  2.5  mg/cm2. 

3.  Results  and  Discussion 

3.1  Crystal  Structure  and  Morphology 

Fig.  la  shows  SEM  image  of  electrodeposited  Sn 
on  copper  current  collector  in  two  different 
magnifications. 

Landscape  image  and  cross  section  view  (Fig.  lb) 
show  very  good  coverage  of  Sn  on  Cu  substrate.  The 
Sn  tetragonal  morphology  can  be  easily  observed  at 
higher  resolution.  The  EDX  spectra  (Figs.  Id  and  le), 
corresponding  to  the  points  ai  and  a2  in  Fig.  la,  show 
the  presence  of  only  elemental  tin.  In  our  hypothesis, 
at  the  solution/current  collector  interface  a  Sn  film  of 
few  m  is  firstly  formed  acting  as  seeds  layer  for  the 
growth  of  highly  ordered  tetragonal  structures.  The 


total  film  thickness  is  about  38  m.  It  has  been  proved 
that  such  rough  morphology  shows  better  cyclic 
properties  with  respect  to  smooth  one  [15].  Fig.  lc 
shows  the  XRD  pattern  of  deposited  tin  film  on 
copper.  There  are  eight  peaks:  Six  of  them  can  be 
indexed  to  tetragonal  tin  [16],  with  a  lattice  parameter 
of  a  =  0.58329(8)  nm  and  c  =  0.31826(7)  nm;  the 
other  two  peaks  can  be  indexed  to  the  copper  current 
collector. 

3.2  Electrochemical  Analysis 

The  Sn  film  electrodes  were  assembled  into 
half-cells  with  Na  metal  as  counter  and  reference 
electrodes.  No  binders  or  conducting  carbon  were 
used.  The  half-cells  were  studied  using  EPS 
(electrochemical  potential  spectroscopy)  [17],  and 
galvanostatic  charge/discharge  cycling.  In  EPS 
measurement,  a  series  of  constant  potential  steps  were 
applied  to  the  electrochemical  cell.  On  each  step,  the 
cell  was  permitted  to  attain  quasi  equilibrium 
conditions  by  letting  the  current  decay  to  a  small  but 
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Fig.  1  (a)  SEM  image  of  electroplated  Sn  film  at  two  different  magnifications,  (b)  cross  section  view  of  the  layer,  (c)  XRD 
pattern  of  Sn  and  Cu  current  collector,  (d)  and  (e)  EDX  spectrums  carried  out  in  sites  a!  and  a2,  respectively. 


Electrochemical  Behavior  of  Electrodeposited  Sn  Films:  Possible  Negative  Electrode 
for  Na*  Rechargeable  Batteries 


361 


finite  value  (C/50,  typically  ~  17  mA/g  in  our 
measurements).  When  small  voltage  steps  (V  =  3  mV 
in  our  study)  are  applied,  the  charge  accumulated  on 
each  step  can  be  integrated  from  the  measured  current 
to  yield  an  accurate  derivative  of  the  voltage-charge 
relation.  This  derivative  relationship  is  a  very  useful 
tool  to  define  the  evolution  of  the  phase  equilibrium 
inside  electrode  materials.  In  fact,  under 
thermodynamic  equilibrium  condition  the 
potential/composition  profile  obeys  Gibbs  phase  rule. 
So,  the  first  derivative  of  quasi  equilibrium 
charge/voltage  curve  can  show  sharp  peaks  for 
two-phase  regions  or  bell-like  peaks  for  single-phase 
regions,  revealing  whether  a  phase  transformation 
occurs.  Fig.  2a  presents  potential  profiles  obtained  by 
EPS  for  first  two  cycles  of  Na-Sn  alloying/dealloying 
between  0.01  V  to  1.00  V.  It  is  obvious  that 
alloying/dealloying  reaction  is  not  completely 
reversible  as  charge-discharge  tails  are  not  fit  to  each 
others,  but  a  significant  amount  of  capacity  (82%  and 
90%  during  first  and  second  cycle,  respectively)  has 
been  recovered.  During  the  first  alloying  reaction,  the 
capacity  was  885  mAh/g,  which  shows  achieving  final 
composition  of  Nai5Sn4  (theoretical  capacity  847 
mAh/g). 

The  larger  observed  capacity  compared  to  the 


expected  one  may  be  attribute  to  the  formation  of 
passivating  layer  on  the  Sn  surface.  The  first  alloying 
reaction  shows  a  potential  profiles,  which  is  different 
from  the  others.  In  fact,  three  plateaux  are  observed  at 
0.160  V,  0.065  V  and  0.040  V,  which  correspond  to 
three  very  sharp  peaks  in  the  differential  capacity 
curve  (Fig.  2b).  According  to  the  equilibrium  phase 
diagram,  it  can  be  assumed  that  this  three  peaks 
correspond  to  the  formation  of  the  NaSn2,  NaSn  and 
Nai5Sn4  phases.  This  behavior  was  already  observed 
in  Sn  micrometric  particles  [18],  showing  successive 
strong  irreversible  de-alloyings.  On  the  contrary,  four 
different  plateaux  (labeled  1  -4  in  Fig.  2)  are  present  in 
our  profiles  during  de-alloying  reaction  in  both  first 
and  the  second  cycle.  These  plateaux  can  be  related  to 
the  different  biphasic  regions  present  in  Na-Sn  phase 
diagram  [13].  Recently,  L.  Ellis  et  al.  [19]  have 
reported  in  situ  X-ray  measurement  for  Sn  sodiation, 
and  demonstrated  the  existence  of  three  crystalline 
and  one  amorphous  phases,  corresponding  to  the  four 
plateaux  in  the  voltage  profile,  close  to  this  one 
reported  in  this  work.  Accordingly,  the  equilibrium 
between  Sn/NaSn3,  NaSn3/NaSn,  NaSn/Na9Sn4  and 
Na9Sn4/Na14Sn5  correspond  to  plateaux  1,  2,  3  and  4, 
respectively.  Fig.  3  shows  cycling  performance  of 
electroplated  tin  film  at  3  different  C-rates.  In  the 
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Fig.  2  (a)  Galvanostatic  sodiation  and  desodiation  curves  for  the  first  2  cycles;  (b)  same  experiment  reported  as  incremental 
capacity  vs.  potential. 
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Fig.  3  Discharge  capacity  vs.  cycle  number  at  three 
different  rates:  C/20  (black  squares),  C/10  (red  circles)  and 
C/5  (bleu  triangles). 

lower  rates  (C/20  and  C/10,  corresponding  to  42.5  and 
85  mA/g,  respectively)  range,  discharge  capacity  little 
increases  after  the  first  cycle,  most  likely  owing  due  to 
improvement  of  ions  diffusion  in  the  solid  state. 

After  the  4th  cycle,  we  have  seen  decrease  in 
discharge  capacity,  probably  due  to  active  material 
fading.  At  higher  rate  (C/5),  the  average  discharge 
capacity  for  the  first  6  cycles  is  475  mAh/g,  which  is 
much  less  than  that  at  lower  rates  and  highlight  that 
sodiation/desodiation  takes  place  partially.  In  addition, 
it  seems  that  C-rate  affects  directly  on  the  stresses 
induced  by  volume  expansion/shrinkage  taking  place 
during  the  alloying/dealloying  process. 

4.  Conclusions 

We  explored  the  possibility  to  use  tetragonal  tin 
layer,  obtained  by  very  cheap  method,  as  electrode  in 
Na+-ion  rechargeable  batteries.  The  films  show  very 
good  crystallinity  after  deposition  without  any  further 
heat  treatments.  Alloing/dealloying  voltage  profiles 
show  similar  behavior  to  that  previously  reported: 
Four  different  bi-phasic  regions  in  the  voltage  range 
between  0.01  V  and  1.00  V.  At  low  current  rates,  the 
specific  capacities  are  in  good  agreement  with  the 
theoretical  value  of  847  mAh/g,  but  their  values  drop 
after  4  cycles,  probably  due  to  active  material  fading. 


At  current  rate  of  170  mA/g  (C/5),  the  specific 
capacity  decreases  to  551  mAh/g,  but  remain  more 
stable  for  at  least  10  cycles.  As  for  similar 
alloying/dealloing  processes  (Si-Li  System),  metal 
nanostructure  thick  film  could  be  used  to  improve 
cyclability. 
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Abstract:  This  paper  discusses  the  vein  rocks  located  in  the  central  part  of  Khibina  Mts..  These  rocks  are  located  on  the  border  of 
two  zones:  ore  and  external.  In  the  ore,  zone  is  present  titanite-nepheline  bodies  whilst  in  the  outer  zone  are  massive-syenites 
"khibinites".  The  veins  can  be  divided  into:  mikrosyenites,  melteigites  and  other  rock  types  (e.g.,  tinguaites).  The  results  of  LA-ICP 
MS  (laser  ablation  of  inductively  coupled  plasma  mass  spectrometry)  indicate  that  these  rocks  are  associated  with  deep  zones  which 
formed  as  a  result  of  injection  of  primary  magma  and  contamination  of  fluids  during  migration,  where  there  have  hydrothermal 
processes.  Rock-forming  minerals,  their  corrosion  and  accessory  minerals  indicate  these  processes,  too. 

Key  words:  Khibina  Massive,  Kola  Peninsula,  vein  rocks,  geochemistry. 


1.  Introduction 

Khibina  Massive  is  a  central  intrusion,  located  near 
the  Imandra,  Umboziero  Gully  Lakes,  in  the  central 
part  of  Kola  Peninsula  on  the  northern  part  of  Baltic 
Shield.  These  intrusion  is  late  Devonian  [1-5], 
surrounded  by  Archean  migmatitic  gneisses  of  Kola 
series  and  metasediments  of  Imandra- Varzuga 
complexes  [6-8].  There  is  a  multi  zonal  massive  with 
characteristic,  specific  mineralization.  In  the  central 
part  is  located  an  ore  zone  with 
titanite-apatite-nepheline  bodies  which  accompanied 
alkaline  syenites-khibinites  from  the  occurrence  zone. 
These  syenites  are  massive  rocks  with  grey  color  and 
visible  alkali-feldspar  with  aegirinites  (Fig.  1).  In  the 
contact  of  these  zones,  they  are  present  a  numerous 
mineral  veins,  represented  by  the  tingualites, 
metleigites  and  microsyenites.  These  rocks  are  studied 
in  the  paper. 
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2.  Methods 

It  was  make  the  observation  in  the  Malaya  Belaya 
Valley  near  Ramzaya  Pass  and  in  other  region  at  the 
Bolshoi  Vudyavr  Lake.  The  150  samples  were  taken 
and  make  the  microscopic  observation  used  a 
polarysed  cross  and  reflected  light  (Leica  DM  2500P) 
and  electron  scanning  microscope  (Hitachi  SU6600) 
with  EDS  (energy  dispersive  X-ray  spectrometry).  In 
selected  samples  was  make  the  LA-ICP-MS  (laser 
ablation  inductively  coupled  plasma  mass  spectrometry) 
analysis  (model  ELAN  DRC  II,  Perkin  Elmer  with 
laser  microsonde  (LA)  LSX-500,  of  CETAC).  The 
microscopic  observation  was  made  in  the  Optical  and 
Electron  Microscpy  Laboratory  in  Geology  and 
Litosphere  Protection  Department,  Earth  Science  and 
Spatial  Management  Faculty  at  the  Maria 
Curie-Sklodowska  University  in  Lublin,  whilst  the 
LAICP  was  make  in  Environment  Analytic  Laboratory 
in  Geochemistry  and  Environment  Protection  Faculty 
at  the  Jan  Kochanowski  University  in  Kielce,  Poland. 
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Sample  places 


Fig.  1    Geology  map  of  the  Khibina  Mts.  with  the  samples  localization. 


3.  Results 

3. 1  Microscopic  Observations 

In  the  occur  of  the  veins,  rocks  are  present  massive 
syenites  with  a  nepheline-aegirine  urtites  and 
apatite-nepheline  ores.  Massive  syenites  usually  have 
large  crystalline  structure,  dense,  disordered  texture, 
rare  radiant  and  directional.  Background  of  massive 
syenites  form  the  automorphic  nepheline  crystals  (up 
to  25%  vol.),  orthoclase  (up  to  24%  vol.)  and  apatite 
(up  to  15%  vol.)  among  which  are  the  needle  of 
aegirine-acmite  (up  to  21%  vol.),  and  xenomorphic 
crystals  of  eudialite  (often  with  numerous  inclusions 
these  last  crystals-to  15%  vol.),  plagioclases  (acid 
varieties  to  10%  vol.).  In  these  rocks,  they  are  also 
subordinate  to  a  maximum  of  15%  vol.:  ilmenite, 
titanite,  magnetite,  arfvedsonite,  riebeckite, 
astrophyllite,  aenigmatite,  lorenzenite  and  others  (Fig.  2). 


Nepheline-aegirine  urtites  are  gery-black-grennish 
rocks  with  the  visible  on  the  background  a  nepheline 
crisstals  with  accompanied  plagioclases,  aegirines 
acmites  with  accesoric  minerals.  It  has  a 
large-crystalls,  pseudotrahyte  structure.  The 
apatite-nepheline  ores  it  are  a  green-brownish  color 
(with  red  tint).  There  are  massive  rocks  with  nepheline 
crystalls  and  tytanitem  with  addition  of  apatite.  In 
these  massive  syenites  is  present  a  numerous  of 
mineral  veins  (Table  1)  such  as:  tinguaites,  melteigites 
and  microsyenites. 

Syenite  veins  have  a  massive  texture,  fine 
crystalline  structures.  These  veins  are  the  largest 
group  of  rocks.  Form  the  veins  of  both  a  sharp  contact 
and  (often)  those  which  have  blurred  boundaries  are 
plastically  deformed,  with  budinage  structures.  They 
are  usually  microsyenites  (ijolites)  with  an  increased 
amount  of  mafic  minerals  (up  to  70%  vol.),  are  rich  in 
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v  *•* 


(a)  (b) 
Khibinite  rock  ((a)  thin  section,  crossed  pollars;  (b)  the  polished  sample  of  rocks) 

t 


(c)  (d) 
Urtite  rock  ((c)  thin  section,  crossed  pollars;  (d)  the  polished  sample  of  rocks) 


(e)  (f) 
Titatnite-nepheline-apatite  ore  mineralization  ((e)  thin  section,  crossed  pollars;  (f)  the  polished  sample  of  rocks) 


(g)  (h) 
Microsyenite  rocks  (thin  section,  crossed  pollars,  sample:  (g)  14CH02,  (h)  19CH02) 


(0  © 

Tinguaite  rocks  (thin  section,  crossed  pollars,  sample:  (i)  23CH02,  (j)  02CH10) 


(k)  (1) 
Melteigite  rocks  (thin  section,  crossed  pollars,  sample:  (k)  08CH02,  (1)  02CH03) 
Fig.  2    Microphotographs  of  typical  rocks  from  studied  area. 
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Table  1    Planimetric  composition  of  studiem  rocks. 
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ore  minerals,  mainly  from  the  group  of  oxides  of  iron 
and  titanium  with  admixtures  of  sulfides  (including 
copper).  Such  melanocratic  variety  of  syenite  more 
often  found  in  repeated  vein  fillings.  The  background 
are  the  crystalls  of  orthoclase  rarely  microclines 
(under  to  55%  vol.),  and  aegirines  (under  53%  vol.) 
with  arfvedsonite  (24%),  ith  addition  of  apatite  (24%), 
nepheline  (39%),  eudialite  (do  48%  obj),  acid 
plagioclases  (Ab,  Lb,  9%  vol.  of  rocks),  and  ilmenite, 
titanite  with  sulphides  (pyrite  and  chalkopyrite).  The 
apatite  is  rich  with  strontium-calcium  compounent.  Near 
these  minerals  is  present  the  inclusions  of  fluorite. 

Tinguaites  hale  a  beautifull  grenish  color, 
nematoblastic,  fibroblastic  structure,  radial,  dense 
texture.  Matrix  consists  of  radial  concentrations  of 
aegirine  and  astrophillite  (up  to  45%  vol.).  Near  these 
minerals  are  orthoclase  (22%),  nepheline  (5%),  and 
titanite  (up  to  3%  vol.).  The  rocks  have  also  magnetite 
and  ilmenite  (up  to  4%  vol.).  These  minerals  form 
narrow  blades  supporting  aegirynes.  These  veins  are 
most  frequently  met  in  the  central  part  of  Khibina.  In 


the  background  of  these  rocks,  they  are  present  a 
radial  conglomerates  of  aegirine  with  astrophyllites 
(under  45%  vol.). 

Melteigites  there  are  veins  have  sharp  contacts  with 
host  rocks.  They  have  lepido-nemato-crystallic 
structures,  massive,  poikilite,  porphyritic  textures. 
Phenocrystalls  of  minerals  are  augite  with  addition  of 
diopsides  (to  24%  vol.),  and  nodules,  usually  made  up 
of  olivine  (rich  in  Fe,  Fig.  4)  with  rims  of  Fe-biotite 
and  phlogopite  (Fig.  3).  Matrix  is  composed  of  apatite 
(up  to  10%  vol.),  nepheline  (up  to  15%  vol.),  orthoclase 
(up  to  15%  vol.)  and  plagioclase  (up  to  9%  vol.).  Less 
frequently  rocks  are  pyroxenites  occurring  in  what  are 
called  pyroxene  (augite)  suns  (79%  vol.).  In 
melteigites,  pyroxene  is  sometimes  characterized  by 
structures  such  as  radial  augite  sun.  The  presents 
pyrites  have  a  small  additin  of  Ni  (0.5  wt.%)  and  Cd 
(0.2  wt.%). 

3.2  LA-CP  Analysis 

In  the  Tables  2  and  3,  it  is  resent  results  of  analysis 
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18ach02(3) 
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(a)  The  fotograph  after  La  ICP  analysis  on  theilmenite  and  augite. 
Fig.  3    Microphotographs  from  electron  microscopy  (SE)  of  samples  rock. 


(b)  Olywine  with  micas  corona. 


MglFe 


(a)  Projection  of  olywine  from  the  melteigite  rocks. 
Fig.  4    EDS  analysis  of  selected  minerals. 


(b)  Pyroxene  projection  from  the  analysed  samples. 


Table  2    Chemical  composition  of  rocks  used  by  the  LA-ICP-MS  (reference  material:  polymetalic  sulphide  MASS-1,  USGS). 


Sample 


Pb 


Cu 


As 


Zn 


Cd 


Ni 


Co 


Ag 


mg/kg  (ppm) 


15,508 

3.3 

24 

208 

23 

84 

11 

<  1 

15,827 

31 

35 

101 

30 

30 

6 

<  1 

21CH02 

15,744 

33 

2 

233 

25 

10 

5 

<  1 

15,867 

37 

32 

218 

24 

231 

25 

<  1 

15,888 

345 

25 

40 

24 

34 

17 

<  1 

15,919 

36 

35 

215 

24 

8 

1 

<  1 

08CH02 

15,895 

32 

36 

234 

24 

8 

1 

<  1 

15,899 

38 

36 

243 

25 

8 

1 

<  1 

<  1 

2 

<  1 

<  1 

<  1 

2 

<  1 

<  1 

14CH02 

<1 

<  1 

<  1 

24 

<  1 

2 

<  1 

<  1 

<  1 

2 

<  1 

30 

<  1 

2 

3 

<  1 

<  1 

2 

<  1 

8 

<  1 

3 

<  1 

<  1 

18CH02 

<  1 

1 

<  1 

20 

<  1 

7 

7 

<  1 
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Table  3  Chemical  composition  of  rocks  used  by  the  LA-ICP-MS  (reference  material:  REE  element  sof  the  glass  matrix 
NIST  612). 


Sample 

r  D 

Mi 
INI 

Ti 
1  1 

J_.ll 

OQ 

mg/kg 

;  (ppm) 

3 

1    /IT  f\ 

l,4zU 

1  1  m 
1,1U3 

i  ^ 
Id 

i 
1 

Z 

Z6 

16CH02 

<  1 

<  i 

26 

o 
0 

37 

1  c 

lo 

7 

<  1 

<  1 

<  1 

199 
V!),ZU  / 

13 
1 

25 

4 

771 
13 1 

<  i 

lo 

is> 

Id 

<  1 

/I  A  1  A 

4,434 

3 

n 

y 

n  1 

y  i 

24 

11 

J 

?R  1R1 

11 

/  o 

15CHCH 

23 

57 

61 

72 

<  l 

42,221 

3 

7 

269 

48 

109 

111 

178 

5 

52,081 

<  1 

3 

103 

3 

<  1 

<  1 

<  1 

<  1 

246 

<  1 

<  1 

1 

3 

2 

2 

<  1 

<  1 

240 

<  1 

<  1 

4 

23CH02 

2 
3 

312 
36 

696 
9 

2,341 
23 

2 

<  1 

871 
640 

<  1 
1 

<  1 
1 

6 
48 

4 

2 

<  1 

1 

1 

214 

<  1 

<  1 

<  1 

36 

7 

4 

<  1 

153 

<  1 

<  1 

4 

selected  rock  samples  from  the  study  area.  In  the  veins 
samples  is  present  a  higher  content  of  Pb  and  Zn 
(samples  21CH02  and  08CH02).  In  samples:  16CH02 
and  15CH03  was  measured  a  high  contents  of  Ti,  Ni 
and  Ce. 

4.  Conclusions 

Describe  veins  rocks  from  Khibina  have  advanced 
syenitization  processes.  The  primitive  mineral 
association  has  a  corrosion  in  rocks  veins.  This  was  a 
product  of  mixed  chemical  composition  of  veins 
materials  with  occur  rocks.  The  melteigite  rock 
probably  is  a  differenciation  product  in  the  intrusive 
magma  and  migration  of  basic  derivates  with 
contamination  of  syenite  fluids.  These  fluids  make 
unstable  the  ol-opx  association  (Figs.  3  and  4).  In  next 
step,  a  cpx  minerals  was  crystallized  which  was  make 
aegirinization  (secondary  processes).  It  can  provide 
the  numerous  structures  disintegrating,  crowns 
reaction  in  olivine  and  pyroxene  around,  and  relatively 
high  content  in  the  background  a  nepheline  and 
feldspar.  In  addition,  vein  rocks  fill  the  various 
discontinuities  accompanying  the  transition  zone 
which  is  undoubtedly  limit  ore  occurrences  and 
khibinites.  This  causes  multistage  migration  of  fluids 
in  zones  splits  and  enrichment  rock  core  in  the  phase 
associated  with  the  crystallization  of  residual  and 


after-magmas  processes  [4].  Hence,  are  visible  in  the 
area  of  various  multiple  fill  the  same  vein  different 
generation  and  the  presence  of  carbonate  minerals 
found  in  many  rocks  in  Khibina  [3,  9,  10]  and  the 
enrichment  of  these  rocks  in  the  elements  cerium 
group.  In  addition,  the  activity  associated  with 
hydrothermal  after-magmas  processes  involves  the 
enrichment  of  these  rocks  Zn-Pb  sulphides  and 
minerals,  for  example,  fluorite. 
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Abstract:  The  aim  of  this  work  is  to  determine  the  porosity  of  concrete  by  means  of  two  NDT  (non-destructive  testing)  using  digital 
images.  In  one  test,  the  images  were  obtained  through  X-ray  microtomography  and  for  another  test  via  digital  scanner.  In  both  of  the 
tests,  the  images  were  processed  using  techniques  of  mathematical  morphology  and  pixel  processing.  For  both  NDT,  it  was  used 
concrete  samples  with  20-30  MPa  of  compressive  strength.  The  results  for  the  porosity  were  compared  with  results  obtained  by  the 
standard  test  proposed  by  NBR  (Norma  Brasileira)  9778  (2005)  showing  the  compatibility  between  two  studied  methods  and  pattern 
method. 

Key  words:  Concrete  porosity,  X-ray  microtomographic,  images  processing,  matematical  morfology. 


1.  Introduction 

The  porosity  of  a  material  can  be  thought  as  a 
property  to  present  pores  or  voids  and  is  symbolized 
by  the  fraction  of  the  total  mass  of  a  sample  that  is 
occupied  by  pores  or  voids.  In  the  concrete  matrix,  the 
porosity  has  a  great  influence  on  their  mechanical 
strength,  workability  and  durability.  By  the  previous 
becomes  important  to  develop  methods  to  regulate  not 
only  the  index  of  voids,  but  besides  the  physical 
characteristics  and  placement  of  the  pores. 

Different  methods  for  determining  the  porosity  of 
concrete  does  not  always  result  in  similar  values.  That 
is  because  traditional  testing  methods  for  determining 
this  index  usually  change  the  structure  of  the  cement 
paste,  especially  when  they  involve  the  addition  or 
emulsion  of  liquid  [1].  An  example  is  the  use  of 
mercury  intrusion  to  study  the  pore  system  of  the 
cement  paste,  admitting  that  the  pores  will  be 
narrowing  with  increasing  depth,  whereas  in  fact, 
some  pores  have  a  narrower  aperture,  thus  distorting 
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the  value  of  the  porosity  determined  by  this  method  [2]. 
Systematization  of  methods  to  determine  the  porosity 
of  some  materials  has  been  the  object  of  the  works 
being  produced  over  time.  Some  of  these  methods  use 
high-quality  images  and  work  with  techniques  and 
programs  of  digital  image  processing. 

In  this  study,  two  NDT  (non-destructive  testing)  to 
obtain  the  porosity  of  concrete  were  used.  In  the  first, 
NDT  X-ray  microtomography  was  used  with  the 
purpose  of  being  acquired  digital  images  of  the 
internal  structure  of  the  test-bodies.  In  the  second, 
NDT  a  commercial  scanner  was  used  to  directly  scan 
the  surfaces  of  cross  sections  of  these  tests-bodies.  For 
both  tests,  the  images  were  digitally  processed  to 
measure  the  porosity  of  the  material.  This  digital 
image  processing  consisted  mainly  in  performing 
mathematical  morphology  operations  and  processing 
of  pixels. 

The  average  surface  porosity  found  in  the  concrete 
were  compared  with  the  volumetric  porosity  detected 
using  the  Brazilian  standard  NBR  (Norma  Brasileira) 
9778  (2005)  [3].  This  method  is  a  pattern  for  this  kind 
of  measures;  tests  with  the  standard  method  were 
performed  in  the  Department  of  Materials  Engineering, 
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Federal  University  of  Sao  Carlos  SP.  In  Section  2,  we 
are  presented  some  basics  related  to  image  processing 
and  mathematical  morphology;  the  Section  3  showing 
the  experimental  procedure  and  commented  the 
performed  tests  and  the  Section  4  presented  the 
obtained  results  and  conclusions. 

2.  Images  Processing 

Understanding  an  image  as  a  compact  way  to 
represent  a  lot  of  information,  in  a  system  of  image 
processing  this  information  may  go  through  diverse 
forms  of  representation  depending  on  the  purpose 
defined  for  the  application.  The  first  step  of  the  system 
is  the  acquisition  of  an  image,  where  one  type  of 
sensor  can  be  used  for  digitizing.  The  sensor  may  be  a 
scanner,  a  TV  camera,  a  photo  camera,  a  CT  scan,  etc.. 
The  nature  of  the  sensor  and  the  image  it  producing  is 
determined  by  the  application. 

In  image  processing,  morphological  study  focuses 
on  the  geometric  characteristics  of  the  images  and 
applies  in  many  areas  of  image  processing,  such  as 
filtering,  segmentation,  enhancement  and  others  [4]. 

The  main  idea  of  mathematical  morphology  consists 
of  sets  of  picture  elements  (pixels),  which  are  placed  in 
groups,  representing  a  2D  structure  (shape),  and  over 
which  mathematical  operations  are  carried  to  highlight 
specific  aspects  of  the  forms  so  that  they  can  be 
recognized  and  measures  [5]. 

A  mathematical  morphology  can  be  translated  as 
the  application  of  mathematical  definitions  and 
principles  emphasizing  specific  aspects  of  the  forms 
allowing  them  to  be  counted  or  recognized.  It  is  a 
methodology  for  the  description  of  operators,  and  its 
base  is  to  extract  from  unknown  image  its  geometry 
and  topology  by  using  the  transformation  of  another 
completely  defined  image,  called  structuring  element 
image.  The  main  applications  of  mathematical 
morphology  are  in  the  process  of  noise  filtering, 
segmentation  parts,  skeletization,  etc.  [6]. 

3.  Experimental  Procedure 

The  test-bodies  were  molded  in  concrete  with  20  MPa 


and  concrete  with  30  MPa  compressive  strength,  the 
materials  were  collected  in  a  cement  truck.  Table  1 
shows  the  dimensions  of  cylindrical  test-bodies  used 
in  each  testing. 

The  bodies  were  originally  molded  with  200  mm  x 
100  mm,  and  which  were  cut  to  build  the  smaller 
bodies. 

3.1  Volumetric  Porosity  by  NBR  9778 

The  procedures  for  performing  this  test  are 
described  in  detail  in  the  Brazilian  standard  NBR  778 
(2005).  The  method  is  based  on  the  technique  of 
mercury  porosimetry.  This  technique  consists 
essentially  in  measuring  the  extent  of  mercury 
penetration  into  an  evacuated  solid  as  a  function  of  the 
applied  hydrostatic  pressure.  The  use  of  these 
techniques  requires  specialized  equipment  and  the 
proper  handling  of  mercury  for  it  to  be  a  contaminant 
element.  Table  2  shows  the  values  of  volumetric 
porosity  obtained  in  the  test  for  each  of  the  studied 
samples. 

3.2  Surface  Porosity  by  Microtomography 

To  accomplish,  the  X-ray  microtomography  was 
used  a  compact  instrument  for  high-resolution 
microtomography  (14.97  um)  from  Sky  Scan,  model 
1172.  This  system  has  a  variable  size  of  focus 
depending  on  the  power  applied.  Table  3  summarizes 
the  main  characteristics  of  this  equipment. 

In  order  to  minimize  distorting  effects  of  CT  due  to 
the  beam  hardening  effect  [7],  two  filters  (Al  +  Cu) 
were  used,  placed  between  the  radiation  source  and 
the  object  to  be  inspected.  Fig.  1  shows  the  CT 


Table  1    Dimensions  for  cylindrical  test-bodies. 


Method 

Diameter  (mm)         Height  (mm) 

NBR  9778 
Micro  CT 
Scanner 

100  200 
25  50 
100  15 

Table  2    Values  of  volumetric  porosity  with  NBR  778. 

Concrete 

Volumetric  porosity  (%) 

20  MPa 
30  MPa 

1.9 
2.1 
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Table  3    Main  features  of  the  Sky  Scan  1172. 


Device 


Characteristics 


X-ray  source 
X-ray  detector 

Maximum  size  of  the  object 

Details  of  detectability 
Reconstruction 
Radiological  safety 


20-100  kV,  10W,  <5  urn 

20-80  kV,  8  W,  <  8  urn 

11   Mp,  CCD  (charge  couple 

device)  of  12-bits  with  optical 

coupling  to  a  scintillator 

27    mm    in    diameter  (single 

sweep) 

<  0.8  urn  at  the  highest 
resolution 

Computer  or  cluster 

<  1  |a,Sv/h  at  any  point  on  the 
equipment  surface 


Fig.  1    CT  scanner  Sky  Scan  1172. 

scanner  used  in  the  research.  Due  to  the  CT  scanner 
used  can  only  measure  samples  of  concrete  with  small 
diameter,  we  used  7  cylindrical  test-bodies  whose 
dimensions  were  specified  in  Table  1 . 

A  program  was  developed  in  MATLAB  to  read 
each  generated  image  and  determine  the  area  of  voids 
to  obtain  the  surface  porosity.  In  this  program,  only 
two  parameters  required  user  intervention  the  contrast 
intensity  and  the  binarization  threshold.  These 
parameters  were  modified  for  each  set  of  processed 
images.  For  each  body-test,  an  average  of  1,000 
images  were  obtained.  Each  image  was  adjusted  in 
contrast,  binarized  and  filtered  using  morphological 
closing  operation  (dilatation  followed  by  erosion)  with 
circular  structuring  element  with  2  pixels  of  radius  to 
remove  isolated  points  interpreted  as  noise.  In  the  2D 
result,  the  empty  spaces  images  were  counted.  The  Fig.  2 
showed  the  sequence  of  operations  of  the  developed 
program. 


Read  Image 
1 


Adjusts  grayscale  image 


Image  Binarization 


Image  Inversion 
I 


Morphological  Closing  Operation 


Sum  White  Pixels  (pores) 
I 


Compute  percentage  of  white  pixels  to  total 
pixels  of  the  image,  calculates  the  total  area  of 

the  sample,  concrete  area,  area  of  voids 
(porosity),  the  surface  porosity  ratio  (%)  and 
pore-size  distribution 


Fig.  2  Flow  diagram  of  software  to  process  tomography 
images. 

The  Eq.  (1)  was  used  to  compute  the  surface 
porosity  of  concrete. 

AT-AC 


P(%) 


*  100 


(1) 


where,  P  represented  the  superficial  porosity  and  AT, 
Ac  are  total  area  and  concrete  area,  respectively. 

Fig.  3  shows  a  cross-section  with  two  cuts  from  an 
image  reconstructed  from  radiographs  obtained. 

Fig.  4  shows  the  processing  stages  of  reconstructed 
slice  acquired  in  microtomography  process. 

Tables  4  and  5  show  the  results  achieved  with  this 
method  for  the  samples  tested  with  compressive 
strength  of  20  MPa  and  30  MPa,  respectively.  The 


r  Irwwrse  [irey 
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Fig.  3    Two  cuts  from  a  reconstructed  images. 
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Table  6    Porosity  values  obtained  with  microCT. 


Fig.  4  Processing  stages  of  tomography  images:  (a) 
Original  image;  (b)  Ajusted  image;  (c)  Binarized  image;  (d) 
Inverted  image  and  (e)  Filtered  image. 

Table  4    Surface  porosity  for  bodies  of  20  Mpa  concrete. 


Test-body  Images 


Surface  porosity 
(%) 


Standard  deviation 

(%) 


1 

1,045 

1.1 

0.57 

2 

1,044 

0.8 

0.42 

3 

1,044 

1.4 

0.97 

4 

1,044 

0.6 

0.52 

5 

1,046 

1.2 

1.27 

6 

1,046 

0.9 

0.77 

7 

1,046 

1.4 

1.06 

Table  5    Surface  porosity  for  bodies  of  30  Mpa  concrete. 


Test-body  Images 


Surface  porosity 
(%) 


Standard  deviation 
(%) 


1 

1,046 

1.0 

0.80 

2 

1,046 

1.2 

0.55 

3 

1,046 

2.3 

1.62 

4 

962 

0.9 

0.57 

5 

1,046 

2.5 

1.26 

6 

1,046 

1.3 

0.67 

7 

1,046 

2.2 

0.95 

Concrete      Number  of  images 
20MPa  7,315 
30  MPa  7,238 


Volumetric  porosity  (%) 

1.0 

1.6 


Tables  4  and  5  columns  represent  the  sample  number, 
the  number  of  images  analyzed,  the  surface  porosity 
and  the  standard  deviation  between  images  of  the 
same  test-body. 

The  standard  deviation  of  about  1%  of  all  samples 
studied  indicate  the  similarity  of  the  porosity  values 
calculated  for  each  slice  of  the  reconstruction 
belonging  the  same  test-body.  Table  6  shows  the  total 
mean  surface  porosity  for  the  two  series  of  test-bodies. 

3.3  Surface  Porosity  by  Scanner 

Forty-two  circular  surfaces  measuring  100  mm  in 
diameter  acquired  from  cylindrical  test-bodies  of  100  mm 
diameter  by  1 5  mm  in  height  were  used  for  each  type 
of  concrete.  On  each  surface  were  acquired  four 
square  images  measuring  25  mm  side,  totaling  168 
images  to  be  analyzed  for  each  particular  type  of 
concrete.  A  shield  made  with  transparent  plastic  with 
openings  in  the  reading  area  was  used  so  the  concrete 
surface,  which  did  not  damage  the  scanner  glass.  Fig. 
5  shows  an  example  of  the  concrete  surface  to  be 
scanned. 


Fig.  5    Concrete  surface  to  be  scanned. 
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The  images  were  scanned  using  a  multifunctional 
EPSON  TX200  and  digitized  with  2400  dpi  of 
resolution.  To  crop  the  interest  areas,  the  tools  available 
on  a  typical  software  was  used  to  manipulate  images 
how  to  Paint  or  Photoshop.  A  program  was  developed 
in  MATLAB  to  read  the  cropped  images  and  process 
them  in  order  to  determine  which  areas  are  empty.  Each 
image  was  read  and  subsequently  separated  into  its 
HSV  (Hue,  Saturation  and  Value  mode  color) 
components,  since  the  images  obtained  by  the  scanner 
are  in  RGB  (Red,  Green  and  Blue  color  mode)  and  not 
grayscale  as  in  the  case  of  images  obtained  by 
microtomography.  HSV  was  chosen  over  the  standard 
RGB  because  this  was  not  possible  to  clearly  identify 
the  pores  in  any  of  its  layers,  R,  G  or  B.  Making  use  of 
the  HSV  standard,  the  layer  V  was  chosen  for 
processing;  due  to  this  layer  had  the  best  representation 
of  pores  in  the  image.  Each  image  was  converted  to 
HSV  pattern,  extracted  V  layer;  then  V  layer  was 
adjusted  in  contrast,  binarized  and  filtered  using 
morphological  closing  operation.  In  the  result  image, 
the  empty  spaces  were  counted.  The  Fig.  6  showed  the 
sequence  of  operations  of  the  developed  program  to 
process  the  cropped  images  acquired  with  the  scanner. 


Read  RGB  Image 


T 


Convert  Image  to  HSV 


Extract  V  Layer 


Adjusts  V  Layer  Contrast 


I 


Binarize  V  Layer 
I 


Invert  V  Layer 


I 


Morphological  Closing  Operation 
I 


Sum  White  Pixels  (pores) 


I 


Compute  percentage  of  white  pixels 
to  total  pixels  of  the  image 


In  this  program,  it  also  the  parameters  contrast 
intensity  and  binarization  threshold  requiring  user 
intervention.  These  parameters  were  adjusted  for  each 
group  of  processing  images.  Fig.  7  shows  the 
processing  stages  of  the  image  acquired  cutting 
scanned  image. 

With  the  same  structure  of  Tables  4  and  5,  Tables  7 
and  8  shows  the  results  achieved  with  this  method  for 
both  types  of  concrete. 
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Fig.  6    Flow  diagram  of  software  to  process  scanned  images. 


(g)  (h) 
Fig.  7    Processing  stages  of  scanned  images:  (a)  Original 
image;  (b)  H  layer;  (c)  S  layer;  (d)  V  layer;  (e)  Ajusted  image; 
(f)  Binarized  image;  (g)  Inverted  image  and  (h)  Filtered 
image. 
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Table  7 

Surface  porosity  for  bodies  of  20  Mpa  concrete. 

Test-body  Images 

Surface  porosity 
(%) 

Standard  deviation 

(%) 

1 

24 

1.08 

0.53 

2 

24 

1.67 

0.82 

3 

24 

1.46 

0.45 

4 

24 

2.28 

1.46 

5 

24 

2.03 

1.81 

6 

24 

1.63 

0.57 

7 

24 

1.79 

0.73 

Table  8 

Surface  porosity  for  bodies  of  30  Mpa  concrete. 

Test-body  Images 

Surface  porosity 
(%) 

Standard  deviation 
(%) 

1 

24 

1.04 

0.48 

2 

24 

1.15 

0.56 

3 

24 

1.65 

0.58 

4 

24 

1.98 

0.71 

5 

24 

1.58 

0.53 

6 

24 

1.64 

0.53 

7 

24 

1.92 

0.71 

Table  9  shows  the  total  mean  surface  porosity  for 
the  two  series  of  test-bodies. 

4.  Analysis  of  Results  and  Conclussions 

Table  10  summarizes  a  comparison  between  the 
results  obtained  by  the  two  methods  discussed  in  this 
work  and  the  results  of  the  test  prescribed  in  NBR 
9778. 

Considering  that,  the  linear  surface  porosity  and 
volumetric  porosity  must  be  equivalent.  The  average 
surface  porosity  found  using  microtomography  and 
scanner  can  be  compared  with  the  volume  porosity 
found  in  the  samples  of  concrete,  prepared  with  the 
same  batch  and  tested  by  the  method  shown  in  NBR 
9778.  Table  10  shows  some  differences  between  the 
values  found  for  the  concrete  porosity  by  different 
methods;  these  possible  differences  were  also 
commented  by  other  authors  [7]. 

Analyzing  the  two  methods  used  to  measure  the 
concrete  surface  porosity,  similar  results  were  found 
in  the  30  MPa  concrete  and  difference  values  for  the 
20  MPa  concrete.  These  differences  could  be  due  to 
the  selection  of  parameters  contrast  intensity  and 
binarization  by  the  program  user.  With  this  selection, 


Table  9 

Porosity  values  obtained  with  scanner. 

Concrete 

Number  of  images     Volumetric  porosity  (%) 

20  MPa 
30  MPa 

7,315  1.0 
7,238  1.6 

Table  10 

Results  summary  to  concrete  porosity. 

Concrete 

.               r             ..      Volumetric  porosity 
Average  surface  porosity  .... 

(  /a) 

MicroCT       Scanner       NBR  9778 

20  MPa 
30  MPa 

1.7                1.0  1.9 
1.6                1.6  2.1 

what  may  happen  is  that  some  pores  are  lost  or  some 
noise  can  be  considered  as  pore. 

It  can  be  seen  that  the  porosity  found  for  20  MPa 
concrete  by  the  method  was  shown  in  NBR  9778  and 
found  by  analysis  of  images  obtained  by 
microtomgraphy  is  lower  than  that  found  in  the  30  MPa 
concrete.  This  may  be  due  to  changes  in  the  water 
cement  ratio  in  the  feeder  of  cement  truck  change  the 
voids  in  concrete  with  higher  strength,  since  that  the 
technical  literature  reports  higher  values  of  porosity 
for  concrete  with  lower  resistance.  In  general,  to 
analyze  the  compressive  strength  of  the  concrete 
should  also  be  considered  the  size  and  distribution  of 
pores;  in  our  samples  the  30  MPa  concrete  has  a 
higher  amount  of  smaller  pores. 

The  use  of  imaging  by  microtomography  or  the  use 
of  a  scanner  (in  both  cases  complemented  with 
techniques  of  image  processing)  is  shown  as  viable 
methods  to  estimate  the  concrete  porosity.  The 
method  using  the  scanner  is  presented  as  an  acceptable 
method  for  their  results,  and  especially  for  its  low  cost 
and  ease  operation. 
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Abstract:  Nanosized  fibrous  cerium(IV)  hydrogen  phosphate  membrane,  Ce(HP04)2'2.9H20  (nCePf),  was  prepared  and 
characterized  by  chemical,  XRD  (X-ray  diffraction),  TGA  (thermogravimetric  analysis),  SEM  (scanning  electron  microscopy)  and 
TEM  (transmission  electron  microscopy).  Novel  supported  nanofibrous  Ce(IV)  phosphate/polyiondole  nanocomposite  membranes 
were  prepared  via  in-situ  chemical  oxidation  of  the  monomer  that  was  promoted  by  the  reduction  of  Ce(IV)  ions  present  in  the 
inorganic  matrix.  The  presence  of  Ce(IV)  ions  allows  redox  reactions  necessary  to  oxidative  polymerization  to  occur.  The  resultant 
material  was  characterized  by  TGA,  elemental  (C,  H,  N)  analysis  and  FT-IR  (Fourier  transform  spectroscopy).  SEM  images  of  the 
resulting  nanocomposite  reveals  a  uniform  distribution  of  the  polymer  on  the  inorganic  matrix.  Amount  of  polyindole  polymer 
present  in  the  composite  is  found  to  be  ~  7.0%. 

Key  words:  Nanofibrous  cerium(IV)  hydrogen  phosphate  membrane,  self  support  polymerization,  indole. 


1.  Introduction 

In  recent  years,  heterocyclic  conducting  polymers 
containing  nitrogen  atoms  like  polyaniline  and 
polypyrrole,  and  their  substituted  derivatives  have 
received  increasing  attention  in  various  fields, 
electronics,  industry  and  others  [1-3].  Their  electrical 
and  electrochemical  properties  show  great  promise  for 
commercial  applications.  However,  among  various 
aromatic  compounds-based  conducting  polymers, 
polyindole  and  its  derivatives  has  been  less 
investigated  although  there  exists  close  structural 
similarities  with  the  polymers  mentioned  above  [4]. 
Polyindole  can  be  obtained  from  chemical, 
electrochemical  and  interfacial  polymerization  of 
indole  [2,  3,  5,  6].  Its  electrical  and  electrochemical 
properties  show  great  promise  for  commercial 
applications  [5-10].  Polyindole  is  an  electro  active 
polymer,  owns  advantages  especially  fairly  good 
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thermal  stability  [5,  6].  Some  studies  shows 
polyindole  has  similar  properties  like  polyaniline, 
based  on  their  high  conductance  and  good 
environmental  stability  [2,  4,  11-13].  Inorganic 
layered  nanomaterials  are  receiving  great  attention 
because  of  their  size,  structure,  and  possible 
biochemical  applications  [14,  15],  that  have  been 
proven  to  be  good  carriers  for  organic  polar  molecules. 
Examples  of  these  are  zirconium  phosphates, 
Zr(HP04)2-H20  (ZrP),  which  are  inorganic  cation 
exchange  materials  with  high  thermal  stability, 
solid-state  ion  conductivity,  resistance  to  ionizing 
radiation,  and  which  are  known  as  hosts  capable  of 
incorporating  different  types  and  sizes  of  guest 
molecules  [15]. 

Researchers  have  been  capable  of  encapsulating 
functional  biomolecules  into  these  inorganic  matrices 
protecting  them  from  interacting  with  environment, 
avoiding  denaturation  and  enhancing  their  shelf  life 
[14,  15].  Crystalline  cerium  phosphates  have  been 
studied  for  a  long  time  as  ion  exchangers,  their 
structures  remains  unknown  until  recently  [16-18]. 
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The  reason  is  that,  the  composition,  the  structure  and 
the  degree  of  crystallinity  of  their  precipitates  results 
from  reaction  of  solutions  containing  a  Ce(IV)  salt 
which  is  mixed  with  a  solution  of  phosphoric  acid  of 
[(PC>4)/Ce(IV)  ratio],  strongly  depend  on  the 
experimental  conditions  such  as  rate  and  order  of 
mixing  of  the  solutions,  stirring,  temperature  and 
digestion  time,  this  also  implemented  on  fibrous 
cerium  phosphate  [19].  To  date,  most  of  the  work  on 
fibrous  cerium  phosphate  was  carried  out  on  its  ion 
exchange  [20],  intercalation  [21]  and  electrical 
conductance  properties  [22].  Studies  on  its 
polystyrene,  polyacrylamide  [23]  and  (polyvinyl 
chloride-based  polyvinyl  (alcohol)  [24]  composites 
have  been  reported.  Nanoscaled  tetravalent  metal 
phosphates  and  their  organic  polymer  composites 
comprise  an  important  class  of  synthetic  engineering. 
However,  research  in  such  area  is  still  terra  incognita 
[25-28].  Nanotechnologies  are  at  the  center  of 
numerous  investigations  and  huge  investments. 
However,  chemistry  has  anticipated  for  long  the 
importance  decreasing  the  size  in  the  search  of  new 
properties  of  materials,  and  of  materials  structured  at 
the  nanosize  in  a  number  of  applications  relate  to  daily 
life.  Organic-inorganic  nanocomposite  membranes 
have  gained  great  attention  recently  [28,  29].  The 
composite  material  may  combine  the  advantage  of 
each  material,  for  instance,  flexibility,  processability 
of  polymers  and  the  selectivity  and  thermal  stability  of 
the  inorganic  filler  [26-30].  Conducting  polymers  are 
interesting  materials  owing  to  their  electrical 
properties  [25,  27].  Belonging  to  this  class  are 
polyaniline,  polypyrrole,  polybenzimidazole  and 
others.  Polyppyrole  polymerization  by  fibrous  ceium 
phosphate  has  been  reported  [31,  32]. 

In  our  laboratory,  we  are  carrying  systematic 
investigations  on  novel  tetravalent  metal 
phosphates/organic  heterocyclic  conducting  polymers 
nanocomposite  membranes.  Recently,  we  have 
reported  [33]  the  preparation  and  characterization  of 
fibrous        cerium  phosphate/polybenzimidazole 


nanocomposite  membrane.  The  present  study 
describes  the  preparation  and  characterization  of  novel 
supported  fibrous  cerium  phosphate/polyindole 
nanocomposite  membranes  via  in-situ  chemical 
oxidation  that  was  promoted  by  the  reduction  of 
Ce(IV)  ions  presenting  in  the  inorganic  matrix. 

2.  Experiments 

Chemicals  Ce(S04)-2.4H20,  H3P04  (85%)  of 
British  drug  house,  PVA  (Polyvinylealcohol),  MWt 
(molecular  weight)  =  25,125  g/mol  of  Aldrich.,  Indole 
of  Reidel  de-Haen,  were  used  as  such.  Other  reagents 
used  were  of  analytical  grade. 

2. 1  Instruments  Used  for  Characterization 

X-ray  powder  diffractometer  Siemens  D-500,  using 
Ni-filtered  CuKa  {X  =  1.54056  A);  TG/DTA 
(thermogravimetric  loss/differential  thermal  analysis), 
DTSIIExtra  6000  Thermogram  and  TG/DTA 
Perkin-Elmer  SII;  Fourier  Transform  IR  spectrometer, 
model  IFS  25  FT-IR,  Bruker;  SEM  (scanning  electron 
microscopy)  Jeol  SMJ  Sm  5610  LV;  TEM 
(transmission  electron  microscopy)  Zeiss  TEM  10CR 
andpH  Meter  WGW  521. 

2.2  Preparation  of  Nanofibrous  Cerium  Phosphate 
Membrane:  Ce(HP04)r2.9H20 

Nanofibrous  cerium  phosphate  membrane  was 
prepared  from  adding  150  mL  of  0.05  M  CeS04-4H20 
in  0.5  M  H2SO4  solution,  dropwise,  to  150  mL  of  6  M 
H3PO4  at  80  °C  with  stirring.  After  complete  addition, 
the  resultant  material  left  to  digest  at  that  temperature 
for  5  h.  To  that,  1  L  of  hot  distilled  water,  (~  60  °C), 
was  added  with  stirring  for  1  h.  The  resulting  fibrous 
cerium  phosphate  was  subjected  to  washing  by  with 
distilled  water  up  to  pH  3,  then  filtered,  washed  with 
ethanol  and  dried  in  air. 

2.3  Exchange  Capacity 

Exchange  capacity  of  nanosized  fibrous  cerium 
phosphate  membrane  was  determined  by  addition  of 
25  mL  of  0.1  M  NaCl  solution  to  100  mg  of  the 
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material,  with  stirring  for  1  h,  then  titrated  with  0.1  M 
NaOH  solution. 

2.4  Thermal  Analysis 

Thermal  analyses  were  carried  out  at  temperature 
range  about  20-775  °C  in  nitrogen  atmosphere,  the 
rate  was  1 0  °C/min. 

2.5  Preparation  of  Nanofibrous  Cerium  Phosphate 
Polyindole  Composite  Membranes 

Cerium  phosphate/polyindole  nanocompsite 
membranes  were  prepared  as  follows: 

(1)  By  immersion  of  fibrous  cerium  phosphate  self 
supported-sheet  160  mg  in  10  mL  ethanol  solution 
containing  200  mg  indole,  at  room  temperature  for  5 
days.  It  was  observed  with  time  the  color  changes 
gradually  (pale  brown,  brown  and  finally  to  dark 
brown).  The  impregnated  sheet  was  removed,  washed 
with  ethanol  and  left  to  dry  in  air. 

(2)  By  immersion  of  fibrous  cerium  phosphate  self 
supported-sheet  160  mg  in  10  mL  ethanol  solution 
containing  200  mg  indole,  at  room  temperature  for  2 
days.  It  was  observed  with  time  the  color  changes 
gradually  (pale  brown  and  finally  to  brown).  The 
impregnated  sheet  was  removed,  washed  with  ethanol 
and  left  to  dry  in  air. 

3.  Results  and  Discussion 

Nanofibrous  cerium  phosphate  membrane,  Ce 
(HP04)2-2.9H20  (nCePf),  was  prepared  and 
characterized  by  chemical,  XRD,  TGA,  FT-IR,  SEM 
and  TEM. 

XRD  of  nCePf  is  shown  in  Fig.  1,  with  dm  =  10.89  A 
Thermogram  of  nCePf  is  shown  in  Fig.  2.  The 
thermal  decomposition  occurs  in  continuous  process, 
the  thermal  analysis  was  carried  out  at  temperatures 
between  10-775  °C,  the  final  product  was  CeP207. 
Loss  of  water  of  hydration  occurs  between  60-200  °C, 
followed  by  POH  groups  condensation.  The  total 
weight  loss  found  to  be  equal  to  19.09%. 

Fig.    3    shows    FT-IR    spectrum    of  fibrous 


Ce(HP04)2-2.9H20,  with  a  trend  similar  to  that  of 
M(IV)  phosphates.  It  consists  of  broad  band  centered 
at  3,350  cm"1  is  due  to  OH  groups  symmetric 
stretching  of  H20,  small  sharp  band  at  1,628  cm"1  is 
related  to  H-O-H  bending,  and  sharp  broad  band 
centered  at  1,045  cm"1  is  corresponds  to  phosphate 
groups  vibration.  The  bands  at  the  region  630-450  cm"1 
are  ascribe  the  presence  of  §(P04). 

SEM  morphology  image  of  the  nanosized  fibrous 
cerium  phosphate  (nCePf)  is  shown  in  Fig.  4.  The 
photograph  shows  its  average  size  is  ~  20.5  nm. 

The  ion  exchange  capacity  of  nanosized  fibrous 
cerium(IV)  phosphate,  Ce(HP04)2-2.9H20,  found  to 
be  equal  to  5.21  meq/g. 

TEM  image  of  the  nanosized  fibrous  cerium 
phosphate  (nCePf)  is  shown  in  Fig.  5.  The  photograph 
shows  its  average  size  is  ~  15  nm. 

New  supported  polyindole  nanocomposite  was 
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Fig.  1    XRD  of  nanofibrous  cerium  phosphate. 


200.0       400.0        600.0  700.0 

Temperature  (°C) 
Fig.  2    TGA  of  nanofibrous  cerium  phosphate. 
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Fig.  3    FT-IR  spectra  of  nanofibrous  cerium  phosphate. 


Fig.  4    SEM  morphology  image  of  (nCePf). 


Fig.  5  Miograph  image  (TEM)  of  (nCePf)  composite  thin 
film  in  PVA  matrix. 

prepared  and  investigated  in  this  work.  Self-supported 
sheet  of  nanosized  fibrous  cerium(IV) 
hydrogenphosphate,  Ce(HP04)2-2.9H20  (nCePf),  was 
used  as  the  inorganic  matrix,  since  the  fibrous 


morphology  present  by  this  material  makes  molding 
possible,  giving  desired  shape  to  the  final 
(nCePf)/polyindole  naoncomposite  membrane.  It  was 
found  that  when  nanosized  fibrous  cerium(IV) 
hydrogen  phosphate  immersed  in  ethanolic  solution  of 
polyindole  the  color  of  self  supported  sheet  gradually 
changes  with  time  to  (pale  brown,  brown  and  finally 
to  dark  brown)  due  to  the  formation  of  nanocomposite 
membrane.  It  is  important  to  note  that  the 
self-suported  sheet  integrity  is  preserved.  This  result  is 
interesting  because  the  shape  integrity  makes  the 
building  of  molded  conducting  device  possible.  The 
resultant  material  was  characterized  by  elemental  (C, 
H,  N)  analysis  TGA  and  by  FT-IR  spectroscopy. 

From  elemental  (C,  H,  N)  analysis,  the  amount  of 
organic  material  present  in  the  (nCePf)  polyindole 
naoncomposite  was  6.81%  in  weight.  The  %  of 
elemental  (C,  H,  N)  related  to  the  organic  polymer 
found  to  be  C  =  5.58%,  H  =  0.41%  and  N  =  0.82%. 
From  TGA  analysis,  the  amount  of  the  organic 
polymer  found  to  be  7%  in  weight.  However,  on 
repeating  the  same  experiment  for  a  period  48h  time 
immersion  of  fibrous  cerium  phosphate  in  indole 
solution,  gave  low  yield  of  polyindole  (3.97%  in 
weight).  It  is  elemental  (C,  H,  N)  analysis,  were  C  = 
3.26%,  H  =  0.327%  and  N  =  0.475%.  So,  the  degree 
of  polymer  loading  can  be  controlled  by  varying  aging 
time. 

Fig.  6  shows  SEM  image  for  (nCePf)/polyindole 
naoncomposite  membrane,  reveal  a  uniform 
distribution  of  the  polymer  on  the  inorganic  matrix, 
coating  the  nCePf  and  filling  the  empty  spaces 
between  the  fibers.  Amount  of  polyindole  present  in 
the  composites  was  found  to  be  7%  in  weight. 

Fig.  7  shows  FT-IR  spectrum  of  (nCePf)/polyindole 
nanocomposite  membrane.  It  consists  of  broad  band 
centered  at  3,425  cm"1  is  due  to  OH  groups  symmetric 
stretching  of  H2O,  small  sharp  band  at  1,632  cm"1  is 
related  to  H-O-H  bending,  and  sharp  broad  band 
centered  at  1,055  cm"1  is  corresponds  to  phosphate 
groups  vibration.  Small  band  at  2,921  cm"1  corresponds 
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Fig.  6    SEM  image  of  (nCePf)/polyindole  nanocomposite. 

to  C-H  bonds,  bands  in  the  region  of  1,557-1,500  cm"1 
are  related  to  stretching  C-C  bonds  characteristic  of 
indole  unites,  bands  in  the  region  1,500-1,400  cm"1 
corresponds  to  C-N  bonds.  The  band  at  742  cm"1  is 
correspond  to  out  of  plane  deformation  of  C-H  bond 
in  the  benzene  ring  [34]. 

Thermogram  of  (nCePf)/polyindole  nanocomposite 
membrane  is  shown  in  Fig.  8,  the  weight  loss  up  to 
180  °C  is  7.35%  due  to  the  removal  of  external  water 
molecules.  The  weight  loss  occurs  between  180-700  °C 
corresponds  to  the  decomposition  of  polyindole  and 
condensation  of  P-OH  groups  of  the  inorganic 
material  to  pyrophosphate,  CeP207.  The  POH  groups 
condensation  found  to  cuoerimose  with  that  of 
polyindole.  The  total  weight  loss  found  to  equal  to 
18.562%.  From  thermal  and  elemental  (carbon, 
hydrogen,    nitrogen)  C,  H,  N  analysis,  the  resultant 
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Fig.     7        FT-IR    spectrum    of  (nCePf)/polyindole 
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Fig.  8    TG/DTA  of  (nCePf)/polyiondole  nanocomposite. 

product  was  formulated  as:  Ce(HPC>4)2-(Pmd)x-2.1 
H20.,  where,  x  =  ~  7.0%. 

4.  Conclusions 

Novel  supported  nanofibrous  cerium 
phosphate/polyindole  nanocomposite  membrane  was 
prepared  via  in-situ  chemical  oxidation  of  the  monomer 
that  was  promoted  by  the  reduction  of  Ce(IV)  ions 
present  in  the  inorganic  matrix.  The  presence  of  Ce(IV) 
ions  allows  redox  reactions  necessary  to  oxidative 
polymerization  to  occur.  The  self-supported  sheet 
integrity  is  preserved.  This  result  is  interesting 
because  the  shape  integrity  makes  the  building  of 
molded  conducting  device  possible.  The  formulation 
of  the  fibrous  cerium  phosphate  nanocomposite  was 
supported  by  thermal,  elemental  (C,  H,  N)  analysis, 
FT-IR  spectrum  and  SEM.  The  degree  of  polymer 
loading  can  be  controlled  by  varying  aging  time. 
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Abstract:  Terms  of  synthesis  were  determined  for  creation  of  new  generation  premixes  and  for  their  testing  in  experiments. 
Heteronuclear  chelate  citrates  of  general  formula:  MI2MI1L2'«H20  (where,  M1  =  Zn,  Co,  Fe,  Mn,  Cu;  Mn  =  Mn,  Zn,  Co,  Cu;  n  =  0/4) 
were  synthesized.  Identity  and  composition  of  synthesized  compounds  were  defined  by  microelemental  analysis,  determination  of 
melting  temperature  and  X-ray  diffraction  analysis.  X-ray  diffraction  method  was  used  also  to  define  crystallinity  of  the  compounds 
and  their  citric  acid  (H4L)  component.  X-ray  amorphous  and  iso-structural  orders  were  also  revealed. 

Key  words:  Bio-metals,  chelate  compounds,  microelemental  analysis,  X-ray  diffraction  analysis. 


1.  Introduction 

Even  the  most  primitive  forms  of  live  organisms  are 
unable  to  provide  vital  processes  without  participation 
of  bio-metals.  In  optimal  doses  they  are  able  to  affect 
actively  growth  and  development,  propagation, 
productivity,  resistance  of  organisms  to  diseases. 
Investigations  of  researchers  prove  that  bio-metals 
play  active  role  in  the  physiological  and  biological 
processes  going  on  in  poultry  and  animals.  Up  to  75 
chemical  elements  are  found  in  their  organisms  and 
namely  this  conditioned  the  interest  shown  to  the 
so-called  "vital  metals"  by  wide  specter  of  researchers. 
We  considered  five  microelements  (Zn,  Co,  Fe,  Mn, 
Cu),  because  these  microelements  belong  to  a  number 
of  very  critical  and  essential  for  life  elements  that 
animal  organisms  have.  They  play  an  important  role  in 
a  water  and  organic  compounds  metabolism, 
absorbtion  and  utilization  of  nutrients.  They 
generateoptimal  conditions  for  the  normal  function  of 
many  organs  including  cardiac  muscule  and  nervous 
system  [1-4].  Their  role  and  significance  are 
exemplified  by  the  fact  that  one  of  the  main  factors 
which  defines  low  quantity  and  quality  indices  of 


foodstuff  is  the  deficiency  of  bio-metals  in  plants,  soil, 
as  well  as  in  agricultural  animals  and  poultry. 
Decisive  role  in  the  resolution  of  this  problem  is 
attributed  to  the  provision  of  live  organisms  with 
optimal  quantity  and  ratio  of  microelements.  Since 
bio-metals  in  live  organisms  fulfill  their  functions  in 
the  form  of  chelates,  at  filling  up  of  deficit  in 
bio-metals  in  the  form  of  chelates,  we  observe  sharp 
increase  in  their  biological  activity  [5-8]. 

2.  Materials  and  Methods 

For  synthesis,  we  used  the  following  original 
compounds:  ZnO,  MnC03,  CoC03Co(OH)2, 
CuC03Cu(OH)2,  Fe(CH3COO)2-4H20,  and  important 
for  vital  processess  oxyacid-citric  acid  (H4L). 
Composition  of  compounds  and  their  identity  were 
defined  by  well-known  methods:  microelemental 
analysis  (CHN  Analyser,  Labertherm  GmbH).  Purity 
was  determined  by  melting  temperature  measuring 
(Dynalon  SMP10)  and  X-ray  diffraction  method  on 
DPOH-3Mtype  diffractometer,  at  copper  anode 
emission.  The  same  method  was  used  to  define  their 
crystallinity,  and  X-ray  amorphous  and  iso-structural 
orders  were  also  revealed. 
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3.  Results 

Investigations  are  continued  at  the  Laboratory  of 
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Problems  Agrarian  Chemistry  for  creation  of  new 
generation  premixes  and  for  their  testing  in 
experiments  [9-17].  With  this  in  view  heteronuclear 
chelate  citrates  of  the  following  general  formula: 

mVmu-L2-«H20  (where,  M1  =  Zn,  Co,  Fe,  Mn,  Cu; 
M11  =  Mn,  Zn,  Co,  Cu;  n  =  0/4)  were  created.  To 
obtain  these  citrates,  we  take  a  mix:  (1)  ZnO:MnbCC>3, 
(2)  CoC03Co(OH)2:ZnO,  (3)  CuC03Cu(OH)2: 
CoC03Co(OH)2,  (4)  CoC03Co(OH)2:MnC03,  (5) 
CoC03Co(OH)2:CuC03Cu(OH)2,  (6) 
Fe(CH3COO)2-4H20:MnC03,  (7) 
Fe(CH3COO)2-4H20:ZnO,  (8) 
Fe(CH3COO)2-4H20:CoC03Co(OH)2  and  (9) 
MnC03:CuC03Cu(OH)2  at  2:1  molar  ration  on 
porcelain  cup,  crush  intensely  till  it  is  reduced  to  fine 
dispersion  condition,  then  we  add  citric  acid  dissolved 
in  minimal  volume  of  hot  water.  We  observe  intense 
isolation  of  C02  and  changes  of  color.  The  obtained 
solutions  are  left  at  room  temperature  till  the  next  day, 
are  filtered  and  the  residue  after  its  washing  by  ether  is 


dried  on  air. 

Identity  and  composition  of  synthesized  chelate 
composites:  Zn2-Mn-L2-2H20,  Co2-Zn-L2, 

Cu2-Co-L2-H20,  Fe2-Mn-L2-2H20,  Co2-Mn-L2-4H20, 
Co2CuL2-4H20,  Mn2CuL2-4H20,  Fe2ZnL2-2H20  and 
Fe2CoL2-4H20  were  defined  by  microelemental 
analysis  (Table  1),  by  melting  temperature  detection 
(Table  2),  by  X-ray  diffraction  method  (Table  3). 

We  also  defined  their  crystallinity  and  revealed 
X-ray  amorphous  and  iso-structural  orders  by  X-ray 
diffraction  analysis. 

Table  3  represents  the  results  of  relative  intensities 
(I/Io),  reflection  angle  (20),  inter-plane  distances  (d) 
obtained  by  deciphering  of  X-ray  diffraction  patterns 
of  the  composites. 

4.  Discussions 

X-ray  diffraction  studies  did  not  reveal  presence  of 
starting  compounds  in  the  samples  in  the  form  of 
separate  phase.  Therefore,  we  can  conclude  that 


Table  1 

Results  of  microelemental  analysis  of  chelates. 

# 

Formula 

Practical  (%) 

Theoretical  (%) 

M1 

Mn 

C 

H 

M1 

M11 

C 

H 

1 

Zn2 

MnL2-4H20 

20.94 

9.05 

24.03 

2.15 

21.69 

9.17 

24.00 

2.34 

2 

Co2 

Zn-L2 

20.87 

10.72 

25.09 

2.17 

21.01 

11.58 

25.66 

2.49 

3 

Cu2 

CoL2H20 

22.11 

9.96 

24.45 

1.91 

21.81 

10.12 

24.73 

2.06 

4 

Co2 

•Mn-L2-4H20 

19.03 

8.73 

20.96 

2.72 

18.91 

8.83 

21.11 

2.89 

5 

Co2 

•CuL2-4H20 

18.21 

9.79 

22.75 

2.63 

18.66 

10.05 

22.79 

2.85 

6 

Fe2- 

MnL2-2H20 

19.33 

9.14 

24.62 

2.15 

19.27 

9.46 

24.77 

2.41 

7 

Fe2- 

ZnL2-2H20 

19.15 

11.13 

24.79 

2.23 

18.94 

10.99 

24.35 

2.37 

8 

Fe2- 

CoL2-4H20 

17.78 

8.93 

22.98 

2.65 

18.03 

9.49 

23.18 

2.89 

9 

Mn2-Cu-L2-4H20 

17.17 

10.09 

22.97 

2.77 

17.63 

10.18 

23.08 

2.89 

Table  2 

Results  of  measurements  of  composites  melting  temperatures. 

# 

Formula 

Melting,  T  (°C) 

Color 

1 

Zn2-Mn-  L2-4H20 

>286 

White 

2 

Co2ZnL2 

>286 

Lilac 

3 

Cu2CoL2H20 

>286 

Dark  light-blue 

4 

Co2MnL2-4H20 

>286 

Grey-lilac 

5 

Co2CuL2-4H20 

>286 

Dark  lilac 

6 

Fe2MnL2-2H20 

>286 

Mustard  yellow 

7 
8 

Fe2ZnL2-2H20 
Fe2CoL2-4H20 

>286 
>286 

Light  brown 
Grey-brown 

9 

Mn2-CuL2-4H20 

>286 

Light  green 
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Table  3    Results  of  X-ray  diffraction  analysis  of  chelate  composites  and  citric  acid. 
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absolutely  new  individual  compounds  are  obtained 
which  are  in  full  conformity  with  the  results  of  the 
above  given  micro  elemental  analysis  and  melting 
point  measurements  (Tables  1  and  2). 

X-ray  diffraction  patterns  of  the  composites 
Fe2-Mn-L2-H20  and  Fe2-Co-L2-4H20  are  characterized 
by  weakly  expressed  diffuse  peaks,  therefore  they 
represent  X-ray-amorphous  compounds.  Other 
compounds  are  characterized  by  clearly  expressed 
crystalline  structures,  various  relative  intensities  (//To), 


reflection  angles  (28),  and  definite  intra-spatial 
distance  (d)  values.  It  should  be  stated  that  X-ray 
diffraction  patterns  of  the  composites  Cu2-Co-L2-H20 
and  Co2-Cu-L2-4H20  are  identical,  which  speaks  of 
their  iso-structure  that  is  conditioned  by 
analogousness  of  their  chemical  composition  and 
geometrical  characteristics  of  crystalline  structure, 
which  can  not  be  said  about  the  composites: 
Zn2-Mn-L2-2H20,  Co2-Zn-L2,  Co2-Mn-L2-4H20, 
Mn2-Cu-L2-4H20,  Fe2-Zn-L2-2H20  (Table  3). 
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5.  Conclusions 

For  the  first  time,  we  synthesized  Heteronuclear 
chelate  citrates  of  general  formula:  M^M'^-nE^O  by 
using  modified  and  simplified  methods.  We  studied 
identity  and  composition  of  synthesized  compounds, 
determined  their  crystallinity,  and  revealed  X-ray 
amorphous  and  iso-structural  orders. 

We  have  conducted  preliminary  studies  on  these 
synthesized  compounds  as  new  poultry  premixes. 
Based  on  good  preliminary  data,  we  concluded  to 
conduct  large-scale  studies  in  this  direction. 
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Abstract:  Oxygen  plays  a  crucial  role  in  oenology.  The  oxygen  leads  to  changes  in  the  chemical  and  sensory  profile  of  wines.  On 
one  hand,  the  impact  of  moderate  oxygen  exposure  of  red  wine  has  a  positive  effect  in  its  colour,  aromatic  bouquet  and  mouth-feel 
properties.  On  the  other  hand,  oxygen  has  a  negative  effect  on  white  wine's  quality,  as  well  as  the  sensory  and  compositional  levels. 
The  purpose  of  this  study  was  to  quantify  the  dissolved  oxygen  in  red  and  white  wine,  during  different  cellar  operations  such  as 
racking,  tartaric  stabilization,  filtering  and  bottling.  The  techniques  that  contribute  most  to  the  enrichment  of  oxygen  in  wine  are:  cold 
tartaric  stabilization  (2.29  mg/L,  white  wines),  bottling  (1.38  mg/L  red  wines  and  1.11  mg/L  white  wines)  and  bag-in-box  filling 
(2.47  mg/L  red  wines;  2.22  mg/L  white  wines).  After  identifying  the  most  critical  technical  operations  in  wine  making,  preventive 
and  corrective  measures  had  to  be  taken  to  reduce  the  dissolved  oxygen  content  of  wines,  as  well  as  preventing  their  depreciation. 

Key  words:  Wine,  oxygen,  oxidation,  dissolution. 


1.  Introduction 

Oxygen  exists  as  a  gaseous  element  representing 
20.9%  of  earth  atmosphere,  under  standard 
temperature  and  pressure,  and  it  is  indispensable  for 
the  majority  of  the  living  organisms. 

The  oxygen  plays  a  crucial  role  in  enology.  Louis  [1] 
had  already  noted  that  oxygen  could  be  a  great  enemy 
of  wine  or  wine  making. 

Oxygen  can  dissolve  in  wine  at  several  stages 
during  winemaking,  bottling  and  aging.  It  plays  a 
beneficial  role  during  the  winemaking  process:  in 
promoting  yeast  multiplication,  in  the  maturation  and 
stabilization  of  wine  colour,  in  the  removal  of  the 
reduced  compounds  taste,  and  attenuation  of  the 
vegetable  character,  etc..  However,  when  the  wine  is 
exposed  to  excessive  amounts  of  oxygen,  it  can  have 
negative  effects  such  as  the  promotion  of  the  growth 
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of  unfavourable  microorganisms  in  wine  (acetic 
bacteria,  Brettanomyces),  the  degradation  of  colour 
and  the  production  of  organoleptic  defects,  etc.. 

The  oxidation  in  white  wine  causes  important 
sensory  alterations  in  the  colour  and  the  aroma.  The 
browning  is  characterized  by  progressive  darkening 
which  can  be  globally  characterized  by  absorbance  at 
420  mn  [2].  The  oxygen  induces  a  loss  of  freshness 
and  fruitiness,  and  develops  an  unpleasant  oxidized 
character,  which  is  generally  described  as 
non-desirable  flavours  of  "honey-like",  "boiled-potato", 
"cooked  vegetable",  "farm-feed"  and  "woody-like" 
[3-7]. 

Phenolics  compounds  are  the  main  oxidation 
substracts.  Various  researchers  have  studied  the 
impact  of  oxygen  exposure  on  the  polyphenolic 
composition  of  red  wines  [8-10].  Some  controlled  and 
continuous  02  is  introduced  into  red  wine 
(micro-oxygenation)  to  improve  wine  quality,  such  as 
colour  stability  and  intensity,  softening  of  astringent 
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tannins,  and  decreases  reductive  and  vegetative 
aromas. 

During  winemaking,  the  different  unit  operations 
may  integrate  some  oxygen  from  the  atmosphere  into 
the  wine.  The  operations  such  as  pumping,  transport, 
racking,  centrifugation,  stabilization,  filtration  and 
bottling  have  been  studied  by  various  authors  [11-17]. 

Various  researchers  have  monitored  the  dissolved 
oxygen  of  wine  at  bottling  and  during  storage  and  the 
influence  of  packaging  in  wine  conservation  [3,  6,  10, 
13,  14,  17-19]. 

The  Esporao  SA  is  one  of  the  biggest  wineries  in 
Portugal.  The  purpose  of  this  study  was  to  quantify 
the  dissolved  oxygen  in  red  and  white  wine,  during 
different  oenological  treatments  such  as  racking, 
tartaric  stabilization,  filtering  and  bottling,  to  identify 
the  critical  points  (the  unit  operations)  that  contribute 
most  to  the  integration  of  oxygen  in  wine  and  to 
establish  preventive  and  corrective  measures. 

2.  Experiments 

2.1  Wines 

The  white  or  red  wines  analyzed  were  produced 
during  2010  vintage  in  Esporao  S.A.  (Alentejo)  from 
different  varieties  (Arinto,  Antao  Vaz,  Perrum; 
Moreto,  Trincadeira  e  Castelao,  Touriga  Nacional, 
Cabernet  Sauvignon,  among  other).  All  wines  were 
obtained  following  the  usual  winemaking  processes. 

2.2  Wine  Treatments 

Red  and  white  wines  were  submitted  to  specific 
oenological  treatments  (Table  1)  with  systems, 
equipment  and  procedures  commonly  used  in  winery. 
The  increase  of  dissolved  oxygen  (OdiSS)  was 
calculated  as  the  difference  between  the  oxygen 
concentration  measured  before  (Obef)  and  immediately 
after  (Oaft)  each  treatment  (Odiss  =  Oafl  -  Obef). 

2.3  Measurement  of  Dissolved  Oxygen 

The  measurement  of  dissolved  oxygen  was  carried 
out  using  a  portable  oximeter  Crison  OXI  45  P. 


2.4  Statistical  Procedures 

(1)  Data  of  the  increase  of  dissolved  oxygen  are 
expressed  as  the  average  and  standard  deviation  of 
measurements  performed  for  the  different  wines. 

(2)  The  influence  of  oenological  treatments  was 
analyzed  by  non-parametric  procedures  by  the 
paired-samples  T-test  (before  and  after  treatment), 
with  the  Wilcoxon  signed  ranks  test,  using  the  SPSS 
statistics  software  (Version  17.0). 

3.  Results  and  Discussion 

The  level  of  dissolved  oxygen  at  a  given  time  in  the 
wine  depends  on  the  dissolution  and  consumption 
kinetics.  When  applied  to  a  wine  movement  (pumping, 
filling  and  stirring)  the  kinetics  of  oxygen  dissolution 
is  higher  than  the  consumption  by  the  wine's 
compounds.  Dissolved  oxygen  reacts  with  the 
constituents  of  the  wine  and  disappears  gradually, 
reaching  very  low  levels. 

After  racking,  the  evolution  of  dissolved  oxygen 
consumption  in  the  wine  was  monitored  over  time 
(Fig.  1).  It  was  found  that  the  kinetics  of  consumption 
varied  depending  on  the  type  of  wine,  having  been 
noted  that  the  white  wine  kinetics  were  slower  than 
the  red  wines  kinetics.  Red  wines  consume  oxygen 
faster  than  white  wines  because  of  its  constitution  of 
phenolic  compounds.  In  the  case  of  red  wines,  in  just 
four  days  phenolic  compounds  consumed  more  than 
half  of  dissolved  oxygen.  For  white  wines,  it  took 
more  than  30  days.  Increasing  temperature  increases 
the  kinetics  of  dissolved  oxygen  consumption. 

The  monitoring  of  oxygen  concentrations  before 
and  after  specific  oenological  treatments  provided  data 
on  the  levels  of  enrichment  of  the  winery  and  its 
variability. 

Statistical  analysis  performed  using  nonparametric 
paired  samples  revealed  that  the  enrichment  in 
dissolved  oxygen  is  significant  for  all  operations 
performed  (Table  2). 

The  increase  of  dissolved  oxygen  for  each  treatment 
is  shown  in  Table  3. 
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Table  1    Description  of  oenological  treatments. 


White  wine 

n 

Red  wine 

n 

Protein  stabilization  (bentonite  addition) 

4 

Racking  (bottom-bottom) 

110 

Continuous  tartaric  stabilization  (cooling  method) 

12 

Tartaric  stabilization  (metatartaric  acid  addition) 

12 

Filtration  (membrane  cartridges,  nominal  porosity  1-0.45  urn) 

29 

Tangential  filtration 

47 

Bottling  (glass  bottles) 

28 

Bottling  (glass  bottles) 

147 

Packaging  (bag-in-box,  5  L) 

9 

Packaging  (bag-in-box,  5  L) 

26 

n  =  number  of  analysed  wines. 


16  - 
14  - 

Red  wine 

13.6 

12  - 
10  - 
8  - 

6  - 

4  - 

2  - 

1.16 

0.72  D4g 

0.39 

0  - 

i 

0                   2  3 
Time  (days) 

Dissolved  oxygen  (mg/L)   ----Temperature  (5C) 

4 

White  wine 


15  -i 
14  - 
12  - 
10  - 

8  - 

6  - 

4  - 

2  - 

0  -- 


14.3 


11.2 


11.7 


12.3 


9.1 


3.07 


2.42 


2.27 


1.89 


1.29 


16  28  50 

Time  (days) 

•  Dissolved  oxygen  (mg/L)  ----Temperature  (K) 


64 


(a) 

Fig.  1    Evolution  of  dissolved  02  in  wine  after  racking:  (a)  red  wine,  (b)  white  wine. 


(b) 


Table  2  The  average  content  of  dissolved  oxygen  (mg/L),  before  and  after  each  oenological  treatment  and  p-value  of 
Wilcoxon  signed  ranks  test. 


Oenological  treatment 

Before  treatment 

After  treatment 

p-value  (Wilcoxon  test) 

Racking 

0.19  ±0.12 

0.65  ±0.34 

0.000 

Tartaric  stabilization 

0.28  ±0.25 

0.38  ±0.26 

0.002 

Red  wine 

Tangential  filtration 

0.14  ±0.49 

0.40  ±0.33 

0.000 

Bottling  (glass  bottles) 

0.49  ±0.42 

1.87  ±0.93 

0.000 

Packaging  (bag-in-box) 

0.98  ±  0.28 

3.45  ±0.86 

0.000 

Protein  stabilization 

1.22  ±0.85 

1.35  ±0.89 

0.042 

Cold  tartaric  stabilization 

0.60  ±0.37 

2.89  ±  1.80 

0.008 

White  wine 

Filtration  (membrane  cartridges) 

0.20  ±0.17 

0.71  ±0.16 

0.000 

Bottling  (glass  bottles) 

0.67  ±0.24 

1.78  ±0.42 

0.000 

Packaging  (bag-in-box) 

0.65  ±0.22 

2.87  ±0.96 

0.008 

Racking  increased  the  dissolved  oxygen  content  by 
on  average  of  0.46  mg/L,  exceeding  that  observed  [11,  15]. 
The  oxygen  enrichment  in  the  racking  occurs 
essentially  at  the  beginning  and  end  of  the  transfer, 
due  to  the  dissolution  of  air  contained  in  the  pipes  and 
tanks.  The  racking  should  be  performed  bottom  to 
bottom  of  the  tanks,  and  it  uses  inert  gas  in  the  pipes 


and  tanks  before  starting  the  transfer  process  [13,  14]. 

For  the  tartaric  stabilization  of  red  wine,  metatartaric 
acid  was  added  just  prior  to  bottling.  The 
homogenization  of  the  wine  was  achieved  with  nitrogen. 
The  addition  of  the  metatartic  acid  leads  to  an  increase  of 
0.10  mg/L  of  dissolved  oxygen,  which  is  a  reduced 
amount  when  compared  to  the  cold  tartaric  stabilization. 
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Table  3    Increase  of  dissolved  oxygen  (mg/L)  in  each  oenological  treatment. 


Oenological  treatment 
Racking 

Tartaric  stabilization 

Red  wine  Tangential  filtration 

Bottling  (glass  bottles) 
Packaging  (bag-in-box) 
Protein  stabilization 
Cold  tartaric  stabilization 

White  wine  Filtration  (membrane  cartridges) 

Bottling  (glass  bottles) 
Packaging  (bag-in-box) 

The  important  enrichment  that  occurred  during  cold 
tartaric  stabilization  of  the  white  wines  was  on 
average  2.21  mg/L  (Table  3),  and  these  results  are 
consistent  with  data  found  by  other  authors  [12,  13, 
15].  The  sharp  temperature  drop  caused  an  increase  in 
the  oxygen  dissolution  rate,  therefore,  it  is  important 
to  protect  the  wine  from  the  oxygen  in  the  operation. 

The  protein  stabilization  of  white  wine  by  the 
addition  of  bentonite  showed  a  small  enrichment  in 
dissolved  oxygen.  It  was  noted  that  this  operation  did 
not  contribute  a  lot  to  the  increase  of  dissolved  oxygen 
in  the  wine,  similar  to  what  was  observed  for  the 
tangential  filtration  of  red  wines,  in  which  the 
dissolved  oxygen  was  lower  than  that  reported  by 
Vidaletal.  [11,  13]. 

The  filtration  by  membrane  cartridges  increased  the 
dissolved  oxygen  content  mainly  due  to  air  inlets  in 
the  circuit,  in  particular  through  the  pipe  connections, 
the  pump  and  filter  body.  The  dissolved  oxygen  in 
wine  is  higher  than  that  reported  in  Refs.  [11,  13,  15], 
these  authors  noted  average  values  of  0.10  mg/L  and 
0.06  mg/L,  respectively.  This  is  due  to  the  length  and 
the  diameter  of  the  tubing  as  well  as  the  turbulent  flow, 
which  influences  the  kinetics  of  oxygen  dissolution. 

Unlike  the  results  of  Vidal  et  al.  [13,  14],  filling  in 
BIB  (bag-in-box)  has  contributed  to  a  greater 
enrichment  in  dissolved  oxygen  than  bottling.  This  is 
due  to  a  problem  with  an  inertized  filling  machine 
used  to  fill  BIB  containers.  To  reduce  the  risk  of  wine 


02D  (mg/L)  (mean  ±  standard  deviation) 

0.46  ±0.31  (n=  110) 

0.10  ±0.06  (n=  12) 

0.26  ±0.30  (m  =  47) 

1.38  ±0.49  (m=147) 
 2.47  ±  0.93  {n  =  26)  

0.13  ±0.05  («  =  4) 

2.29  ±  1.78  (n=  12) 

0.51  ±0.18  (n  =  29) 

1.11  ±0.47  (n  =  28) 
 2.22  ±  0.77  (n  =  9)  

oxidation  in  BIB,  the  concentration  of  SO2  was 
increased  up  to  a  maximum  of  50  mg/L,  (1  mg/L 
oxygen  oxidizes  4  mg/L  SO2),  so  that  the  SO2  protects 
the  wine  from  oxidation  and  browning  [19]. 

Bottles  were  filled  using  a  filling  machine  inertized 
with  nitrogen.  However,  this  operation  was 
responsible  for  some  increase  in  dissolved  oxygen 
content,  both  in  white  wine  and  in  red  wine.  These 
results  are  in  agreement  with  those  previously 
described  [13,  14,  18].  The  filling  is  the  main  critical 
point  for  wine  enrichment  in  oxygen. 

4.  Conclusions 

The  rationalization  of  the  operating  conditions  and 
the  reasonable  use  of  neutral  gases,  especially  at  the 
beginning  and  the  end  of  operation,  are  necessary  to 
control  and  decrease  significantly  the  oxygen  pick  up 
in  all  the  operations. 

The  most  important  enrichments  occur  especially 
during  the  operations  of  cold  tartaric  stabilization  and 
during  the  bottling  and  packaging,  especially  if  no 
protection  is  taken  to  protect  the  wine  from  the 
oxygen  in  the  air. 

Controlling  the  oxygen  ingress  during  filling  (bottle 
or  BIB)  is  an  important  way  of  improving  the  quality 
and  the  shelf-life  of  wine. 
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Abstract:  The  Curaua  (Ananaserectifolius  LB  Smith)  is  a  species  belonging  to  the  family  Bromeliaceae.  Occurs  in  the  states  of  Para, 
Amazonas,  Amapa,  Acre,  Mato  Grosso  and  Goias,  has  two  varieties,  purple  and  white,  and  its  fibers  are  used  in  automotive  and 
textile  industries  due  to  its  strength,  softness  and  lightweight.  Currently,  only  the  fiber  is  used  in  industry,  the  rest  is  considered  waste. 
However,  this  residue  contains  compounds  with  important  properties  to  be  discovered  and  studied.  An  enzyme  complex  found  in  this 
residue  is  bromelain,  a  group  of  proteolysis'  enzymes  with  application  in  several  areas  such  as  in  the  food,  pharmaceutical  and 
cosmetics  industries.  Today,  its  use  is  aiming  the  pharmaceutical  industry,  the  production  of  ointments,  gels,  creams  and  lotions 
because  they  offer  a  wide  range  of  therapeutic  efficacies:  antiedemas,  anti-inflammatory,  antithrombotic  and  fibrinolytic  activities. 
This  work  studied  the  enzyme  purification  and  recovery  presented  in  the  leaves  of  Curaua  by  means  of  an  ATPS  (aqueous  two-phase 
system)  PEG  4000/potassium  phosphate.  The  protein  content  was  measured  by  Bradford  reagent  and  the  enzymatic  activity  was 
measured  by  the  Biuret  reagent.  Batch  assays  were  performed  aiming  the  enzyme  extraction  and  recuperation,  using  the  partition 
coefficient  as  indicator.  It  was  used  pH  7.0,  8.0  and  9.0,  varying  the  proportional  composition  between  the  polymeric  and  the  saline 
phases  (tie-lines).  The  white  variety  of  Curaua  was  used.  The  best  purification  factor  was  4.53  and  about  113.04  |rg/mL  of  total 
protein  measured  by  Bradford. 

Key  words:  Curaua,  purification,  bromelain. 


1.  Introduction 

Enzymes  are  organic  molecules  present  in  cells  of 
living  organisms,  with  the  specific  function  of 
catalyzing  chemical  reactions.  The  enzyme  starts 
increasing  velocity  of  a  chemical  reaction  by  acting  as 
a  catalyzer  [1].  A  catalyzer  does  not  modify  the 
balance  of  a  reaction,  however,  it  reduces  the  amount 
of  energy  required  to  make  the  reaction  to  happen. 

Economically  viable,  the  use  of  enzymes  in  various 
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industrial  sectors  has  been  increasing  over  the  years  [2-3]. 
Among  enzymes  with  industrial  use,  the  bromelain 
stands  out  due  to  its  diversity  in  applications:  leather 
processing  in  textile  industries;  softening  fibers  and 
detergent  production  [4];  in  food  industry  providing 
meat  softening  [5];  and  in  breweries  it  is  used  to 
purify  beers  hydrolyzing  tannin-protein  complexes 
produced  during  fermentation  [5]. 

Bromelain  is  also  used  in  the  pharmacological 
research  area  due  to  its  anti-inflammatory  effect  [6] 
and  its  potential  to  help  treating  health  problems  such 
as  angina,  indigestion  and  breathing  problems  [7].  In 
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therapies  against  cancer,  it  is  used  to  increase  lysis  of 
cancer  cells  [8]. 

Curaua  (Ananaserectifolius)  is  a  little  known  and 
studied  species.  Its  cultivation  started  at  Lago  Grande 
de  Curuai,  city  of  Santarem  in  Para  state,  and  has  been 
expanding  to  other  regions  of  the  same  state. 
Researches  have  shown  that  Curaua  fiber  presents 
excellent  quality,  being  compared  to  fiberglass  due  to 
its  resistance,  softness  and  low  weight  [9,  10].  In 
automotive  industry,  it  has  been  used  to  produce 
sunshade  for  trucks  [11]. 

Currently,  the  separation  and  purification  process  of 
bioproducts  is  a  very  important  segment  for  industries, 
as  it  can  represent  from  80%  to  90%  of  the  production 
cost.  Therefore,  the  development  of  an  efficient  and 
low  cost  process  is  highly  relevant  [12]. 

The  liquid-liquid  extraction  has  been  arousing 
interest  as  an  intermediate  separation  step  that  either 
replaces  more  expensive  methods  of  separation  or 
reduces  the  number  of  separation  steps  required  by  the 
process.  This  separation  process  is  based  on  the 
distribution  of  the  solute  between  the  phases  and 
partial  miscibility  of  the  liquids  [13]. 

The  liquid-liquid  extraction  process  consists  of 
transferring  a  substance  from  a  liquid  mixture  to 
another  immiscible  (or  partially  miscible)  liquid  phase 
by  putting  them  in  contact.  This  process  is  widely 
used  in  chemical  and  pharmaceutical  industries,  such 
as  in  the  recovery  of  antibiotics  or  organic  acids  from 
fermentation  broths. 

Nevertheless,  the  application  in  the  purification  of 
proteins  is  still  limited  mainly  due  to  the  possibility  of 
protein  denaturation  if  in  contact  with  organic  solvents, 
yielding  a  useless  product  [14-17].  The  ATPS  (aqueous 
two-phase  systems)  are  created  by  the  addition  of 
aqueous  solutions  of  two  hydrophilic  polymers  such 
as  PEG  (polyethylene  glycol)  and  dextran  or  a  polymer 
and  a  salt,  such  as  PEG  and  potassium  phosphate. 

These  substances  provide  the  appropriate  and 
convenient  method  for  extracting  substances  of 
biological  origin,  as  the  formation  of  phases  between 


70%  and  90%  of  water  provides  a  bland  environment 
to  work  with  biologically  active  compounds, 
preserving  their  molecular  stability  and  thus,  allowing 
their  processing  on  this  medium  [18]. 

The  pH  in  the  medium,  besides  influencing  in  the 
maintenance  of  enzymatic  activity,  also  affects  the 
charge  distribution  on  the  surface  of  the  protein.  Thus, 
adjusting  the  pH  of  the  medium  to  values  close  to  the  IP 
(isoelectric  point)  of  the  protein  (pH  value  where  the 
molecule  presents  electric  charge  near  zero),  the  effects 
that  act  over  the  partition  in  ATPS  are  reduced  [19]. 

Ferreira  [20]  analyzed  the  purification  of  the 
bromelain  enzyme  extracted  from  the  stem  of  the 
pineapple  (Ananascomosus)  in  ATPS  (aqueous 
two-phase  system)  PEG/potassium  phosphate  and  got 
better  results  using  polyethylene  glycol  4000,  rather 
than  polymer  with  lower  molecular  weight  (PEG 
1500). 

Thus,  this  work  had  as  objective  obtaining  PF 
(purification  factors)  of  the  bromelain  enzyme  derived 
from  white  variety  of  Curaua  (Ananaserectifolius  LB 
Smith)  through  the  ATPS  PEG  4000/potassium 
phosphate. 

2.  Materials  and  Methods 

PEG  with  molecular  weight  of  4000  was  obtained 
from  Sigma  (Switzerland).  Milli-Q-quality  distilled 
water  was  used.  Potassium  hydrogen  phosphate, 
comassie  brilliant  blue  G,  di-sodium  hydrogen 
phosphate,  sodium  hydroxide  were  obtained  from 
Synth  (Sao  Paulo,  Brazil). 

The  following  equipments  were  used:  Mars 
electronic  scale  AL  200;  an  Arno  multiprocessor;  pH 
meter  Analyser  pH  300;  thermostatic  bath  FANEM; 
automatic  micropipettes  and  spectro-photometer 
UV/VIS  Spectronic  2 ID. 

The  concentration  values  of  PEG  and  phosphate  for 
each  tie  line  were  reproduced  according  to  Ref.  [20]. 

Curaua's  leaves  were  collected  from  field  at  Sao 
Manoel  Experimental  Farm  in  Botucatu,  Sao  Paulo. 
Upon  arriving  on  lab,  they  were  washed  with  distilled 
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water,  dried  with  paper  towel  and  stored  in  plastic 
bags,  refrigerated  until  use. 

2. 1  Samples  Preparation 

Buffer  and  Curaua' s  leaves  were  weighted  in  1:1 
proportion;  then,  they  were  grinded  using  a  high 
efficiency  multiprocessor  and  the  extract  was  filtered 
using  nylon  filter  removing  fibers  and  other  present 
particles. 

2.2  Preparation  of  Reactive  Biuret 

It  was  dissolved  in  500  mL  of  distilled  water,  1.5  g 
of  copper  pentahydrate  and  6.0  g  of  tartarate  of  sodium 
and  potassium  sulphate,  and  it  was  added  300  mL  sodium 
hydroxide  solution  10%  under  constant  agitation,  and 
finally  added  distilled  water  for  a  1  L  solution. 

2.3  Determination  of  Bromelain  Activity 

In  order  to  determine  the  activity  of  bromelain,  we 
used  the  Biuret  reagent  and  a  solution  to  5  g/L  BSA 
(bovine  serum  albumin),  this  method  is  based  on  the 
reaction  of  the  reactive  Biuret,  consisting  of  a  mixture 
of  copper  and  sodium  hydroxide  with  a  complexing 
that  stabilizes  copper  in  solution.  Readings  were  taken 
in  a  spectrophotometer  at  540  nm,  and  finally 
calculated  the  enzymatic  activity,  according  to  Eq.  (1). 


AE  = 


Vreactor(l)xl06(M^) 


(1) 


MMBSAxVenzyme  (ml)  x  timereaction  (min) 

where,  Vreactor  is  the  reactor  volume;  MMBSA  is  the 
molar  mass  of  PEG  and  Venzyme  is  the  enzyme  volume 
used  for  the  analysis. 

Total  protein  in  samples  of  bottom  and  top  phases 
were  determined  by  the  Bradford  method  [21]. 

Table  1  Partition  of  the  bromelain  enzyme  derived  from 
phosphate  in  pH  7.0. 


2.4  Preparation  of  Buffers 

Buffer  solutions  were  prepared  in  pH  7.0,  8.0  and 
9.0,  from  a  solution  of  mono  potassium  phosphate 
dibasic  to  15%  (w/w). 

For  this,  it  was  prepared  standard  solutions  of 
monobasic  potassium  phosphate  (solution  A)  and 
bibasic  potassium  phosphate  (solution  B)  mixing  them 
into  a  beaker  containing  silver  chloride  electrode  for 
measuring  the  pH,  until  the  desired  pH. 

2.5  Determination  of  the  Purification  Factor 

To  calculate  the  PF  it  was  used  the  Eq.  (2). 

PF  =  — —  (2) 

crude 

where,  EAcrude  is  crude  enzymatic  activity  (U)  and  EA 
is  enzymatic  activity  of  top  and  bottom  phase. 

3.  Results  and  Discussions 

Tables  1  -3  demonstrate  the  influence  of  the  pH  and 
the  concentration  of  the  PEG  4000/potassium 
phosphate  system  regarding  the  purification  of  the 
bromelainenzyme  derived  from  white  variety  of 
Curaua. 

The  ATPS  that  presented  better  result  in  the 
purification  of  the  bromelain  enzyme  derived  from 
white  Curaua  was  the  PEG  4000/potassium  phosphate 
system  in  pH  7.0  ("tie  line"  2  and  3);  the  PEG 
4000/potassium  phosphate  system  in  pH  9.0  ("tie  line" 
1,  2  and  3)  also  presented  considerable  values. 

These  results  are  based  on  the  values  of  enzymatic 
activity  and  purification  factor. 

Analyzing  the  behavior  of  the  three  "tie  lines"  of 
the  ATPS  PEG  4000/potassium  phosphate  in  pH  7.0, 

leaves  of  white  Curaua  using  the  ATPS  PEG  4000/potassium 


Tie  line 

Proteins  (jig/mL) 

EAQJ) 

PF 

1 

Top  phase 
Bottom  phase 

18.73 
32.57 

0.004 
0.004 

2.81 
2.82 

2 

Top  phase 
Bottom  phase 

22.02 
35.75 

0.006 
0.006 

3.70 
3.80 

3 

Top  phase 
Bottom  phase 

30.62 
54.72 

0.007 
0.007 

4.53 
4.29 
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Table  2  Partition  of  the  bromelain  enzyme  derived  from  leaves  of  white  Curaua  using  the  ATPS  PEG  4000/potassium 
phosphate  in  pH  8.0. 


Tie  line 

Proteins  ()ig/mL) 

EA  (U) 

PF 

1 

Top  phase 

21.64 

0.004 

1.36 

Bottom  phase 

47.88 

0.004 

1.30 

2 

Top  phase 

22.23 

0.008 

2.23 

Bottom  phase 

37.74 

0.006 

1.84 

3 

Top  phase 

32.32 

0.005 

1.41 

Bottom  phase 

87.74 

0.004 

1.28 

Table  3  Partition  of  the  bromelain  enzyme  derived  from  leaves  of  white  Curaua  using  the  ATPS  PEG  4000/potassium 
phosphate  in  pH  9.0. 


Tie  line 

Proteins  (|ig/mL) 

EA(U) 

PF 

1 

Top  phase 

35.39 

0.008 

2.39 

Bottom  phase 

113.04 

0.007 

2.23 

2 

Top  phase 

28.97 

0.008 

2.47 

Bottom  phase 

74.14 

0.008 

2.42 

3 

Top  phase 

30.65 

0.007 

2.22 

Bottom  phase 

68.16 

0.008 

2.31 

it  is  possible  to  notice  the  increasing  trend  of  the 
purification  factors  as  the  length  of  the  mooring  lines 
increases;  however,  the  ATPS  PEG  4000/potassium 
phosphate  in  pH  8.0presented  the  inverse  behavior; 
higher  values  for  the  purification  factors  were 
obtained  in  lower  concentrations  of  the  components  of 
the  system. 

Sbruzzi  [22]  showed  in  her  study  about  the  costs  of 
this  process  (ATPS),  using  the  fiber  of  curaua,  purple 
and  white  variations,  and  the  estimated  sale  price, 
showing  its  economic  viability,  whether  it  is  for  use  as 
a  pre-process  in  traditional  chromatographic  systems 
or  direct  commercialization,  as  an  alternative  for  the 
bromelain  of  pineapple. 

4.  Conclusions 

The  values  of  the  purification  factor  of  the 
bromelain  enzyme  derived  from  white  Curaua  were 
affected  by  concentrations  of  the  polymeric  and  saline 
phases,  as  well  as  by  the  pH  value  in  the  medium. 

The  best  result  was  obtained  working  with  the 
lowest  pH  studied  in  this  work  but  with  highest 
polymeric  concentration  ("tie  line"  2  and  3)  of  the 
ATPS  PEG  4000/potassium  phosphate. 

The  technique  of  the  aqueous  two-phase  PEG 


4000/potassium  phosphate  systems  proved  being 
promising  in  the  purification  of  the  bromelain  enzyme 
derived  from  the  white  variety  of  the  Curaua. 

It  can  be  used  to  extract  this  proteolytic  enzyme 
during  the  extraction  of  the  white  Curaua  from 
automotive  industry,  for  instance,  maximizing  profit, 
reducing  environmental  impact  and  the  waste. 
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Abstract:  Sulfate  reducing  bacteria  were  isolated  from  the  soil  sample  of  Roaster-Acid  unit  of  Debari  Zinc  Smelter  of  Hindustan 
Zinc  Ltd.,  India  and  were  adapted  for  effluent  water  of  Zinc  hydrometallurgy  plant  to  sustain  zinc  concentration  of  about  500  ppm 
and  iron  concentration  of  about  160  ppm  with  high  acidity  (pH  <  1.0).  Nutrient  broth,  whey,  ethanol  and  sucrose  were  tested  as 
carbon  source  for  SRB  (sulfate  reducing  bacteria)  out  of  which  whey  was  found  as  the  best  electron  donor  to  expedite  the  rate  of 
biological  treatment  with  SRB.  A  two  stage  anaerobic  bioreactor  was  developed  with  a  treatment  capacity  of  30  L  effluent  water  in  4 
days.  The  concentration  of  heavy  metals  in  the  treated  water  was  within  the  permissible  limits.  The  suggested  methodology  generates 
no  solid  waste,  is  environment  friendly  and  may  be  commercially  implemented  for  treatment  of  zinc  plant  effluent. 

Key  words:  Sulfate  reducing  bacteria,  anaerobic  treatment,  isolation,  hydrometallurgy  effluent. 


1.  Introduction 

Industrial  wastewaters  especially  originating  from 
acid  mine  drainage  and  hydrometallurgical  industries 
contain  high  concentration  of  heavy  metals  as  sulfates 
with  free  acidity.  The  improper  treatment  and/or 
discharge  of  this  acidic  water  can  have  serious 
environmental  impacts  and  can  unbalance  sulfur 
cycle  [1].  Conventionally,  the  effluent  water  is  being 
neutralized  with  lime  followed  by  alum  flocculation. 
The  method  is  expensive  and  results  in  high  total 
dissolved  solids  in  treated  water  along  with  generation 
of  residual  sludge,  which  is  dumped  in  SLF  (secured 
land  fill)  sites  and  needs  a  large  capital  expenditure 
resulting  in  permanent  loss  of  fertility  of  costly  land. 
Membrane  separation  is  also  being  used  as  an 
alternative  route  in  effluent  treatment  to  reduce  solid 
waste  load.  However,  relatively  high  cost  of  frequent 
membrane  replacement  and  energy  consumption 
restricts  commercial  use  of  this  technique  [2].  Recently, 
biotechnology  has  shown  new  opportunities  and 
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challenges  in  many  sectors  especially  in  environmental 
remediation  of  industrial  wastes  (solid/liquid)  to 
achieve  better  environment  management  and  pollution 
control. 

Present  research  is  conducted  with  an  aim  of  in-situ 
precipitation  of  heavy  metals  as  metal  sulfides  from 
effluent  of  zinc  hydrometallurgical  plant.  The 
advantage  of  sulfide  precipitation  include:  better  water 
quality  due  to  lower  solubility  of  metal  sulfides 
compared  with  their  corresponding  hydroxides; 
selective  precipitation  of  high  grade  metal  sulfide 
concentrates  that  can  be  reprocessed  to  earn  revenue; 
reduction  in  sludge  volume;  and  toxicity  and 
opportunities  in  reagent  saving.  A  better  overall 
environment  solution  can  be  attained  due  to  metal 
recycle  and  effective  solid-liquid  separation.  Under 
anaerobic  conditions,  SRB  (sulfate-reducing  bacteria) 
oxidize  simple  organic  compounds  completely  to  CO2 
by  utilizing  sulfates  as  electron  acceptor  and  generate 
sulfides  and  alkalinity  [3].  SRB  utilize  a  number  of 
volatile  fatty  acids  and  hydrogen  as  electron  donors  to 
expedite  its  activity  [4].  These  sulfides  subsequently 
precipitate    metals    as    metal    sulfides    in  in-situ 
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conditions  with  saving  of  costly  reagents  like  lime  and 
sodium  sulfide.  Metal  sulfides  can  also  be  precipitated 
using  biogenic  H2S  in  ex-situ  conditions  [5,  6].  Hybrid 
H2S  bioreactor  concept  has  also  been  introduced  by 
controlling  methanogenic  reactions  of  typical 
anaerobic  sludge  reactor  [7].  In  case  of  anaerobic 
bacteria,  sulfates  are  the  main  source  of  H2S 
production  and  its  concentration  results  in  major 
impact  on  SRB  growth  [8].  SRB  consumes 
approximately  one  to  two  moles  of  protons  per  mole 
of  sulfate  that  is  reduced  and  produces  approximately 
two  equivalents  of  alkalinity  per  mole  of  reduced 
sulfate  [9].  The  electron  accounting  can  be  made  by 
assuming  complete  oxidation  of  organic  compounds  [5]. 
SRB  are  notable  for  their  end  product,  hydrogen 
sulfide,  which  is  far  more  noticeable  substance  than 
sulfates  due  to  its  chemical  properties  and 
physiological  effects  [10].  SRB  could  be  carbon 
limited  in  organic  medium  containing  cow  manure 
and  saw  dust  since  these  materials  are  largely 
polymers  of  cellulose  and  hemicellulose,  which  are 
insoluble  organic  matters  [11,  12]. 

This  paper  includes  isolation,  screening  and 
adaptation  of  SRB  from  hydrometallurgical  smelter  of 
Hindustan  Zinc  Limited  at  Udaipur  in  India. 
Evaluation  of  culture  performance  on  bioreactor  of  30  L 
effluent  treatment  capacity  has  also  been  analyzed  and 
discussed. 


2.  Experimental  Materials  and  Methods 

2.1  Effluent  Water 

Effluent  water  from  DZS  (Debari  Zinc  Smelter) 
was  collected  in  sealed  poly  vinyl  containers  and  the 
heavy  metals  were  analyzed  by  Atomic  Absorption 
Spectrophotometer  (Varian  220),  while  acidity  and 
sulfates  were  analyzed  by  conventional  titrimetry  and 
gravimetric  methods.  The  analysis  of  untreated 
(collected  from  DZS)  and  treated  effluent  (anaerobic 
bioreactor)  is  shown  in  Table  1 . 

2.2  Isolation  and  Selection  of  SRB 

To  isolate  SRB,  soil  samples  were  collected  from 
different  locations  of  DZS  of  Hindustan  Zinc  Ltd.  The 
locations  were:  roaster,  acid  plant,  effluent  plant, 
oxidation  pond,  cow  dung  sample  and  concentrate 
storage  yard.  The  samples  were  serially  diluted  for 
10"10  dilutions  with  proper  mixing  and  poured  in 
Petridishes  filled  with  Postgate  media  [13]  using  agar 
as  solidifying  medium.  The  sterile  conditions  were 
maintained  throughout  the  activities.  Petridishes  were 
incubated  under  anaerobic  conditions  at  30  °C  for  48  h. 
The  anaerobic  conditions  in  the  chamber  were 
maintained  by  creating  vacuum  and  filling  it  with  N2 
gas.  The  sterility  inside  incubation  chamber  was 
maintained  by  cleaning  the  chamber  with  ethanol  and 
further  N2  gas  was  passed  through  air  filter. 


Table  1 

Chemical  analysis  of  untreated  and  treated  effluent. 

Elements 

Untreated  effluent  (collected  from  DZS) 

Treated  effluent  (using  anaerobic  bioreactor) 

Concentration 

Concentration 

Zn 

440  ppm 

4.8  ppm 

Cu 

1  ppm 

<  0.5  ppm 

Fe 

160  ppm 

1.8  ppm 

As 

1  ppm 

<  0.5  ppm 

Sb 

<  1  ppm 

<  0.5  ppm 

Mn 

20  ppm 

<  0.5  ppm 

Cd 

30  ppm 

<  0.5  ppm 

Ni 

1  ppm 

<  0.5  ppm 

Pb 

5  ppm 

<  0.5  ppm 

SO42" 

8.76  gpl 

0.26  gpl 

Acidity 

4.40  gpl 

pH  =  6.8 
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Developed  bacterial  colonies  were  streak  plated  and 
incubated  under  anaerobic  conditions  for  48  h.  The 
streak  plating  was  repeated  five  times  to  screen  the 
anaerobic  SRB.  All  experiments  were  conducted  in 
triplicates  (45  sets  for  each  step)  to  isolate  the  purer 
SRB  culture.  The  screened  colonies  were  inoculated 
separately  in  sealed  flasks  containing  Postgate  media 
and  incubated  under  anaerobic  conditions.  The 
maximum  grown  cultures  were  taken  out  (minimum 
one  from  each  collection  site)  and  TSI  (triple  sugar 
iron)  test  was  conducted  to  confirm  the  activity  of 
isolates.  Selection  of  SRB  was  done  by  ascertaining 
its  efficiency  towards  sulfate  reduction  in  ten  times 
diluted  effluent  water  of  Debari  Zinc  Smelter.  400  mL 
of  this  diluted  effluent  was  taken  in  500  mL  conical 
flasks  and  inoculated  with  respective  strains  (10%  v/v). 
The  flasks  were  incubated  for  twenty  days  and  sulfates 
were  analyzed  on  alternate  days.  This  was  done  to 
supplement  the  conclusion  drawn  from  TSI  test  for 
selection  of  SRB.  The  batches  were  named  as  Ctrl,  St  1, 
St  2,  St  3,  St  4,  St  5  and  St  6  for  control  (inoculated 
with  methanol  in  place  of  bacterial  strain),  roaster, 
acid  plant,  ETP  (effluent  treatment  plant),  cow  dung, 
oxidation  pond  ETP  and  concentrate  storage  locations. 

2.3  Adaptation  of  SRB 

Selected  SRB  was  developed  in  Postgate  media  and 
was  used  as  inoculum  for  acclimatization  with  DZS 
effluent.  DZS  effluent  was  diluted  for  ten  times  in  a 
way  so  that  the  effluent  and  nutrient  broth  (Postgate 
media)  ratio  can  be  maintained  as  8:2  (v/v).  400  mL  of 
diluted  effluent  was  taken  in  to  500  mL  conical  flasks 
and  was  inoculated  (10%  v/v)  with  isolated  SRB 
culture  separately.  Sterile  N2  gas  was  purged  into  each 
conical  flask  to  expel  the  dissolved  oxygen  out.  The 
flasks  were  sealed  and  incubated  for  20  days.  This 
treated  effluent  was  further  used  as  inoculum  for 
diluted  DZS  effluent  and  similar  procedure  was 
followed  for  the  raw  effluent.  Flasks  were  harvested 
on  alternate  days  and  were  analyzed  for  zinc,  Iron  & 
pH.  The  treatment  of  effluent  without  dilution  was 


repeated  two  times  to  acclimatize  bacteria  and  after 
which  samples  were  analyzed  daily. 

2.4  Microorganism 

Best-acclimatized  culture  was  named  as 
SRB-CRDL  and  was  maintained  in  Postgate  media 
under  anaerobic  conditions.  The  strain  was 
subcultured  at  every  20  days. 

2.5  Selection  of  Carbon  Source  as  Electron  Donor 

It  was  essential  to  observe  the  effect  of  electron 
donor  as  well  as  selection  of  its  cheaper  source  since; 
SRB  could  mediate  sulfate  reduction  to  a  greater 
degree  when  simple  carbon  and  energy  sources  were 
available  to  catalyze  the  reactions  [9].  Nutrient  broth 
(2:8),  whey  (1:9),  ethanol  (1:9)  and  sucrose  (1:9) 
were  used  as  carbon  source  in  effluent  water. 
SRB-CRDL  was  inoculated  separately  in  500  mL 
sealed  conical  flasks  and  flasks  were  incubated  in  the 
similar  fashion  after  purging  sterile  N2  as  discussed 
earlier.  The  flasks  were  harvested  for  every  alternate 
day  and  were  analyzed  for  zinc,  iron,  manganese, 
lead,  cadmium  and  pH  value.  The  data  collected  was 
used  not  only  for  selection  of  electron  donor  but  also 
for  optimization  of  days  for  maximum  heavy  metals 
removal  and  to  later  design  the  size  of  the 
bioreactors. 

2.6  Semi  Batch  Type  Anaerobic  Bioreactors 

Two-stage  semi  batch  type  bioreactors  were 
installed  to  increase  the  treatment  capacity  of  the 
system.  In  first  stage,  two  bioreactors  of  10  L  capacity 
were  run  for  4  days  and  5  L  from  each  was  then 
transferred  to  another  anaerobic  bioreactor  of  equal 
capacity  for  second  stage  treatment,  while  first  stage 
bioreactors  were  again  leveled  up  with  fresh  effluent. 
The  liquid  streams  were  transferred  through  peristaltic 
pumps  under  anaerobic  and  sterile  conditions.  The 
experimental  set  up  was  continued  for  more  than  30 
days.  In  steady  state  conditions  total,  30  L  effluent 
was  treated  in  4  days. 
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3.  Results  and  Discussion 

3.1  Isolation  and  Selection  ofSRB 

The  results  of  TSI  test  for  isolated  bacterial  strains 
are  given  in  Fig.  1  and  Table  2  which  shows  that  the 
SRB  strain  isolated  from  Acid  Plant  of  Debari  Zinc 
Smelter  is  the  best  out  of  the  strains  isolated  from  other 
places.  Fig.  2  shows  sulfate  reduction  in  ten  times 


diluted  effluent  water  using  different  isolated  cultures. 
After  20  days  of  microbiological  treatment,  experiment 
confirms  that  the  strain  isolated  from  Acid  Plant  is  the 
best  out  of  other  strains.  This  supports  the  result 
concluded  by  TSI  test.  The  bacterial  strain  collected 
from  the  Roaster-Acid  soil  was  more  adaptable  for 
sulphate  ions  and  hence,  it  provided  the  best  results  in 
terms  of  sulphates  reduction  from  the  effluent. 


Fig.  1    TSI  test. 


Fig.  2    Reduction  in  sulfate  concentration  in  ten  times  diluted  DZS  effluent  w.r.t.  different  strains. 
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Table  2    Analysis  of  TSI  test. 


Strain  isolate  No. 

Ctrl 

St  1 

St  2 

St  3 

St  4 

St  5 

St  6 

Strain  collection  place 

Control 

Roaster 

Acid 
Plant 

ETP  (effluent  treatment 
plant) 

COW 

dung 

Oxidation  Concentrate  Remarks 
pond  ETP  storage 

Glucose  fermentation 

X 

vv 

V 

vv 

V 

V 

Red  color  turns 
yellow  in  butt 

Lactose/ sucrose 
fermentation 

X 

X 

V 

X 

V 

X 

Red  color  turns 
yellow  on  slant 

H2S  gas  production 

X 

V 

vv 

X 

VV 

V 

V 

Blackish  color  in 
butt  after  yellow 

Gas  production  (H2/CO2) 

X 

X 

X 

X 

X 

X 

Creaks  or  shifting 
of  slant 

3.2  Adaptation  ofSRB 

Fig.  3  represents  the  adaptation  graph  of  the 
selected  SRB.  The  adaptation  was  started  with  10 
times  dilution  for  20  days  while  second  peak  in  the 
graph  is  for  5  times  dilution  batch  and  was  run  for 
another  20  days,  which  is  followed  by  adaptation  of 
SRB  without  any  dilution  of  effluent  water.  The 
similar  experiment  is  repeated  for  two  more  times. 
After  92  days  of  adaptation  experiments,  the  bacteria 
were  adapted  to  reduce  heavy  metals  concentration  in 
acidic  effluent  water  to  permissible  limits  as  per 
environment  norms.  The  adapted  bacteria  are  capable 
of  reducing  heavy  metals  with  major  constituent  of 


zinc  500  ppm  concentration  to  less  than  5  ppm  and  pH 
of  less  than  1.0  to  almost  neural  pH  within  12  days  of 
biological  treatment.  Only  Zn  and  Fe  were  considered 
in  establishing  the  parameters  since  these  are  the 
major  constituents  of  the  effluent  water. 

3.3  Selection  of  Carbon  Source  as  Electron  Donor 

It  was  required  to  expedite  the  heavy  metal  removal 
to  handle  larger  effluent  volume.  SRB  could  mediate 
sulfate  reduction  to  a  greater  degree  when  simple 
carbon  and  energy  sources  were  available  to  catalyze 
the  reactions  [9].  Table  3  represents  reduction  in 
heavy  metals  concentration  (zinc,  iron,  cadmium, 
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Fig.  3    Adaptation  of  SRB  isolate. 
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Table  3    Effect  of  electron  donors  on  percentage  reduction  of  heavy  metals  and  sulfate  concentration. 


Zinc  Iron 


Days 

Nutrient  broth  Whey 
(2:8)  (1:9) 

Ethanol 
(1:9) 

Sucrose 
(1:9) 

Nutrient  broth  Whey 
(2:8)  (1:9) 

Ethanol  (1:9) 

Sucrose 
(1:9) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

43.06 

49.16 

10.44 

15.20 

36.57 

50.87 

26.29 

29.89 

4 

66.79 

74.19 

28.35 

27.88 

65.64 

68.53 

41.06 

40.16 

6 

79.14 

88.07 

51.65 

57.33 

83.74 

86.27 

64.64 

66.49 

8 

91.89 

97.40 

62.66 

76.11 

90.24 

95.60 

73.71 

81.53 

10 

97.15 

98.77 

71.65 

86.28 

94.08 

97.20 

86.09 

88.52 

12 

98.54 

99.30 

87.23 

93.44 

94.25 

97.53 

87.48 

92.68 

Cadmium  Manganese 


Days 

Nutrient  broth  Whey 

Ethanol 

Sucrose 

Nutrient  broth 

Whey 

Ethanol 

Sucrose 

(2:8) 

(1:9) 

(1:9) 

(1:9) 

(2:8) 

(1:9) 

(1:9) 

(1:9) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

33.19 

35.89 

30.30 

32.43 

38.81 

43.90 

34.20 

40.99 

4 

79.14 

84.18 

72.43 

81.13 

87.11 

91.90 

84.09 

91.32 

6 

98.54 

99.88 

97.99 

99.98 

98.75 

99.87 

98.32 

99.19 

8 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

10 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

12 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

Days 

Sulfate 

Nutrient  broth  (2:8) 

Whey  (1:9) 

Ethanol  (1:9) 

Sucrose  (1:9) 

0 

0 

0 

0 

0 

2 

17.78 

21.20 

15.09 

18.78 

4 

31.32 

39.98 

30.39 

37.18 

6 

46.63 

58.71 

46.90 

51.18 

8 

61.29 

78.88 

58.09 

68.81 

10 

70.37 

89.13 

70.87 

79.72 

12 

90.56 

97.12 

90.18 

93.87 
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Fig.  5    Reduction  in  iron  concentration  in  two-stage  anaerobic  bioreactor  using  whey  (1:9). 
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and  manganese)  and  sulfates  in  presence  of  nutrient 
broth  (2:8),  whey  (1:9),  ethanol  (1:9)  and  sucrose  (1:9) 
respectively.  After  12  days  of  biological  treatment, 
99.30%  and  97.53%  reduction  in  zinc  and  iron, 
respectively  was  observed  using  whey  as  electron 
donor  while  for  nutrient  broth,  ethanol  and  sucrose 
zinc  reduction  is  98.54%,  87.23%  and  93.44%, 
respectively  and  for  iron  it  is  94.25%,  87.48%  and 
92.68%,  respectively.  Sulfate  reduction  in  all  the 
electron  donor  sources  is  almost  90%  other  than  whey 
(97.12%).  After  6  days,  treatment  in  all  the  cases, 
cadmium  and  manganese  have  reached  to  a  level  in 
traces.  Whey  was  selected  as  best  source  for  electron 
donor  considering  the  reduction  in  heavy  metal 
concentration  and  its  cost. 

3.4  Two-Stage  Anaerobic  Bioreactors 

Figs.  4  and  5  show  the  treatment  pattern  of  two 
stage  anaerobic  bioreactors  for  zinc  and  iron  reduction. 
Upper  solid  lines  represent  first  stage  reactor  and 
lower  solid  lines  are  for  second  stage  reactor  treatment.  In 
both  first  stage  anaerobic  reactors,  after  4  days  of 
treatment  period  the  zinc  and  iron  concentration 
reduced  up  to  124.4  ppm  and  45.0  ppm,  respectively. 
5  L  from  total  10  L  was  transferred  to  second  stage 
treatment  process,  where  Zn&Fe  reduced  to  23.6  ppm 
and  8.5  ppm,  respectively,  while  the  first  stage  reactor 


which  has  been  leveled  up  with  fresh  effluent  reduced 
has  reduced  these  metals  to  67.4  ppm  and  24.4  ppm, 
respectively.  This  is  again  transferred  to  second  stage 
reactor  and  reduced  to  9.3  ppm  and  5.6  ppm.  The 
similar  process  is  being  continued  for  36  days  and 
under  steady  state  conditions  the  zinc  and  iron 
concentration  in  treated  effluent  was  4.9  ppm  and  1.8 
ppm  with  treatment  capacity  of  30  L  effluent  in  4  days. 
The  steady  state  condition  was  achieved  after  20  days 
of  treatment  period. 

4.  Conclusions  and  Recommendations 

Sulfate  reducing  bacteria  isolated  from  Acid  Plant 
of  Debari  Zinc  Smelter,  Hindustan  Zinc  Ltd.,  Udaipur, 
India  was  found  the  best  isolate. 

It  was  acclimatized  for  high  heavy  metals 
concentration  (zinc  500  ppm,  iron  ~  160  ppm) 
under  highly  acidic  solution  (pH  <  1.0). 

Whey  was  found  as  the  best  electron  donor  in  the 
effluent  treatment. 

Two-stage  anaerobic  bioreactor  with  treatment 
capacity  of  30  L  effluent  has  successfully  reduced 
heavy  metals  as  metal  sulfides  and  attained  almost 
neutral  pH  within  short  duration  of  4  days.  The  pilot 
runs  on  the  process  confirm  that  it  may  be 
commercialized  for  treatment  of  effluents  of  zinc 
plants. 
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Abstract:  To  investigate  the  influence  of  pressure  and  temperature  on  the  jet  velocity  of  a  three-dimensional  flow  was  the  main  goal 
of  this  study.  Using  a  precipitation  chamber  with  approximate  capacity  of  600  mL,  it  was  studied  the  thermodynamic  behavior  of 
supercritical  carbon  dioxide  mixture,  dichloromethane  and  grape  seed  extract  via  SAS  (supercritical  antisolvent  process).  For  the 
numerical  solution,  the  Navier-Stokes  equations  were  used  along  with  the  model  of  turbulence  k-e  and  Peng-Robinson  equation  of 
state  with  quadratic  mixing  rules  of  Van  der  Waals.  The  method  of  Chung  was  employed  to  determine  the  viscosity,  thermal 
conductivity  and  mass  diffusivity  of  the  flow  numerically  solved  through  commercial  code  based  on  CFD  (computational  fluid 
dynamics).  Simulations  for  pressures  between  80  bar  and  160  bar  and  temperatures  between  308.15  K  and  318.15  K  showed  large 
variations  in  the  jet  velocity,  an  important  property  in  the  dynamic  mixing  process  that  involves  the  size,  size  distribution  and  particle 
morphology. 

Key  words:  Supercritical  fluid,  CFD,  mathematical  modeling. 


Nomenclature 


k  Turbulent  kinetic  energy  (m2-s~2) 

p  Pressure  (bar) 

R  Universal  gas  constant  (J-moF'-K"  ) 

R  Radial  direction  (m) 

T  Temperature  (K) 

t  Time  (s) 

v  Velocity  (m-s"1) 

x  Axial  direction  (m) 

p  Density  (kg-m"3) 

a  Prandtl  number  (-) 

co  Fugacity  (-) 

A  Thermal  conductivity  (W-nf'-K"1) 

x  Tensor  stress  (pa) 

e  Rate  of  turbulent  kinetic  energy  dissipation  (m2-s"3) 

v  Molar  volume  (m^mol"1) 
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Ph.D.  Student,  research  fields:  computational  fluid  dynamics 
and  supercritical  technology.  E-mail:  flaviareitz@gmail.com. 


1.  Introduction 

Among  the  existing  methods  of  preparing 
nanoparticles,  SAS  (supercritical  antisolvent 
precipitation)  is  becoming  increasingly  important  due 
to  the  large  amount  of  chemicals  that  can  be  processed 
with  this  technology  [1-3].  Above  the  critical  point  of 
the  mixture,  there  condition  of  complete  miscibility, 
typically  in  the  range  of  8-16  MPa,  where  CO2  is 
completely  miscible  with  many  organic  solvents  and 
supercritical  turbulent  mixing  is  the  most  important 
issue.  The  SAS  process  involves  mixing  CO2  and 
organic  solution  (solute  dissolved  in  an  organic  solvent) 
from  two  outputs  of  a  coaxial  nozzle  system,  to  obtain 
the  desired  conditions  to  produce  thermodynamic 
supersaturation  and  particle  formation,  and  the 
resulting  PSD  (particle  size  distribution)  is  affected  by 
blending  [4].  While  blending  effects  are  well  known 
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to  those  based  on  precipitation  in  liquids  [5],  such 
effects  are  not  very  often  reported  in  supercritical 
systems. 

Despite  a  number  of  studies  that  support  this 
evidence  under  the  experimental  point  of  view,  little 
research  is  available  in  the  context  of  numerical 
modeling.  A  mathematical  model  for  mass  transfer 
between  a  droplet  of  organic  solvent  and  a  compressed 
antisolvent  was  developed  by  Werling  and 
Debenedetti  [6],  while  Chavez  et  al.  [7]  have  proposed 
a  methodology  to  calculate  the  time  scales  involved.  A 
more  elaborate  model,  solved  by  CFD  (computational 
fluid  dynamics)  was  proposed  by  Martin  and  Cocero  [3], 
in  which  a  turbulence  model  was  solved  in 
conjunction  with  mass  balances  and  took  into  account 
the  effect  of  mass  transfer  in  the  system.  The  authors 
employed  a  locking  to  solve  the  PBE  (population 
balance  equations)  [8].  The  results  were  strongly 
affected  by  the  value  of  solid-liquid  surface  tension. 
Forecasts  obtained  from  models  are  valuable  for 
understanding  the  process  behavior  at  a  large-scale 
industrial  level.  Henczka  et  al.  [9]  studied  numerically 
supercritical  antisolvent  precipitation  of  paracetamol 
with  carbon  dioxide  as  antisolvent  and  took  into 
account  the  heat  of  mixing  in  the  studied  system.  The 
PSD  was  calculated  using  the  standard  method  of 
moments,  and  the  interaction  between  mixing  and 
precipitation  was  taken  into  account  with  the  P-PDF 
approach  along  with  the  turbulent  model  of  Baldyga 
mixer  [10].  This  first  attempt  showed  the  importance 
of  including  mixing  effects  in  the  model,  and 
highlighted  some  limitations  of  the  approaches 
available.  Since  the  micromixing  models  available 
[1-8],  needed  to  explain  the  interactions  between 
turbulence  and  particle  formation  are  valid  only  for 
standard  fluids,  and  its  extension  to  supercritical  fluids 
is  not  simple  [11].  It  is  known  that  the  fluid  properties 
change  drastically  near  and  above  the  critical  point, 
then  the  dynamic  behavior  of  the  fluid  also  changes 
and  this  directly  influences  the  mixing  at  all  scales. 


A  large  number  of  publications  on  the 
micronization  of  different  materials  with  the  SAS 
process  appears  concerned  with  the  experimental 
analysis  of  the  effects  of  operating  conditions  on  the 
size  distribution  and  morphology  of  the  particles,  and 
how  this  type  of  analysis  is  often  difficult,  because 
the  interactions  between  these  parameters.  This  study 
aims  to  provide  a  mathematical  formulation  capable 
of  describing  the  influence  of  pressure  and  temperature 
on  the  jet  velocity  of  a  three-dimensional  flow 
phenomena  that  influence  the  turbulent  mixing 
observed  in  the  SAS  process.  Once  validated,  the 
effects  of  thermodynamic  properties  can  be  used  for 
modeling,  scale-up  and  optimization  of  systems,  with 
acceptable  precision  for  engineering  calculations. 

2.  Governing  Equations 

2.1  Thermodynamics 

The  density  of  the  mixture  is  described  by 
Peng-Robinson  equation  of  state  with  the  quadratic 
mixing  rules  of  Van  der  Waals: 

v-bm  v(v+bm)+bm(v-bm) 
where,  p  the  pressure,  v  molar  volume,  T  temperature, 
am  and  bm  the  constants  of  the  equation  for  the  case  of 
mixtures.  These  constants  are  obtained  by  mixing  rule 
is  given  by  am  =  a(T)  =  Yli'Zixixjaij  and  bm  = 
b(J)  =  TiiTi)  xixjaij  >  where,  x,  is  the  molar 
concentration  of  the  component,  i  is  combination 
rulesfraction:     ai;-  =  Ja^aj(l  —  and     b^  = 


This  formulation  allows  to  study  the  importance  of 
solving  the  resulting  system  of  equations  in 
non-isothermal  regime. 

2.2  Hydrodynamics 

The  conservation  equations  are  expressed  here  to  a 
differential  control  volume  fixed  in  three  dimensional 
space  of  the  flow  of  a  mixture  in  the  supercritical  state 
and  considering  the  turbulent  flow,  a  variable  can  be 
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considered  as  the  sum  of  an  average  value  plus  a 

fluctuation    0  =  0+  ^  =  0  +  0  .  The  equation  of 

continuity  in  Cartesian  coordinates  for  a  system  with 
variable  density  in  steady  state  is  given  by: 


di 


p  u, .   =  0 


(2) 


where,      with  i  e  {l,  2,  3},  represents  the  velocity 

vector  mean  component  u  to  the  position  vector 
components  respectively. 

The  equation  of  motion  that  takes  into  account  the 
density  and  viscosity  are  variable. 

f—  ~  ~  \     _  r~ 


di 


pui  Uj 


=  djP+dj 


Ti,  -pui  U- 


(3) 


where,  —  pu .  u  .  is  the  Reynolds  tensor,  resolved  by 
the  turbulence  model  and  the  two  equations 


( 

<  2 

dut 

dx. 

du, 
+ — - 

duk  > 

V 

J 

(4) 


Considering  the  viscous  dissipation  of  energy  Eq.  (4) 
and  species  conservation  Eq.  (5)  are  given  by: 


pun 


dx 
and 

d  ' 


d_ 

dx, 


■  + 


pUj  h 


J 


(5) 


dx . 


pujy^puj  yt 


d 


d)X  ; 


dx. 


(6) 


j  J 


where,  ^  is  the  enthalpy,  the  mass  fraction  yt  of 
component  i  and  Dm  the  molecular  diffusion 
coefficient. 

2.3  Model  of  Turbulence 

The  model  of  turbulence  k-s  for  two  equations  was 
used  in  this  study,  and  consists  of  two  semiempirical 
Eqs.  (7)  and  (8)  for  the  turbulent  kinetic  energy  k  and 
for  dissipation  rate  of  turbulent  kinetic  energy  s: 


d  (     ~  ^ 
pkuj 


dxj 


dXj 


<u  + 


Mi 


•kj 


dxj 


-k 


+  Pk -ps 

III  IV 

(7) 


df  ~} 

psuj 

V  J 


dx 


dx.. 


JU+ 


dx, 


(8) 

where,  e  is  the  dissipation  per  unit  mass.  In  Eqs.  (7) 
and  (8),  the  term  /  is  the  convective  transport  of  k  or  e; 
the  term  //  corresponds  to  diffusive  transport;  the  term 
///  is  the  production  rate;  and  the  term  IV  is  the 
destruction  rate.  Pu  is  a  term  of  turbulence  production 
due  to  viscous  forces. 

f     -  -T\      ~         /  „  \ 


Pk=jutVu: 


Vu+Vu 


Vu 


3juT.u+pk 


(9) 


where,  u  is  the  mean  velocity  vector.  The  empirical 


constants  of  the  model  are       =  0.09; 


Cn  =  1-44, 


C£2  =  1.92  ,  (jK  =  1  ,  oe  =  1.3  .  The  turbulent 
viscosity    can   be    calculated    with   the  equation 

_  cvPk2 
Ik  - 

This  model  of  two  equations  has  been  used  in  the 
description  of  turbulent  jet  under  supercritical 
conditions  [12]. 

3.  Methodology 

The  objective  of  this  study  is  to  investigate  the 
effects  of  temperature  and  pressure  in  jet  velocity  of  a 
three-dimensional  flow  in  SAS  process  and  the  work 
of  Franceschi  [13]  was  chosen  as  a  reference  to 
validate  the  model.  The  experimental  apparatus  consists 
of  a  tubular  precipitator  of  approximately  600  mL 
volume.  The  device  design  is  shown  in  Figs.  1  and  2. 
Table  1  summarizes  the  experimental  conditions 
investigated.  Under  these  conditions,  there  is  complete 
miscibility  between  CO2  and  dichloromethane  (single 
phase  region). 

The  simulations  were  made  from  mixed  mesh  with 
tetrahedral  elements  and  hexahedral  elements  and  a 
total  number  of  1.028  x  106  cells.  The  mesh  and 
geometry  details  associated  with  it  are  shown  in  Figs. 
1  and  2.  The  details  of  the  inner  cell  distribution  and 
refining  regions  at  the  exit  of  the  capillary  and  close  to 
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the  exit  of  the  mixing  tube  can  be  seen  in  Fig.  2  which 
shows  a  cutting  plane  at  the  center  of  the  chamber. 

The  chamber  has  an  internal  diameter  of  80  mm  and 
inside  length  of  120  mm,  capillary  length  coaxial 


r 


injection  of  98  mm  and  the  diameter  of  the  inlet  of 
CO2  is  0.97  mm  and  the  solution  inlet  diameter  is  0.1 
mm,  and  the  capillary  tube  starts  at  10  mm  inside  the 
chamber.  Contains  output  by  the  lid,  this  being  a  tube 


Inletsolution  and 

CO, 


Outlet 


Fig.  1    Mesh  with  detail  geometry  with  coaxial  input  and  C02  solution. 


Fig.  2  Aspects  of  the  internal  elements  of  the  mesh. 
Table  1    Operating  conditions. 


Simulation 

Temperature  (K) 

Pressure  (bar) 

1 

308.15 

80 

2 

313.15 

80 

3 

318.15 

80 

4 

308.15 

120 

5 

313.15 

120 

6 

318.15 

120 

7 

308.15 

160 

8 

313.15 

160 

9 

318.15 

160 

AAAAAAAAAAA** 

AAAAAAAAAAAA 


!  A  A  A  AAA  AAA AAA 

;  A  A  A  A  A  AAAAAAA 
AAAAAAA AAAA 

aaajjaaaaaaa. 

:i  AAAAAAAAAAAA 
* AAA AAA A A A A A A 
A A A A A A AAAAAAA 
A A A A A A AAA A AAA 
A A A A AAAAA A AAA 


'0 

j\A  .4 
AAM 
AAA 
AAA 
AAA 
A  A/ 

aa4 
AAA 
AA4 
A  *i 
*>A# 
SA* 


of  length  164  mm  and  internal  diameter  6.34  mm 
starting  at  35  mm  from  the  bottom  of  the  chamber  and 
the  end  of  the  outlet  tube  diameter  is  3.87  mm. 

To  solve  the  system  of  Eqs.  (l)-(8),  we  used  the 
ANSYS  FLUENT  14.5  software,  which  makes  use  of 
the  finite  volume  method.  The  segregated  solver  was 
employed  through  the  SIMPLEC  algorithm,  whose 
goal  is  to  solve  the  linear  systems  in  each  equation 
one  by  one,  updating  its  coefficients  due  to 
nonlinearities  and  the  coupling  between  variables. 
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As  convergence  criteria  was  considered  the  residual 
sum  below  10"5  for  all  variables  and  imposed  the  closure 


of  mass  balance 


™COl  +msol-mout 


<io-n 


4.  Results  and  Discussion 


The  results  obtained  from  simulations  performed  by 
the  software  for  carbon  dioxide  mixture, 
dichloromethane  and  the  grape  seed  extract,  since  the 
presence  of  the  grape  seed  extract  in  the  system  does 
not  change  the  phase  behavior  of  the  binary  solvent 
system  plus  antisolvent  [13]  so  that  still  considered 
the  system  in  a  single  phase. 

4. 1  Effect  of  Pressure 

Experimental  studies  have  shown  that  the  pressure 
directly  affects  the  size,  size  distribution  and  morphology 
of  the  particle,  mostly  related  to  temperature,  antisolvent 
flow  rate  and  the  concentration  of  the  solute  [13]. 

Analyzing  Figs.  3-5,  there  is  the  influence  on  the 
magnitude  of  the  velocity  obtained  by  the  centerline  of 
the  precipitation  chamber.  As  the  pressure  decreases, 
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the  jet  velocity  increases,  which  results  in 
homogeneous  and  nearly  spherical  particles,  as 
demonstrated  experimentally  by  Petit-gas  et  al.  [14] 
for  various  mixtures  in  the  SAS  process. 

Under  these  operating  conditions  and  according  to 
the  results,  the  most  advantageous  option  would  be  the 
use  of  lower  pressure,  as  experimentally  obtained  for 
Franceschi  [13]  in  a  precipitation  chamber. 

4.2  Effect  of  Temperature 

With  respect  to  the  influence  of  temperature  on 
the  velocity  of  the  jet  can  be  seen  in  Figs.  6-8  that 
their  maximum  values  promoted  increase  of  jet 
velocity  and  thus  increase  the  yield  of  precipitated 
particles,  as  described  by  Bashipour  and  Ghoreishi  [15] 
when  extract  P-carotene  from  the  aloe  vera  plant. 

In  addition  to  these  operating  conditions,  other 
factors  may  also  influence  the  velocity  of  the  jet  and 
hence  the  particle  size.  Cases  such  as  the  size  of  the 
precipitation  chamber,  flows  and  concentrations  of  the 
components  of  the  mixing  chamber  must  also  be  analyzed. 
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Fig.  3    Values  of  the  magnitude  of  the  velocity  along  the  axis  of  the  chamber  with  a  temperature  of  308.15  K. 
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Fig.  4    Values  of  the  magnitude  of  the  velocity  along  the  axis  of  the  chamber  with  a  temperature  of  313.15  K. 
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Fig.  5    Values  of  the  magnitude  of  the  velocity  along  the  axis  of  the  chamber  with  a  temperature  of  318.15  K. 
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Fig.  6    Values  of  the  magnitude  of  the  velocity  along  the  axis  of  the  chamber  with  a  pressure  of  80  bar. 
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Fig.  7    Values  of  the  magnitude  of  the  velocity  along  the  axis  of  the  chamber  with  a  pressure  of  120  bar. 
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Fig.  8    Values  of  the  magnitude  of  the  velocity  along  the  axis  of  the  chamber  with  a  pressure  of  160  bar. 


5.  Conclusions 

A  model  used  to  describe  the  velocity  of  the  jet  in 
the  precipitation  of  the  grape  seed  extract  employing 
the  SAS  process  was  developed.  The  model  includes 
the  density  of  the  mixture  via  Peng  and  Robinson  [  1 6] 
equation  of  state  with  quadratic  mixing  rules  of  Van 
der  Waals  for  the  analysis  of  the  mixture  behavior.  For 
the  investigated  variables,  pressure  and  temperature, 
results  showed  they  influence  the  velocity  of  the  jet 
solution  which  experimentally  impacts  the  size  and 
morphology  of  the  particles  with  a  stronger  effect  at 
lower  pressure  and  higher  temperature  operation. 

The  predictions  matched  well  with  the  experimental 
data  in  all  cases  studied,  and  although  they  have  not 
yet  combined  with  a  population  balance  that 
determines  the  size  and  size  distribution  of  particles, 
the  model  is  able  to  explain  the  experimental  trend  in 
the  jet  velocity  of  a  three-dimensional  flow  and  offers 
a  valuable  insight  into  the  process.  The  next  steps  of 
this  work  will  be  toward  another  validation  of  the 
model  in  a  wide  range  of  operating  conditions.  In 
particular,  the  cases  in  which  the  mixture  affects 


precipitation  will  be  investigated.  A  follow  up  study 
will  also  examine  the  range  of  Reynolds  numbers  for 
which  the  effects  of  turbulence  can  be  optimized, 
resulting  in  a  final  product  of  good  quality. 
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Abstract:  The  qualitative  and  quantitative  assessment  of  gas  flow  has  become  increasingly  relevant  in  the  use  of  everyday  systems. 
The  micro  flow  sensor,  developed  by  Innovative  Sensor  Technology  AG  (Switzerland),  is  by  principle  a  calorimetric  flow  sensor 
produced  as  a  micro  system  on  a  glass  substrate  by  means  of  photolithography  and  glass  etching  technology.  These  structures  are 
arranged  as  a  platinum  micro  heater  and  sensor  in  a  Wheatstone  bridge.  The  subsequent  etching  process  produces  an  exposed  area  of 
polyimide  membrane  that  is  only  a  few  microns  thick  and  includes  the  resistive  sensor  structure  as  the  active  area.  In  addition,  the 
RTD  (resistance  temperature  detector)  technology  included  on  the  sensor  allows  for  the  implementation  of  a  variety  of  electronic 
biasing  and  signal  processing  modes.  Since  the  sensor  can  be  powered  and  the  bridge  can  be  measured  in  both  CTA  (constant 
temperature  anemometer)  and  calorimetric  mode,  new  possibilities  are  presented  for  both  low  and  high  flow  rates  with  regard  to 
temperature  compensation,  self-calibration  and  self-monitoring. 

Key  words:  Calorimetric  flow  sensor,  pick  &  place  assembly,  CTA-mode,  low  and  high  flow  rates. 


1.  Introduction 

The  sensors  are  characterized  by  a  very  small  time 
response  and  low  power  consumption,  achieved  through 
the  principle  of  MEMs  (Micro-Electro-Mechanical 
Systems)  hot  film  anemometry.  Despite  an  extremely 
thin  membrane,  the  glass  carrier  produces  a  robust 
sensor  and  provides  easy  handling.  The  heater/sensor 
structure,  in  combination  with  the  full  bridge  circuitry 
"on  chip",  ensures  simple  signal  evaluation  and  offers 
the  ability  to  detect  the  flow  direction.  The  thin  film 
sputtering  and  etching  processes  of  1ST  AG 
(Innovative  Sensor  Technology  AG)  allow  for  low 
initial  development  costs  for  application-specific, 
custom  sensor  designs  and  support  cost  effective 
volume  production.  A  wide  range  of  sensor 
dimensions,  resistance  values,  structural  formations 
and  mounting  configurations  are  possible.  As  a  result, 
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cost  effective  solutions  including  a  sensor  on  PCB 
(printed  circuit  board)  with  electronic  signal 
processing  and  flow  channel  can  be  easily  developed 
in  a  timely  manner. 

2.  Sensor  Advantages 

Compared  with  its  silicon-based  counterparts,  the 
micro  flow  sensor  (Fig.  6)  from  1ST  AG  offers  major 
advantages  relating  to  dynamic  measuring  range  and 
assembly  of  the  sensing  element  into  customer 
systems  for  design-in  processes  (Fig.  1).  Due  to  the 
sensor's  ability  to  operate  in  a  combination  of  CTA 
(constant  temperature  anemometer)  and  calorimetric 
mode,  the  flow  signal  saturation  can  be  shifted  from  1  m/s 
(only  calorimetric  mode)  to  50  m/s  (both  modes 
combined).  These  flow  velocities  translate  to  a 
volumetric  flow  range  from  0.001  L/min  full  scale  to 
10  L/min  full  scale  (depending  on  the  flow  channel 
cross-sectional  area)  with  a  possible  resolution  of  1% 
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of  the  full  scale  value. 

The  response  time  of  10  ms  and  warm  up  time  of 
100  ms  allows  for  time -based  low  power  heater 
management.  With  the  platinum  technology  of  all 
resistive  structures,  a  heater  temperature  regulation 
loop  can  be  realized  in  order  to  maintain  a  constant 
temperature  and  to  protect  the  sensor  against  damage 
from  overheating.  Especially  during  power  down 
cycles,  the  resistive  structures  can  be  used  to  measure 
the  gas  temperature.  As  a  result,  it  is  possible  to 
maintain  active  temperature  compensation  by  tracking 
the  heating  power  of  the  sensor's  micro  heater. 

3.  Sensor  Assembly 

The  sensors  can  be  produced  economically  and  are 
well-suited  for  pick  &  place  assembly.  It  can  be 
handled  similar  to  an  IC-die  or  SMD  (surface  mount 
device)  component,  assembled  by  standard  pick  & 
place  machines  or  by  hand.  Electronic  connections  to 
a  PCB  can  be  achieved  through  wire  bonding  or  flip 
chip  technology.  The  sensor  can  be  assembled  over  a 
milled  PCB  channel  or  directly  embedded  into  a  PCB 
cavity  (Fig.  1).  The  goal  of  IST's  sensor/electronic 
interface  development  is  to  provide  customers  with  an 
easy-to-use  measurement  concept  for  realizing  their 
own  proprietary  solutions  down  to  the  component 
level  of  flow  measurement  technique. 

4.  Calorimeter  Mode 

The  most  common  used  measuring  mode  for  micro 
machined  flow  sensors  is  the  calorimetric  mode.  Here 
the  temperature  distribution  along  a 
sensor-heater-sensor-structure  depending  on  gas  flow 
is  used  to  provide  an  electrical  differential  signal  by 
thermal  imbalance  of  upstream  and  downstream 
sensor  elements  (e.g.  thermopiles)  (Fig.  2). 
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Fig.  2  Example  of  symmetrical  and  imbalanced 
temperature  distribution  on  a  3  mm  sensor  area. 

In  case  of  a  flow  rate  near  by  200  mL/min,  a 
temperature  difference  of  some  degree  centigrade  (°C) 
between  up-  and  down-stream  sensor  appears  and 
causes  a  bridge  signal  of  about  40-50  mV.  Here,  the 
user  is  faced  with  premature  saturation  effects, 
preventing  high  flow  applications.  This  can  be 
eliminated  by  using  both  flow  measuring  modes  at 
one  sensor  structure  (calorimetric  and  anemometric 
mode).  The  anemometric  mode,  especially  the 
CTA-mode,  is  responsible  for  compensating  the  loss 
of  heating  power  caused  by  gas  flow.  So  an  extended 
measuring  range  of  calorimetric  mode  can  be  achieved. 
And  additionally  the  output  of  the  CTA  can  be  used 
up  to  flow  rates  of  10  L/min  (Fig.  3). 

5.  CTA-Mode 

The  most  popular  anemometer  is  the  CTA.  It  is 
function  is  well  described  by  the  King's  Law: 

»  =  0.3-0.5  0-) 
By  converting  and  simplifying  this  equation  we 
obtain  this  formula: 


U  =  U0-  jl  +  k-v" 


(2) 


Fig.  1    Assembly  example  of  micro  flow  sensor. 


where,  U:  CTA-output;  U0:  free  convection  offset;  k: 
fluidic  depend  constant;  y  :  fluid  velocity. 

Thus,  the  output  is  determined  by  an  offset  (Uo),  an 
offset  dependent  slope  of  the  curve,  and  the  fluidic 
dependent  gain  value  (k). 

Here,  Uo  represents  the  value  of  constant 
temperature  difference  (AT)  between  the  heater  and 
fluid.  Strictly  speaking,  the  controller  of  a  CTA  keeps 
the  resistive  structure  at  a  constant  temperature.  So, 
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Fig.  3    Example  of  extended  flow  range  by  using  CTA-output. 

different  passivation  thicknesses  and  flow  element 
surfaces  impact  the  CTA  characteristics,  even  a 
deviation  in  the  heat  transfer  coefficient  ( a ).  The 
characteristics  depend  not  only  on  sensor 
alignment/orientation  and  fluid  type,  but  also  on  fluid 
temperature  and  sensor  contamination  (such  as  dust). 
So  the  equation  for  a  still  fluid  is: 

a0  -A-AT-U^l  +  k-v")  (3) 

a0-A-AT~  Ul  (4) 

That  means  each  deviation  of  CTA-characteristic, 
caused  by  free  convective  parameters  like  &o  can  be 
compensated  by  a  (70-adjustment  at  v  =  0  ■  In  order 
to  calibrate  the  A>value  (fluidic-caused  errors)  more 
than  one  point  is  necessary.  But,  fluidic-caused 
errors  mostly  depend  on  mechanical  alignment 
deviations  of  the  sensor  mounting.  If  mounting  errors 
can  be  determined  by  module  construction,  a  one 


point  calibration  provides  good  results  for  many 
applications. 

5.  Chip  Design  and  Circuit 

To  get  a  high  flow  rate  sensor  circuit  a  microflow 
sensor  has  to  be  implemented  into  a  CTA-bridge. 

The  complete  4-resistor  sensor  bridge  acts  as  the 
heater  of  the  CTA  (Rsl,  -,  Rs4).  The  outputs  "Flow 
Dir"  provide  once  the  calorimetric  mode  and  also  the 
signal  for  flow  direction. 

4-resistor  +  Temp,  sensor 


Pt  1200 


Using  as  full  bridge  in  a  CTA 
Flow  direction  sensitive 
Fig.    4       Four  resistors 
compensation. 


layout    with  temperature 
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Fig.  5    Sensor  layout  and  CTA  circuit. 
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Fig.  6    Microflow  sensor  element  3  mm  x  5  mm  x  0.5  mm. 
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Abstract:  A  three-layer  waveguide  structure  sensor  consists  of  LHMs  (left-handed  materials)  film  surrounded  by  dielectric  cladding 
and  antiferromagnetic  substrate  is  proposed.  LHMs  known  as  MTMs  (metamaterials)  have  simultaneous  negative  permeability  and 
permittivity.  The  dispersion  relation  for  the  structure  is  derived  for  TE  guided  modes.  Two  ranges  of  frequencies  are  chosen  such  that 
Voigt  permeability,  pv,  either  negative  or  positive.  The  sensitivity  is  proven  to  be  affected  by  different  parameters  including  the  film 
thickness,  LHM  parameters,  and  Voigt  frequency. 

Key  words:  MTMs,  LHMs,  antiferomagntic,  optical  waveguide  sensors,  sensitivity. 


1.  Introduction 

A  sensor  in  general  recognizes  and  senses  the 
presence  of  an  analyte  (the  specie  to  be  detected).  The 
ideal  sensor  not  only  has  to  respond  to  low 
concentrations  of  analytes,  but  also  must  have  the 
ability  to  discriminate  among  species.  Sensors  have  a 
broad  range  of  applications  in  medical,  pharmaceutical, 
environmental,  defense,  bioprocessing  or  food 
technology. 

The  need  for  fast,  easy  to  use,  compact  and  cheap 
sensors  is  the  motivation  of  this  study.  In  this  paper, 
the  proposed  structure  is  a  three-layer  system  consists 
of  new  man-made  materials  named  metamaterials  as 
film  sandwiched  between  an  antiferromagnetic 
substrate  and  dielectric  cladding.  Antiferromagnetic 
metal  is  studied  extensively  because  its  negative 
dielectric  constant  influences  the  dispersion  relation  in 
a  qualitative  manner. 

Ritchie  [1]  first  study  led  to  the  fast  growing  field 
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of  plasmonics  owing  to  their  potential  application  in 
areas  as  diverse  as  biosensing  [2].  Several  studies  is 
conducted  to  understand  both  surface  magnetic 
polaritron  (plasmons)  and  guided  waves  in  magnetic 
(anti/ferromagnetic)  layered  waveguides  [3-10]  where 
magnetic  layered  intercalated  with  dielectric  layers  or 
MTMs  (metamaterials)  layers. 

LHM  (left-handed  materials)  or  MTMs  are 
materials,  which  have  been  fabricated  in  the 
laboratory.  Veselago  [11]  theoretically  predicated  the 
existence  of  LHM  with  both  negative  permeability  (pi 
<  0)  and  negative  permittivity  (e  <  0)  and  verified 
experimentally  by  Pendry  et  al.  [12,  13]  and  Shelby  et 
al.  [14].  These  materials  have  unique  electromagnetic 
properties  made  them  attractive  to  several  applications 
including  antenna  design,  superlenses,  higher 
resolution  magnetic  resonance  imaging  [15,  16], 
sensors  [17-19]  and  isolators  [20,  21]. 

The  dispersion  relation  for  the  proposed  structured 
is  derived  for  TE  (transverse  electric  fields) 
propagating  modes  only.  The  homogenous  sensitivity, 
S,  is  defined  as  the  change  in  the  effective  refractive 
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index  with  cladding  index  [22],  is  then  calculated  to 
measure  the  ability  of  the  proposed  structure  to  detect 
and  sense  small  changes  at  the  cladding.  Next  section 
is  devoted  to  cover  the  background.  Numerical 
calculations  and  result  discussions  are  given  in 
Section  3.  The  paper  is  concluded  in  the  Section  4. 

2.  Theoretical  Background 

Fig.  1  displays  the  three-layered  waveguide  sensor. 
The  guiding  layer  consists  of  LHM  sandwiched 
between  a  semi-infinite  linear  dielectric  cladding  and 
a  semi-infinite  AFM  (antiferromagnetic)  substrate 
layer.  While  the  LHM  film  has  both  negative 
permittivity  s3  and  negative  permeability  fi3,  the  linear 
dielectric  cladding  has  positive  permittivity  si  and 
positive  permeability  /uj.  The  AFM  substrate  has 
positive  permittivity  S2  and  tensor  permeability  fi2 
defined  by  Eq.  (1). 

Mn      M2     0  ^ 
M02  ~ J  M12  M11 
0  0 


value  of  the  Voigt  permeability: 
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CO02  =  Y2H0,  coa  =  l2Ma,  com  =  4ny2Ms2,  co  is  the  circular 
frequency,  H0  is  the  external  magnetic  field  that  is 
applied  in  the  z  direction  and  perpendicular  to  the 
direction  of  propagation  of  the  waves,  Ms2  is  the 
unscreened  saturation  magnetization  in  the  direction 
of  Ho,  y2  is  the  absolute  value  of  the  gyromagnetic 
factor  in  the  antiferromagnet,  Ha  is  the  internal 
anisotropy  field  confined  in  the  x  direction,  and  He  is 
the  exchange  field. 

The  bulk  polaritons  for  antiferromagnet  in  an 
external  magnetic  have  two  frequency  gaps  [23,  24]. 
The  two  gaps  correspond  to  the  regions  of  negative 


(6) 


1  3 

where,  piv  has  two  poles  at  co  and  co  and  two  zeros  at 
co2  and  co4  given  by: 

(L3)2  2   ,      2  _/.     2     2,.  2,1  \2\"2 

of      =  ml  +  co02  +  con  ma  +1  4»02a»1  +  4a>m  maa02  +  (mm  coa  )  I 

(7) 

(2  4)       _  I    2  \1/2 

®    '     =+»02  +2a)ma)a)  (8) 

The  surface  polaritons  for  the  semi-infinite 
antiferromagnet  in  the  Voigt  configuration  have  two 
branches  for  each  direction  of  propagation.  The 
transverse  fields  TE  are  assumed  to  have  the 
oscillation  form  as  follows: 

E=(0   0   Ez)ej{°*-ky)  (9) 


H={HX    Hy  0)eJ 


(at  -ky ) 


(10) 


where,  k  =  Nco  I  c  is  the  propagation  constant,  TV  is 
the  effective  refractive  index,  and  c  is  the  speed  of 
light.  Applying  Eqs.  (9)  and  (10)  to  Maxwell's 
equations  and  applying  boundary  conditions  at  x  = 
-a/2  and  x  =  a/2,  we  obtain  the  dispersion  equation  as 
follows: 

iF3  (F2-Fl)cot(jBa)-(FlF2-F32)  =  0  (11) 

where,  F,  =  -jXh  F2  =  -j[i02  \Mi2~^2fin],  F3  =jl3,  Xj  = 
Nla,,  «/=(^2-«/)1/2,  n,=  {s,)m,  h  =  [Njun  ~ 
a2ni2\l[Nni2  —  azun],  a2=(N2  —  s2[iv)m,  fiv  =  (fin2  — 
Hi22)//un,  h  =  Nfi3/a3,  a3  =  (N2  —  s3/u3)m,  and  =ja3. 


x  ▲ 


x  -a/2 


x  =-al2 


Dielectric    e1>0,  fa  >  0 


LHM       £3  <  0,  n3  <  0 


-►y 


AFM      £2,  H2  (w) 


Fig.  1  The  proposed  three  layered  waveguide.  The  film  is 
LHM  with  £3  <  0,  fi3  <  0  surrounded  by  AFM  substrate 
and  dielectric  cladding. 
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The  homogenous  sensitivity  due  to  guided  waves,  S, 
is  calculated  from  the  dispersion  Eq.  (11)  by 
differentiating  the  effective  refractive  index  N  with 
respect  to  the  cladding  refractive  n ; 
be  expressed  as  Eq.  (12). 

dN 


S!m.  So,  S,  can 


S(a>)  = 


dn. 


(12) 


The  calculation  of  S  is  tedious;  hence,  it  is  not 
reproduced  here.  The  values  of  the  refractive  index  N 
for  surface  waves  are  first  numerically  calculated 
form  Eq.  (11)  and  accordingly  the  sensitivities  S  is 
calculated  numerically  from  Eq.  (12).  The  behavior 
of  our  system  will  depend  on  the  values  of  oo\  coi,  co2, 
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5" 

a  =  o.05 

3  - 

20  -10 

0 

-3  - 
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10            q  20 
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-30  -20  -1 
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Q  =  5 


and  co  .  Moreover,  it  will  be  affected  by  the  values 
of  MTM's  parameters  and  the  thickness  of  the 
film. 

3.  Results  and  Discussion 

In  this  work,  we  studied  the  sensitivity  of  the  sensor. 
Moreover,  we  test  the  effect  on  the  sensitivity  for  piv  < 
0  and  fiY  >  0.  In  the  numerical  calculation,  we  defined 
the  reduced  frequencies,  Q  =  colcom,  Q0  =  0002/00 /,  Qa  = 
co  J  coj,  Qm  =  cojcoj  and  wave  number,  Q  =  kc/com  = 
NQ. 

The  sensitivity  is  calculated  for  both  cases  when  /uv 
<  0  and  fxv  >  0. 


0  =  0.1 


Q=\5 


Q 


(c)  (d) 
Fig.  2    The  sensitivity  due  to  guided  waves  S  is  calculated  as  function  of  the  wave  number  Q  for  /#„  >  0  at  different  values  of 

reduced  frequency  Q  =  0.05,  0.1,  5  and  15.  The  film  thickness  a  =  0.05  cm. 
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Q  =  0.05 


10  -10 


0  =  0.1 


J  10 

Q 


Q=\5 


Fig.  3  The  sensitivity  due  to  guided  waves  S  is  calculated  as 
reduced  frequency  Q  =  0.05,  0.1,  5  and  15.  The  film  thickness, 

3. 1  Results  for  /uv>  0 

In  this  case,  the  reduced  frequency  is  tuned  such 
that  fiv  >  0.  The  sensitivity  S  is  plotted  in  Fig.  1  as 
function  of  wave  number  Q  at  the  chosen  set  of  values 
of  reduced  frequencies  Q  =  0.05,  0.1,  5  and  15.  The 
calculation  has  performed  using  the  following  values: 
a  =  0.05  cm,  Qm  =  5,  Qm-Qa  =  0.1,  n}  =1,  s2  =  1.0,  s3 
=  -1,  fi3  =  -4,  Qo  =  0.3,  fi02  =  2. 

From  Fig.  2,  we  can  see  that  S  varies  as  the  wave 
number  changes.  The  sensitivity  goes  to  zero  at  the 
zeros  of  wave  number  Q  which  corresponds  to  the 
zeros  of  co  and  infinities  at  the  poles  of  co  as  stated  in 
Eqs.  (7)  and  (8).  The  position  of  the  zeros  and 
infinities  changes  as  Q  changes.  Fig.  2  gives  us  a 


function  of  the  wave  number  Q  for  fiv  >  0  at  different  values  of 
a  =  0.15  cm. 

promising  sensor  for  the  stated  values.  The  system  is 
very  sensitive  with  changes  occurs  at  the  cladding 
layer.  In  this  case,  the  maximum  value,  S  =  ±7,  occurs 
at  Q  =  0.05  and  at  Q  =  ±20. 

To  study  the  effect  of  the  film  thickness,  the  same 
calculation  is  repeated  keeping  all  parameters 
unchanged  except  the  film  thickness  value  is  changed 
to  a  =  0.15  cm.  The  sensitivity,  S,  is  plotted  in  Fig.  3 
versus  the  wave  number,  Q,  at  Q  =  0.05,  0.1,  5  and  15. 

We  notice  that  the  maximum  sensitivity,  S  =  ±3, 
occurs  at  Q=  0.05  at  Q  =  ±3.  The  figure  exhibits  the 
zeros  and  poles  which  is  associated  with  co.  We  see 
the  effect  of  film  thickness  appears  at  the  limits  of  the 
propagation  constant,  Q,  which  decreases  dramatically 
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with  increasing  a.  Also,  the  maximum  value  of  S 
decreases  with  respect  to  the  previous  calculations. 

The  effect  of  changing  MTM's  parameters  are  done 
by  repeating  the  previous  calculations  keeping  all 
parameters  unchanged  and  swamping  the  values  of  S3 
and  n£  i.e.,  s3  =  -4,  pi3  =  -1.  Fig.  4  exhibits  the 
sensitivity,  S,  versus  the  wave  number,  Q,  at  a  =  0.05 
cm  and  £  =  0.05,  0.1,  5  and  15. 

The  maximum  sensitivity  S  =  ±100  occurs  at  Q  = 
0.05  and  Q  =  ±10.  From  previous  calculation,  we  can 
control  the  sensitivity  of  the  sensor  by  changing  the 
value  of  the  film  thickness  and  the  MTM's  parameters. 


(a) 


3.2  Results  for  piv<  0 

The  reduced  frequency  is  tuned  such  that  piv  <  0. 
The  sensitivity  S  is  displayed  in  Fig.  5  as  function  of 
wave  number  Q  at  chosen  set  values  of  reduced 
frequencies  Q  =  3.765,  3.7625,  6.8  and  6.95  that 
satisfied  the  condition  /uv  <  0.  The  calculation  has 
performed  using  the  following  values:  a  =  0.05  cm, 
Qm  =  5,  Qm-Qa  =  0.1,  ni  =  1,  s2  =  1.0,  s3  =  -1.0,/^  =  -4, 
Q0  =  03,fi02  =  2. 

The  sensitivity,  as  presented  in  Fig.  5,  has  its 
maximum  value,  S  =  7,  at  Q  =  6.7625  and  Q  =  10. 
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200  - 

fi=0.05 

~  1 
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,  ,  ,  , 

Q  10 

-200  - 

-400  - 

s 

4  ■ 


£=15 


(c)  (d) 
Fig.  4    The  sensitivity  S  is  calculated  as  function  of  the  wave  number  Q  for  fiv  >  0  at  different  values  of  reduced  frequency  £2 
=  0.05,  0.1,  5  and  15.  The  film  thickness,  a  =  0.05  cm  and  e3  =  -4.0,  fi3  =  -1. 
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£=6.7625 


Q  =  6.765 


(c)  (d) 
Fig.  5    The  sensitivity  S  is  calculated  as  function  of  the  wave  number  Q  for  pv  <  0  at  different  values  of  reduced  frequency  Q 

=  6.7625,  6.765,  6.8  and  6.95. 


4.  Conclusions 

The  homogenous  sensitivity  of  the  proposed  sensor 
is  studied.  The  sensor  is  a  three-layer  waveguide 
consists  of  a  LHM  film  surrounded  by  AF  substrate 
and  dielectric  cladding.  The  sensitivity  is  calculated 
for  range  of  frequencies  at  which  the  Voigt  frequency 
is  negative  and  positive.  Moreover,  we  studied  the 
effect  of  the  film  thickness  and  the  LHM  parameters. 
Results  show  that  the  sensitivity  of  the  proposed 
sensor  can  be  tuned  by  changing  the  value  of  the  film 
thickness  and  the  LHM's  parameters.  The  results 
might  be  of  interest  to  the  sensors  fabricators. 
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Abstract:  It  is  discovered  a  new  three-,  four-component  Petasis,  Passerini,  Hantzsch,  Kabachnic-Fields,  Ugi  reactions  with  of  arsine, 
stibine  and  bismuthine  in  organometallic  chemistry.  Modifications  were  replaced  to  a  nitrogen  atom  of  classical  reactions  of  atoms  of 
phosphorus,  arsenic,  antimony  and  bismuth.  It  has  been  proposed  a  new  mechanism  for  possible  reactions. 

Key  words:  Petasis,  Passerini,  Hantzsch,  Kabachnic-Fields,  Ugi  reactions,  phosphorus,  arsenic,  antimony,  bismuth. 

1.  Introduction 

In  recent  years,  it  has  been  actively  investigated 
nucleophilic  reaction  organilgalogenidami  phosphine 
and  with  electrophilic  geterilalkenami  allowing  to 
obtain  primary,  secondary  and  tertiary  phosphines, 
selectively  and  with  high  yield.  Most  of  these 
reactions  are  realized  under  mild  conditions  at  a 
pressure  atmlosfernom  phosphine. 

Although  the  chemistry  of  organophosphorus 
compounds  devoted  a  considerable  number  of 
monographs,  but  the  reaction  of  organic  compounds  of 
arsenic,  antimony  and  bismuth  are  very  limited.  The 
aim  of  this  work  searches  for  new  three-, 
four-component  reactions  Petasis,  Passerini,  Hantzsch, 
Kabachnic-Fields,  Ugi  reactions  with  organic 
compounds  of  phosphorus,  arsenic,  antimony  and 
bismuth. 

The  purpose  of  this  work,  search  for  new  three-, 
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four-component  reactions  Petasis,  Passerini,  Hantzsch, 
Kabachnic-Fields,  Ugi  reactions  with  organic 
compounds  of  phosphorus,  arsenic,  antimony  and 
bismuth. 

2.  Theory 

It  is  known  that  the  Petasis  reaction  is  a  MCR 
(multi  component  reaction)  that  enables  the 
preparation  of  amines  and  their  derivatives  such  as 
a-amino  acids. 

OH  0  R5  RVR= 

R'-%H+    If-V  +   r'-^h  "  R'"Vr 

R 

(1) 

OH  ^  I      QH  R«  RVR' 

0  R' 

(2) 

The  reaction  is  also  referred  to  as  the  Boronic  Acid 
Mannich  reaction,  since  it  proceeds  via  an  imine  with 
the  organic  ligand  of  the  boronic  acid  acting  as  the 
nucleophile,  similar  to  the  role  of  the  enolizable 
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ketone  component  in  the  original  Mannich  reaction. 
One    proposal    is    that    the    transfer    is  actually 


intramolecular,  and  takes  place  via  the  adduct  pictured 
below: 


HO,  OH 
HO-B- 


RVR" 


The  direct  reaction  with  glyoxylic  acid  merits 
particular  mention,  since  it  leads  to  interesting, 
unnatural  a-amino  acids  in  a  single  step,  while 
avoiding  the  appearance  of  toxic  byproducts  such  as 
seen  with  the  Strecker  Synthesis, 
o 

H ' 


OH 


R  OH 


OH 


RVR 

r'-^cooh 


(4) 


This  reaction  can  be  carried  out  with  secondary 
amines,  sterically  hindered  primary  amines, 
hydrazines  or  anilines  in  dichloromethane  at  room 
temperature.  The  range  of  potential  nucleophilic 
partners  includes  alkenylboronic  acids,  and 
electroneutral  as  well  as  electron-rich 
(hetero-)arylboronic  acids.  The  conversion  of 
electron-poor  boronic  acids  can  be  effected  at  elevated 
temperatures  MW  (microwave)  in  suitable  solvents. 

It  is  known  that  Passerini  reaction  three-component 
between  a  carboxylic  acid,  a  carbonyl  compound  such 
as  a  ketone  or  aldehyde,  and  an  isocyanide,  offers 
direct  access  to  a-hydroxy  carboxamides.  The 
Passerini  reaction  proceeds  rapidly  if  the  reaction  is 
performed  in  aprotic  solvents  at  room  temperature. 
High  yields  are  obtained  with  high  concentrations  of 
the  starting  materials  in  the  reaction  mixture.  From 
these  findings,  it  is  assumed  that  the  Passerini 
reaction  does  not  follow  an  ionic  pathway.  Hydrogen 
bonding  is  believed  to  play  a  crucial  role  in  the 
formation  of  the  presumed  cyclic  transition  state  for 
this  reaction. 

It  is  known  that  Hantzsch  Dihydropyridine 
(Pyridine)  Synthesis  allows  the  preparation  of 
dihydropyridine  derivatives  by  condensation  of  an 


-OH 

I  OH 
H 


(3) 


aldehyde  with  two  equivalents  of  a  P-ketoester  in  the 
presence  of  ammonia.  Subsequent  oxidation  (or 
dehydrogenation)  gives  pyridine-3,5-dicarboxylates, 
which  may  also  be  decarboxylated  to  yield  the 
corresponding  pyridines. 


0  0  \ 

R"^R"       R'"^"OH  I 

FT 


0  R"  R"  R"" 


R-CHO      +  2 


0 

0  0 


+  NH, 


OR" 


R"OOC 


COOR" 


(5) 

This  reaction  allows  the  preparation  of 
dihydropyridine  derivatives  by  condensation  of  an 
aldehyde  with  two  equivalents  of  a  P-ketoester  in  the 
presence  of  ammonia.  Subsequent  oxidation  (or 
dehydrogenation)  gives  pyridine-3,5-dicarboxylates, 
which  may  also  be  decarboxylated  to  yield  the 
corresponding  pyridines. 

It  is  known  that  Kabachnik-Fields  reaction 
three-component  coupling  of  a  carbonyl,  an  amine  and 
a  hydrophosphoryl  compound  leads  to 
a-aminophosphonates.  The  Kabachnik-Fields  reaction 
is  very  important  in  drug  discovery  research  for 
generating  peptidomimetic  compounds. 
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The  pathway  of  the  Kabachnik-Fields  reaction 
depends  on  the  nature  of  the  substrates.  The  amine 
and  hydrophosphoryl  compound  form  a  complex  in 
which  either  one  of  the  partners  may  react  with  the 
carbonyl  compound.  Often,  the  basicity  of  the  amine 
determines  the  reaction  pathway.  Weakly  basic 
amines  such  as  anilines,  which  can  act  as  proton 


donors,  favour  the  formation  of  an  imine,  whereas 
alkylamines  such  as  cyclohexylamines  do  not  form 
imines. 
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It  is  known  that  the  U-4CC  (Ugi  four-component 
condensation)  between  an  aldehyde,  an  amine,  a 
carboxylic  acid  and  an  isocyanide  allows  the  rapid 
preparation  of  a-aminoacyl  amide  derivatives.  The  Ugi 
reaction  products  can  exemplify  a  wide  variety  of 
substitution  patterns,  and  constitute  peptidomimetics 
that  have  potential  pharmaceutical  applications.  This 
reaction  is  thus  very  important  for  generating 
compound  libraries  for  screening  purposes. 
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The  mechanism  is  believed  to  involve  a  prior 
formation  of  an  imine  by  condensation  of  the  amine 
with  the  aldehyde,  followed  by  addition  of  the 
carboxylic  acid  oxygen  and  the  imino  carbon  across  the 
isocyanide  carbon;  the  resulting  acylated  isoamide 


rearranges  by  acyl  transfer  to  generate  the  final 
product. 

A  distinctive  feature  of  the  Ugi  reaction  is  the 
possibility  of  introducing  a  condensation  at  the  same 
time  four  or  more  components  under  very  mild 
conditions,  and  ease  of  implementation  and  high 
yields  of  the  reaction  products. 

A  special  activity  izonitrils  linked  due  to  the 
presence  of  divalent  carbon  atom  adjacent  to  the 
nitrogen  lone  pair  of  electrons  carrying  the  imparting 
of  nucleophilic  carbene  carbon  character. 

3.  Results  and  Discussion 

Comparison  of  the  three-,  four-component  Petasis, 
Passerini,  Hantzsch,  Kabachnic-Fields,  Ugi  reactions 
with  organic  compounds  of  arsenic,  antimony  and 
bismuth  showed  that  they  all  substitution  reactions  of 
amines  with  a  variety  of  organic  compounds  which 
have  a  labile  hydrogen  atom. 

We  proposed  a  general  view  of  the  proposed  new 
three-,  four-component  Petasis,  Passerini,  Hantzsch, 
Kabachnic-Fields,  Ugi  reactions  by  replacing  nitrogen 
atoms  arsenic,  antimony  and  bismuth. 

We  propose  to  extend  and  modify  Petasisa 
substitution  reaction  with  amines  on  the  nitrogen  atom 
of  phosphorus,  arsenic,  antimony  and  bismuth. 
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R1-E(OR)2  +  R2-C(=0)-R3  +  R4EHR5 
-ER4R5 


R1R2R3C 


R'-E(OR)2  +  R2-C(=0)-COOH  +  R3EHR4 
R'R2C(C00H)-ER3R4 


where,  E  =  N,  P,  As,  Sb,  Bi;  R1,  R2,  R3,  R4,  R5  =  H, 
CH3,  Ar. 

We  propose  modified  Hantzsch  condensation 
reaction  of  aldehyde  with  two  equivalents  of  a 
(3-ketoester  in  the  presence  of  ammonia,  phosphine, 
arsine,  stibine,  and  bismuthine.  After  oxidation  (or 
dehydrogenation)  gave  pyridin-3,5-dicarboxylic  acid, 
which  can  be  decarboxylated  to  give  the 
corresponding  pyridine,  phosphoridine,  arsinidine, 
stibinidine  and  vismutinidine. 

We  propose  to  expand  and  change 
Kabachnik-Fields  reactions  in  connection  with  the 
replacement  of  amines  on  the  phosphines,  arsines, 
stibines  and  vismutines. 

R-CO-R1    +    R2E!H2    +    H-E2(=0)(OR3)2  -> 
RR1  C(EHR3)-E(=0)(OR3)2 
where,  Eb  E2  =  N,  P,  As,  Sb,  Bi;  R,  R1,  R2,  R3 
Ar. 


Alk, 


We  propose  to  expand  the  Ugi  reaction  and  receive 
bimetallic  organic  compounds  of  arsenic,  antimony 
and  bismuth: 

REH2  +  R'R2C=0  +  E2HC=X  +  HOCOR4 
R-C(=0)-E2(R)-CH(R)-C(=0)-EHR. 
where,  X  =  O,  S,  Se;  E  =  P,  As,  Sb,  Bi;  R1,  R2,  R3,  R4 
=  Alk,  Ar. 

Thus,  we  first  proposed  to  extend  and  modify  the 
new    three-,    four-component    Petasis,  Passerini, 


(9) 

Hantzsch,  Kabachnic-Fields  Ugi  reactions  with  arsine, 
stibine  and  bismuthine  in  organometallic  chemistry. 
We  proposed  a  new  mechanism  for  the  possible 
reactions. 

4.  Conclusions 

The  data  presented  in  this  paper  show  that  the 
organic  compound  of  arsenic,  antimony  and  bismuth 
are  convenient  starting  materials  for  the  formation  of 
the  C-As,  C-Sb,  C-Bi  in  obtaining  various 
organometallic  compounds  which  are  more  widely 
used  as  catalysts,  phannaceuticals  preparations. 

We  have  opened  a  new  three-,  four-component 
Petasis,  Passerini,  Hantzsch,  Kabachnic-Fields,  Ugi 
reactions  with  organic  compounds  of  phosphorus, 
arsine,  stibine  and  bismuthine  in  organometallic 
chemistry.  We  proposed  a  new  mechanism  for  the 
possible  reactions. 
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